Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 126,947 G→T R418L (CGT→CTT)  aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,9470GT98.7% 256.6 / NA 79R418L (CGT→CTT) aceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base T (34/43);  minor base C (1/0);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAA  >  NC_000913/126807‑127088
                                                                                                                                            |                                                                                                                                             
tGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTggg                                                                                                                                       <  1:1006476/149‑1 (MQ=255)
 gCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGt                                                                                                                                      >  2:878863/1‑149 (MQ=255)
      agcagcTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCACGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAAACTCTGGTGCGAACCAGAGCCCTAACTGGGTAAtaa                                                                                                                                 >  2:665619/1‑147 (MQ=255)
       gcagcTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCGTCCAGAAATTCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGAt                                                                                                                                >  2:1332130/1‑149 (MQ=255)
       gcagcTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTCGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAAGGAt                                                                                                                                >  1:1263544/1‑149 (MQ=255)
        cagcTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATc                                                                                                                               <  1:1078146/149‑1 (MQ=255)
         agcTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATcc                                                                                                                              <  1:277073/149‑1 (MQ=255)
          gcTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAAGCTGAGCCTTAACTGGGTAATGATccc                                                                                                                             >  1:1114646/1‑149 (MQ=255)
           cTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCg                                                                                                                            <  2:506730/149‑1 (MQ=255)
            cccGGCAGGGACTGGCGGGGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGc                                                                                                                           <  1:430048/148‑1 (MQ=255)
                    ccaCTGGCGGTTGTATCCCTTGCATTCTGCCTTGGCCGAAGGTGGACTGCAGCAAGTTTTGTGAAATCGAATAATTGGAGCTGGGCCGCATCCATAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACt                                                                                                                   <  1:633956/147‑1 (MQ=255)
                    cGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATGCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACt                                                                                                                   >  2:919606/1‑149 (MQ=255)
                    cGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACt                                                                                                                   <  1:311463/149‑1 (MQ=255)
                     gACTGGCGGTGGTATGCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACta                                                                                                                  >  2:1205442/1‑148 (MQ=255)
                     gACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATTATCCCGCATGTTACTc                                                                                                                  <  1:880115/149‑1 (MQ=255)
                     gACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGCTCCCGCATGTTACTc                                                                                                                  >  1:550624/1‑149 (MQ=255)
                       cTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCAc                                                                                                                >  2:1283465/1‑149 (MQ=255)
                        tGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCCCt                                                                                                               >  1:44498/1‑149 (MQ=255)
                        tGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACt                                                                                                               >  1:1197507/1‑149 (MQ=255)
                          gCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTc                                                                                                             >  1:1329073/1‑149 (MQ=255)
                          gCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTc                                                                                                             >  1:600681/1‑149 (MQ=255)
                          gCGGTGGTATCCATGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTc                                                                                                             <  2:1006736/149‑1 (MQ=255)
                            ggtggtATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGa                                                                                                           <  1:149978/149‑1 (MQ=255)
                            ggtggtATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGa                                                                                                           <  2:561351/149‑1 (MQ=255)
                               ggtATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACaa                                                                                                        >  1:427300/1‑149 (MQ=255)
                                   tCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAAcc                                                                                                    <  1:973108/149‑1 (MQ=255)
                                    cccTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACACAAcca                                                                                                   >  2:407596/1‑148 (MQ=255)
                                     ccTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTCCTCACTTCGACAAAACCGa                                                                                                  >  1:1343985/1‑149 (MQ=255)
                                      ctcgCATGCGGCCGTGGCCGAAGGTGGACTTCAGCAAGTATGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTACTTATCCCGCATGTTACTCACTTCGACAAAACCGat                                                                                                 <  2:867908/146‑1 (MQ=255)
                                       tGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGata                                                                                                >  1:714612/1‑149 (MQ=255)
                                         gCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAACACCGATATc                                                                                              >  1:603347/1‑149 (MQ=255)
                                             gctgcCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCg                                                                                          <  1:665619/149‑1 (MQ=255)
                                               tgcCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAg                                                                                        >  2:320752/1‑149 (MQ=255)
                                               tgcCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAAAGGTCTCGCATGTTAGTCACTTCGACAAAACAGGTATCCCCGAg                                                                                        >  2:108028/1‑149 (MQ=255)
                                                 ccGTGGCCGACGCTGGACTTCAGCAAGAATGCGTCATCCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGCCAAAACCGATATCACCGAGtt                                                                                      >  2:303844/1‑149 (MQ=255)
                                                   ggggCCGAAGGTGGACTTAAGCAAGGTTGGTGAACTCGAAGAAGTGGAACTGGGCCGCATCCAGAAAACCTCTGGTGCGAACCGGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTgg                                                                                    <  1:1045494/147‑1 (MQ=255)
                                                    tGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGa                                                                                   <  1:673705/149‑1 (MQ=255)
                                                    tGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAAATGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGGGCCTTAACTGGGTAATGATCCCGCATGTTACTCACCTCGCCCCAACCGATCTCACCGAGTGGGa                                                                                   >  2:1355754/1‑149 (MQ=255)
                                                     ggCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGaa                                                                                  <  2:61218/149‑1 (MQ=255)
                                                     ggCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGaa                                                                                  <  1:1086660/149‑1 (MQ=255)
                                                      gCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAg                                                                                 <  2:138469/149‑1 (MQ=255)
                                                       ccGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGc                                                                                <  2:872325/149‑1 (MQ=255)
                                                        cGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCg                                                                               >  2:375039/1‑149 (MQ=255)
                                                        cGAAGGTGGACTTCAGCAAGTTCGCCGACCCCGCCGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTCACTCACTCCGCCCACACCGACACCACCGAGTTGGAACCg                                                                               >  2:992543/1‑149 (MQ=255)
                                                         gAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGCGTTgaggaacg                                                                              >  2:377161/1‑142 (MQ=255)
                                                          aaGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGtt                                                                             <  1:120599/149‑1 (MQ=255)
                                                           aGGTGGAGTTCATCAAGTTCGTTTAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGCGTAATGATCCCTCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTc                                                                            <  2:1009982/149‑1 (MQ=255)
                                                            ggTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTcc                                                                           <  2:833851/149‑1 (MQ=255)
                                                              tGGACTTCAGCAAGTTTGGTGAAGTCGAAGGTGTGGAACCGTGTCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATACCACCGAGTTGGAAGCGTTCCGt                                                                         <  1:33607/149‑1 (MQ=255)
                                                                 aCTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTaaa                                                                      <  1:1229645/149‑1 (MQ=255)
                                                                    tcagtcaGTTGGGTGAAATGGAAGACGTGGAGCTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGCTCCCGCATGTCACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAAcag                                                                   <  2:617425/143‑1 (MQ=255)
                                                                    tCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAAcag                                                                   <  1:1205442/149‑1 (MQ=255)
                                                                     cAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAAcagc                                                                  <  2:5858/149‑1 (MQ=255)
                                                                      gtcAAGTTTTGTGGAAACGGAGAAGTGGAACTGGGCCGCATCCAGAAAATCCCTGGTGCGAACCTTAGCCTTAACTGAGTAATGATCCCTCATTGTACCCACTTTGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAAcagca                                                                 <  2:653047/147‑1 (MQ=255)
                                                                         aaGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGaa                                                              <  2:1092174/149‑1 (MQ=255)
                                                                           gTTTGGTTAAATCGCAGACTTGGATCGGGGCCGCATACAGGAAACATCTGGTGCGATCTTGAGCCTTAAGTGGGTGATGGTCCCTCGTGTTACTCAGTGCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACg                                                            <  2:123106/149‑1 (MQ=255)
                                                                                gTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACgaagaa                                                       <  2:746373/149‑1 (MQ=255)
                                                                                   aaaTCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGGACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGCGTTGGAACCGTTCCGTAAACAGCAGAACGAAGCAgcg                                                    >  2:58162/1‑149 (MQ=255)
                                                                                    aaTCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAgcgg                                                   <  1:878863/149‑1 (MQ=255)
                                                                                         aagaagTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTCACTGGGGAATGATCCCGCACGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACCGCAGAACGAAGAAGAGGCGaaa                                              >  2:148119/1‑149 (MQ=255)
                                                                                              gTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCATTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTaa                                         >  1:1237385/1‑149 (MQ=255)
                                                                                                 gcaCTGGGCCGCATCCAGAAAACCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTCCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCt                                      <  1:1332130/147‑1 (MQ=255)
                                                                                                 gAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACACAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCt                                      >  2:25850/1‑149 (MQ=255)
                                                                                                         ccGCATCCAGAACATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACACCAGAACGCAAAAGCGGCGAAACGTAAGCTGGATGGGa                              >  2:1054431/1‑149 (MQ=255)
                                                                                                         ccGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGa                              <  1:1283465/149‑1 (MQ=255)
                                                                                                             aTCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGAt                          <  2:427300/149‑1 (MQ=255)
                                                                                                                  gAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGtt                                                                                      >  1:1293208/1‑84 (MQ=255)
                                                                                                                  gAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGtt                                                                                      <  2:1293208/84‑1 (MQ=255)
                                                                                                                   aaaaTCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCAcccc                    >  1:277674/1‑149 (MQ=255)
                                                                                                                    aaaTCTCTGGGGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCg                   <  1:159089/149‑1 (MQ=255)
                                                                                                                      aTCTCTGGTGCGAACCTGAGCCTTAACTGGGGAATTATCCCGCATGTTCCTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGt                 <  1:320752/149‑1 (MQ=255)
                                                                                                                           tGGCGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACGTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCt            <  1:1355754/149‑1 (MQ=255)
                                                                                                                              tGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTtca         <  2:166995/149‑1 (MQ=255)
                                                                                                                               gCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTtcat        <  2:1263544/149‑1 (MQ=255)
                                                                                                                               gCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGCTGGAAGCGTTCCGTAGACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACGCCGGTTGTCTtcat        >  1:1051006/1‑149 (MQ=255)
                                                                                                                                cgcaCCTGACCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCATCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTtcatc       <  2:1010771/146‑1 (MQ=255)
                                                                                                                                cGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGTAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGCGAACATCACCCCGGTTGCCTtcatc       >  2:537995/1‑149 (MQ=255)
                                                                                                                                     cTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATgga  >  1:897979/1‑147 (MQ=255)
                                                                                                                                     cTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGaa  >  2:667766/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                             
TGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTCATCATGAA  >  NC_000913/126807‑127088

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: