Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,933,049 A→G 5.4% V519A (GTG→GCG)  edd ← 6‑phosphogluconate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,933,0490AG 5.4% 229.7 / 2.4 75V519A (GTG→GCG) edd6‑phosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (30/40);  new base G (2/2);  total (32/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTCGCGAGCAGCCAGTTCCGCTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGA  >  NC_000913/1932907‑1933193
                                                                                                                                              |                                                                                                                                                
ttCGCGAGCAGCCAGGTCCGCTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTAc                                                                                                                                            <  4:295422/149‑1 (MQ=255)
 tCGCGAGCAGCCAGTTCCGCTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTAcc                                                                                                                                           <  5:18809/149‑1 (MQ=255)
  cgcgAGCAGCCAGTTCCGCGTCGTCTACCAGCAGCGTCAGTGCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTAcc                                                                                                                                           <  6:415131/148‑1 (MQ=255)
   gcgAGCAGCCAGTACCGCTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCt                                                                                                                                          <  8:397421/148‑1 (MQ=255)
     gAGCAGCCAGTTCCGCTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGaa                                                                                                                                       <  5:39948/149‑1 (MQ=255)
             aGTTCCGCTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGGAGCGCCGga                                                                                                                               <  6:325017/149‑1 (MQ=255)
             aGTTCCGCTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCg                                                                                                                                 >  2:121866/1‑147 (MQ=255)
              gTTCCGCTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTAc                                                                                                                                            <  7:168739/135‑1 (MQ=255)
              gTTCCGCTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTAc                                                                                                                                            >  8:168739/1‑135 (MQ=255)
               ttCCGCTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGGCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGgaga                                                                                                                             <  5:387232/149‑1 (MQ=255)
                 ccGCTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCGGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGgagcgg                                                                                                                           >  6:218179/1‑146 (MQ=255)
                   gCTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAgtcg                                                                                                                         >  1:78452/1‑149 (MQ=255)
                    cctCGGCTACCAGCAGCGGCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTACTGGGGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAgtcat                                                                                                                        <  7:340028/147‑3 (MQ=255)
                    cTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAgtcgt                                                                                                                        <  3:102741/149‑1 (MQ=255)
                    cTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAgtcgt                                                                                                                        <  2:78452/149‑1 (MQ=255)
                    cTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGGGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAgtcgg                                                                                                                        <  5:218179/149‑2 (MQ=255)
                        gTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCTATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCAt                                                                                                                    <  7:147319/149‑1 (MQ=255)
                           tACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCgg                                                                                                                 >  3:351860/1‑149 (MQ=255)
                           tACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCgg                                                                                                                 >  6:103197/1‑149 (MQ=255)
                             ccAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCg                                                                                                                  <  8:66493/146‑1 (MQ=255)
                             ccAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGt                                                                                                                >  2:212394/1‑148 (MQ=255)
                             ccAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGt                                                                                                                >  4:388516/1‑148 (MQ=255)
                                   gcGTCAGGTAGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGtaactaac                                                                                                         <  6:46820/149‑1 (MQ=255)
                                     gTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCGCTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACg                                                                                                        >  1:375784/1‑148 (MQ=255)
                                     gTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCGCTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACg                                                                                                        >  1:369841/1‑148 (MQ=255)
                                            cGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAAtttt                                                                                                 <  3:209940/148‑1 (MQ=255)
                                             gCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCATCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGTCACTTTACCTTAAGCGCCGGATAGTCGTCCATCGGTAACTAACGCAATTTTGa                                                                                               <  2:276291/149‑1 (MQ=255)
                                               ctgtctgtCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGaaa                                                                                             >  7:21085/1‑149 (MQ=255)
                                               ctgtctgtCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGaaa                                                                                             <  4:75664/149‑1 (MQ=255)
                                               ctgtctgtACATTCACTCGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGaaa                                                                                             <  1:174748/149‑1 (MQ=255)
                                                 gtctgtcCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAaca                                                                                           >  2:182025/1‑149 (MQ=255)
                                                            cacaCGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCaat                                                                                 <  3:388516/148‑1 (MQ=255)
                                                            cacaCGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAAAGCAATTTTTCAACACCGGTCCaat                                                                                 >  8:223721/1‑148 (MQ=255)
                                                             acacGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCaataa                                                                               <  7:376572/149‑1 (MQ=255)
                                                              cacGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCaataa                                                                               >  8:42128/1‑148 (MQ=255)
                                                               acGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGCGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAataata                                                                             >  5:209353/1‑149 (MQ=255)
                                                               acGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAataata                                                                             >  7:401892/1‑149 (MQ=255)
                                                                 gAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATaca                                                                           >  3:19005/1‑149 (MQ=255)
                                                                  aaTGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATacac                                                                          <  5:103197/149‑1 (MQ=255)
                                                                   atgctgTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACAcc                                                                         <  4:40338/149‑1 (MQ=255)
                                                                   atgatgTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGATACGTGGATAGCTGACGGCACTTTACCTGCAGCGCAGGAGAGTCGTCCCTCGGTTACTAACGCAATTTTGTATCCACCGTACAATAATacag                                                                          >  4:224269/1‑147 (MQ=255)
                                                                    tgatgTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACAcc                                                                         >  1:420391/1‑148 (MQ=255)
                                                                             cGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTggcggc                                                               <  1:212394/149‑1 (MQ=255)
                                                                              gTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATTATACCCCAAGTGGCGGCa                                                              >  6:142341/1‑149 (MQ=255)
                                                                                 gcgcACTTATGCCTCCAGCGCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGa                                                           <  3:315896/149‑1 (MQ=255)
                                                                                    cACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGtt                                                        <  7:223721/149‑1 (MQ=255)
                                                                                     aCTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGttt                                                       <  3:165437/149‑1 (MQ=255)
                                                                                     aCTTTTGCCAGCAGCCCGCCATCATAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGACGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGtt                                                        >  6:402093/1‑148 (MQ=255)
                                                                                      cTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAAt                                                                                                    <  3:194340/103‑1 (MQ=255)
                                                                                      cTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAAt                                                                                                    >  4:194340/1‑103 (MQ=255)
                                                                                      cTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTa                                                      <  2:420391/149‑1 (MQ=255)
                                                                                      cTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTa                                                      <  4:333679/149‑1 (MQ=255)
                                                                                       ttttGCCAGCATCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGAAGGCTCTTTACCTGAAGCGCCGGAGAGTCTTGCATATTTAACTAACGCAATTTTGAAACACCTGTCCAATAATACACCAAGTGGCGGCATGAGTTTAt                                                     <  7:84321/149‑1 (MQ=255)
                                                                                         ttGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTAt                                                     <  4:351860/147‑1 (MQ=255)
                                                                                           gCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTaa                                                 >  6:412681/1‑149 (MQ=255)
                                                                                            ccAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAgtcgt                                                                                                                        >  6:157771/1‑77 (MQ=255)
                                                                                            ccAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAgtcgt                                                                                                                        <  5:157771/77‑1 (MQ=255)
                                                                                             cagcagCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAAt                                                >  7:323903/1‑148 (MQ=255)
                                                                                               gcagcCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCt                                             <  8:323903/149‑1 (MQ=255)
                                                                                                  gcCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCGCTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGc                                          <  2:369841/149‑1 (MQ=255)
                                                                                                  gcCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCGCTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGc                                          <  2:375784/149‑1 (MQ=255)
                                                                                                          tCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCAGGAGTTTATGTAATTCTGGCATTCCGtt                                  <  2:331149/149‑1 (MQ=255)
                                                                                                                gCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCttt                            <  5:412681/149‑1 (MQ=255)
                                                                                                                 cTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCtttt                           <  5:244515/149‑1 (MQ=255)
                                                                                                                     tGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGcc                                                                                                                                  <  2:169079/42‑1 (MQ=255)
                                                                                                                     tGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGcc                                                                                                                                  >  1:169079/1‑42 (MQ=255)
                                                                                                                       gtgtTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGc                        <  1:121866/146‑1 (MQ=255)
                                                                                                                       gtgtTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAAATTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCATTTGGcccc                     >  6:289698/1‑149 (MQ=255)
                                                                                                                               tGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACg             <  1:5222/149‑1 (MQ=255)
                                                                                                                                  aTAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAACCACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGa           >  4:123732/1‑148 (MQ=255)
                                                                                                                                    aGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGacaa        >  7:7420/1‑149 (MQ=255)
                                                                                                                                      cTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACCCCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCGTTTGGCCCCTGATGACGGacaac       >  7:86005/1‑148 (MQ=255)
                                                                                                                                        gACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACaacaac    >  1:420900/1‑149 (MQ=255)
                                                                                                                                          cGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGa  >  2:132198/1‑149 (MQ=255)
                                                                                                                                          cGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGa  <  8:21085/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
TTCGCGAGCAGCCAGTTCCGCTTCGTCTACCAGCAGCGTCAGTTCGCCTGTCTGTCCATTCACACGAATGATGTCCCCGTCGCGCACTTTTGCCAGCAGCCCGCCATCGTAGGCTTCTGGTGTTACGTGGATAGCTGACGGCACTTTACCTGAAGCGCCGGAGAGTCGTCCATCGGTAACTAACGCAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGA  >  NC_000913/1932907‑1933193

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑