Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,381,335 A→C 6.9% R328R (CGT→CGG frdA ← fumarate reductase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,381,3350AC 6.9% 290.5 / 10.7 100R328R (CGT→CGGfrdAfumarate reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (43/50);  new base C (3/4);  total (46/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

GGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTT  >  NC_000913/4381188‑4381481
                                                                                                                                                   |                                                                                                                                                  
ggCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACg                                                                                                                                                  <  2:521312/150‑1 (MQ=255)
   gAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTc                                                                                                                                               >  2:895689/1‑150 (MQ=255)
    aaCAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCa                                                                                                                                              >  2:554302/1‑150 (MQ=255)
    aaCAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCa                                                                                                                                              >  2:544366/1‑150 (MQ=255)
    aaCAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCa                                                                                                                                              <  1:895689/150‑1 (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                                                                                                                                          >  2:808542/1‑150 (MQ=255)
        gACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCa                                                                                                                                          >  1:1571450/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:656669/150‑1 (MQ=255)
                 tGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTtct                                                                                                                                 >  1:1445966/1‑150 (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                                                                                                                >  1:494673/1‑150 (MQ=255)
                                  tGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAg                                                                                                                >  2:260428/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)
                                           tCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATaa                                                                                                       <  2:946240/150‑1 (MQ=255)
                                            cGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATaaa                                                                                                      <  2:1608954/150‑1 (MQ=255)
                                            cGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATaaa                                                                                                      <  2:274072/150‑1 (MQ=255)
                                            cGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATaaa                                                                                                      >  1:720366/1‑150 (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                                                                                 <  1:1438579/150‑1 (MQ=255)
                                                                 gtgCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTgg                                                                                 >  2:1046356/1‑150 (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)
                                                                                  ggAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTa                                                                >  2:1024683/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:544366/150‑1 (MQ=255)
                                                                                                    ggATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAg                                              <  1:554302/150‑1 (MQ=255)
                                                                                                     gATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTa                                                                >  1:630181/1‑131 (MQ=255)
                                                                                                     gATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTa                                                                <  2:630181/131‑1 (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:808542/150‑1 (MQ=255)
                                                                                                           cGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgag                                       >  1:1558043/1‑150 (MQ=255)
                                                                                                           cGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgag                                       <  2:71916/150‑1 (MQ=255)
                                                                                                           cGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgag                                       <  1:728128/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)
                                                                                                                                                   |                                                                                                                                                  
GGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTT  >  NC_000913/4381188‑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: