Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,381,331 G→T 5.9% P330T (CCG→ACG)  frdA ← fumarate reductase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,381,3310GT 5.9% 347.7 / 7.5 102P330T (CCG→ACG) frdAfumarate reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (48/48);  new base T (2/4);  total (50/52)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

Only 100 of 102 total aligned reads displayed.

CCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTT  >  NC_000913/4381183‑4381478
                                                                                                                                                    |                                                                                                                                                   
cccACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGc                                                                                                                                                    >  2:182138/1‑150 (MQ=255)
cccACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGc                                                                                                                                                    >  2:51972/1‑150 (MQ=255)
 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:554302/1‑150 (MQ=255)
         aaCAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCa                                                                                                                                           >  2:544366/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:1019673/1‑150 (MQ=255)
             gACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCa                                                                                                                                       >  1:1571450/1‑150 (MQ=255)
             gACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCa                                                                                                                                       >  2:808542/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                                                                                                               <  1:990399/150‑1 (MQ=255)
                                     tCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCa                                                                                                               <  2:1376468/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:1338852/150‑1 (MQ=255)
                                              tttCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGAt                                                                                                      <  2:1015129/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                                                                        <  2:1260635/150‑1 (MQ=255)
                                                                            tCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATgg                                                                        <  1:1051436/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:579567/150‑1 (MQ=255)
                                                                                                   ttAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGc                                                 <  2:530634/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                                    <  1:728128/150‑1 (MQ=255)
                                                                                                                cGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgag                                    <  1:808542/150‑1 (MQ=255)
                                                                                                                cGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgag                                    >  1:1558043/1‑150 (MQ=255)
                                                                                                                cGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgag                                    <  2:71916/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:1456622/150‑1 (MQ=255)
                                                                                                                                             aTGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCata       <  2:777710/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)
                                                                                                                                                    |                                                                                                                                                   
CCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTT  >  NC_000913/4381183‑4381478

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: