Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,381,333 A→T 5.9% L329Q (CTG→CAG)  frdA ← fumarate reductase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,381,3330AT 5.9% 323.4 / 7.3 101L329Q (CTG→CAG) frdAfumarate reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (45/50);  new base T (2/4);  total (47/54)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 102 total aligned reads displayed.

CCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTT  >  NC_000913/4381184‑4381481
                                                                                                                                                     |                                                                                                                                                    
ccACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCa                                                                                                                                                      >  1:1260635/1‑150 (MQ=255)
 cACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAg                                                                                                                                                     >  1:1177772/1‑150 (MQ=255)
    ggCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACg                                                                                                                                                  <  2:521312/150‑1 (MQ=255)
       gAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTc                                                                                                                                               >  2:895689/1‑150 (MQ=255)
        aaCAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCa                                                                                                                                              <  1:895689/150‑1 (MQ=255)
        aaCAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCa                                                                                                                                              >  2:544366/1‑150 (MQ=255)
        aaCAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCa                                                                                                                                              >  2:554302/1‑150 (MQ=255)
         aCAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCAt                                                                                                                                             <  1:1223821/150‑1 (MQ=255)
           aGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATTGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGTCTGCCGTTCAtct                                                                                                                                           <  2:778305/150‑3 (MQ=255)
           aGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGTCTGCCGTTCAtct                                                                                                                                           <  1:530125/150‑3 (MQ=255)
            gACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCa                                                                                                                                          >  2:808542/1‑150 (MQ=255)
            gACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCa                                                                                                                                          >  1:1571450/1‑150 (MQ=255)
            gACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCa                                                                                                                                          >  2:1019673/1‑150 (MQ=255)
            gACCTTTAAGGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCa                                                                                                                                          >  1:579567/1‑150 (MQ=255)
             aCCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGTCTGCCGTTCAtctgc                                                                                                                                         <  1:230495/150‑5 (MQ=255)
              ccTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGTCTGCCGTTCAtctgcg                                                                                                                                        <  1:579890/150‑6 (MQ=255)
                     tGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTtct                                                                                                                                 >  1:1445966/1‑150 (MQ=255)
                     tGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTtct                                                                                                                                 <  1:656669/150‑1 (MQ=255)
                       cGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCg                                                                                                                               <  2:1091788/150‑1 (MQ=255)
                         ggTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGcc                                                                                                                             <  1:1012384/150‑1 (MQ=255)
                            ttCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAg                                                                                                                          <  1:587851/150‑1 (MQ=255)
                                 aGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGAc                                                                                                                     <  2:1306995/150‑1 (MQ=255)
                                    tCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCa                                                                                                                  <  2:1376468/150‑1 (MQ=255)
                                    tCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCa                                                                                                                  <  1:990399/150‑1 (MQ=255)
                                    tCTGATCGGTTTCGATACCGCCCATGGTGTAGTATGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCACATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCa                                                                                                                  <  1:296817/150‑1 (MQ=255)
                                     cTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCaa                                                                                                                 >  2:458649/1‑150 (MQ=255)
                                      tGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAg                                                                                                                >  2:260428/1‑150 (MQ=255)
                                      tGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAg                                                                                                                >  1:494673/1‑150 (MQ=255)
                                      tGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCAGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAg                                                                                                                >  2:1121375/1‑150 (MQ=255)
                                        aTCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTc                                                                                                              <  1:465701/150‑1 (MQ=255)
                                        aTCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTc                                                                                                              <  1:51972/150‑1 (MQ=255)
                                             tttCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGAt                                                                                                         <  2:1015129/150‑1 (MQ=255)
                                             tttCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGAt                                                                                                         <  2:1338852/150‑1 (MQ=255)
                                                cGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATaaa                                                                                                      <  2:1608954/150‑1 (MQ=255)
                                                cGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATaaa                                                                                                      >  1:720366/1‑150 (MQ=255)
                                                cGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATaaa                                                                                                      <  2:274072/150‑1 (MQ=255)
                                                 gATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAAc                                                                                                     >  1:301058/1‑150 (MQ=255)
                                                  aTACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAAcc                                                                                                    <  1:458649/150‑1 (MQ=255)
                                                  aTACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGCTCCCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAAcc                                                                                                    <  2:723228/150‑1 (MQ=255)
                                                      cGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACAt                                                                                                >  2:930716/1‑150 (MQ=255)
                                                       gCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATc                                                                                               >  1:978675/1‑150 (MQ=255)
                                                         ccATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGc                                                                                             >  1:913714/1‑150 (MQ=255)
                                                           aTGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCa                                                                                           <  2:1445966/150‑1 (MQ=255)
                                                                     gtgCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGTGATAAACCACATCGCCACGCGGCGTgg                                                                                 <  1:182138/150‑1 (MQ=255)
                                                                     gtgCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTgg                                                                                 >  2:1046356/1‑150 (MQ=255)
                                                                     gtgCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTgg                                                                                 <  1:1438579/150‑1 (MQ=255)
                                                                        cGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGCCGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGgaga                                                                              >  1:1033560/1‑150 (MQ=255)
                                                                        cGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGgaga                                                                              >  1:1456622/1‑150 (MQ=255)
                                                                         ggTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGAt                                                                             >  2:1565614/1‑150 (MQ=255)
                                                                           tCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATgg                                                                           <  1:1051436/150‑1 (MQ=255)
                                                                           tCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATgg                                                                           <  2:1260635/150‑1 (MQ=255)
                                                                            cGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCTCCACGCCGCGTGGCGATgtg                                                                          <  2:720366/150‑3 (MQ=255)
                                                                            cGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCTCCACGCCGCGTGGCGATgtg                                                                          <  1:1507672/150‑3 (MQ=255)
                                                                              gACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCTCCACGCCGCGTGGCGATgtggt                                                                        <  2:189904/150‑5 (MQ=255)
                                                                                 gttccGGATTCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAAGGGCAGACGTTCATGCAGTTTTTTCTCGCGGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGGTGc                                                                     <  2:133792/147‑1 (MQ=255)
                                                                                 gTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGc                                                                     >  1:265660/1‑150 (MQ=255)
                                                                                  tACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGcc                                                                    <  1:1565614/150‑1 (MQ=255)
                                                                                    ccGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCtt                                                                  >  2:512659/1‑150 (MQ=255)
                                                                                    ccGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCtt                                                                  >  1:777710/1‑150 (MQ=255)
                                                                                      ggAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTa                                                                >  2:1134762/1‑150 (MQ=255)
                                                                                        aaTCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACg                                                              <  1:260428/150‑1 (MQ=255)
                                                                                            ggTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCAt                                                          >  2:140892/1‑150 (MQ=255)
                                                                                                cTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGt                                                      >  2:131855/1‑150 (MQ=255)
                                                                                                 tttAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTg                                                     <  1:1046356/150‑1 (MQ=255)
                                                                                                  ttAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGc                                                    <  2:530634/150‑1 (MQ=255)
                                                                                                  ttAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGc                                                    <  2:579567/150‑1 (MQ=255)
                                                                                                    aaCCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCa                                                  >  2:1256972/1‑150 (MQ=255)
                                                                                                        ggATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAg                                              <  1:554302/150‑1 (MQ=255)
                                                                                                        ggATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAg                                              <  1:544366/150‑1 (MQ=255)
                                                                                                         gATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTa                                                                <  2:630181/131‑1 (MQ=255)
                                                                                                         gATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTa                                                                >  1:630181/1‑131 (MQ=255)
                                                                                                             gACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCTTg                                         <  2:256348/150‑1 (MQ=255)
                                                                                                             gACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTg                                         <  2:1177772/150‑1 (MQ=255)
                                                                                                              aCGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTga                                        >  2:819210/1‑150 (MQ=255)
                                                                                                               cGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgag                                       <  2:71916/150‑1 (MQ=255)
                                                                                                               cGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgag                                       <  1:728128/150‑1 (MQ=255)
                                                                                                               cGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgag                                       >  1:1558043/1‑150 (MQ=255)
                                                                                                               cGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgag                                       <  1:808542/150‑1 (MQ=255)
                                                                                                                 ccAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgagag                                     >  1:1363599/1‑150 (MQ=255)
                                                                                                                    aCGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACt                                  >  2:772510/1‑150 (MQ=255)
                                                                                                                         cGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTc                             >  2:801296/1‑150 (MQ=255)
                                                                                                                           cTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTcgc                           >  1:794331/1‑150 (MQ=255)
                                                                                                                           cTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTcgc                           >  1:1354845/1‑150 (MQ=255)
                                                                                                                           cTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTcgc                           >  1:453308/1‑150 (MQ=255)
                                                                                                                            ttccGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTcgcg                          <  2:499704/147‑1 (MQ=255)
                                                                                                                                      tCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGt                <  1:1573255/150‑1 (MQ=255)
                                                                                                                                        gCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTc              >  1:1630871/1‑150 (MQ=255)
                                                                                                                                        ctgcaTGAACGTCTGCCGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTc              >  2:882930/5‑150 (MQ=255)
                                                                                                                                         cAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTcc             >  2:157791/1‑150 (MQ=255)
                                                                                                                                          gcaTGAACGTCTGCCGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCa            >  2:727882/3‑150 (MQ=255)
                                                                                                                                          aGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCa            <  2:494673/150‑1 (MQ=255)
                                                                                                                                            aTGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCata          <  2:777710/150‑1 (MQ=255)
                                                                                                                                            aTGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCata          <  2:1456622/150‑1 (MQ=255)
                                                                                                                                             tGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCatat         >  2:1606096/1‑150 (MQ=255)
                                                                                                                                              gAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCatata        >  1:1608985/1‑150 (MQ=255)
                                                                                                                                               aaCGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCatatat       >  2:1542858/1‑150 (MQ=255)
                                                                                                                                                aCGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATAtt      <  2:1571450/150‑1 (MQ=255)
                                                                                                                                                 cGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATAttt     >  1:34074/1‑150 (MQ=255)
                                                                                                                                                    cAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGtt  <  2:265660/150‑1 (MQ=255)
                                                                                                                                                    cAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGAACCAGTTCCATATATTTGtt  <  2:180439/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTT  >  NC_000913/4381184‑4381481

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