Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 517,383 A→T 6.2% intergenic (+21/+42) fetB → / ← ybbN iron export ABC transporter permease; peroxide resistance protein/DnaK co‑chaperone, thioredoxin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913517,3830AT 6.2% 340.2 / 9.9 112intergenic (+21/+42)fetB/ybbNiron export ABC transporter permease; peroxide resistance protein/DnaK co‑chaperone, thioredoxin‑like protein
Reads supporting (aligned to +/- strand):  ref base A (62/43);  new base T (3/4);  total (65/47)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

Only 100 of 126 total aligned reads displayed.

CAGATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCTGC  >  NC_000913/517234‑517532
                                                                                                                                                     |                                                                                                                                                     
cAGATTATGGTGCCCTTTATGCTGCTCTCACCCGCCAGCGTGTCGTCCATTATTGCCTGCTATTTAACCTATCGTCAGTTTTATAATTCGCGCCCCCAGTTGGTGGGGACGCAATTGAAGAAGAAAAGATGCGATGCCGGATGCgcagca                                                                                                                                                       <  2:116252/150‑1 (MQ=255)
 aGATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGAGGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcat                                                                                                                                                      >  2:927995/1‑150 (MQ=255)
 aGATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcat                                                                                                                                                      >  1:977193/1‑150 (MQ=255)
 aGATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcat                                                                                                                                                      >  2:506764/1‑150 (MQ=255)
 aGATTATGGTGACCTTTATGCAGCTCTCCACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGGGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcat                                                                                                                                                      <  1:503575/150‑1 (MQ=255)
  gATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcatc                                                                                                                                                     >  2:1734966/1‑150 (MQ=255)
         gTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCg                                                                                                                                              >  1:1354695/1‑150 (MQ=255)
           gACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGATGCTGCGCATCCGGc                                                                                                                                            <  1:1736385/150‑1 (MQ=255)
           gACCTTTATGCTGCTCTCAACCGCCAGCTTGGCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGGGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGATGCTGCGCATCCGGc                                                                                                                                            <  2:365697/150‑1 (MQ=255)
             ccTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGt                                                                                                                                          >  1:570307/1‑150 (MQ=255)
              cTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTTTCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGtt                                                                                                                                         <  2:512323/150‑1 (MQ=255)
              cTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGtt                                                                                                                                         <  1:1648053/150‑1 (MQ=255)
              cTTTATGCTGCTCTCAACCGCCAGCTTCGCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGtt                                                                                                                                         <  2:238415/150‑1 (MQ=255)
               tttATGCTGCGCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTGTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGACGCCGGATGCGATGCTGCGCATCCGGcatcg                                                                                                                                        <  2:1537461/150‑5 (MQ=255)
                   tgctgcTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGt                                                                                                                                    >  2:1303406/1‑150 (MQ=255)
                   tgctgcTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGGGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGt                                                                                                                                    <  1:369726/150‑1 (MQ=255)
                    gctgctCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATTCCGGTTGCGCAGCATCGCATCCGGCGTTGTgggt                                                                                                                                   >  2:765654/1‑148 (MQ=255)
                      tgctCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATATTTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGAGGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCa                                                                                                                                 <  2:426723/150‑1 (MQ=255)
                        ctctCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCAtg                                                                                                                               >  1:366571/1‑150 (MQ=255)
                          ctcAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAAATCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGAGCAGCATCGCATCCGGCGTTGAGGTTCAtgtg                                                                                                                             >  2:1698766/1‑150 (MQ=255)
                            cAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGcc                                                                                                                           <  2:1447821/150‑1 (MQ=255)
                              aCCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCTATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGa                                                                                                                         >  2:832909/1‑150 (MQ=255)
                              aCCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGa                                                                                                                         >  1:1163285/1‑150 (MQ=255)
                               ccGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGAt                                                                                                                        >  2:1409186/1‑150 (MQ=255)
                                    aGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACg                                                                                                                   <  1:710782/150‑1 (MQ=255)
                                    aGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACg                                                                                                                   <  1:1845798/150‑1 (MQ=255)
                                    aGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGATCATGTGCCGATCAACg                                                                                                                   <  2:463833/150‑1 (MQ=255)
                                    aGCTTGTCGACCATTATTGCCTGCTATTTAACATATCGTAAGTTTTATACTTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACg                                                                                                                   <  2:919556/150‑1 (MQ=255)
                                    aGCTTGTAGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACg                                                                                                                   <  1:277321/150‑1 (MQ=255)
                                     gCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACgg                                                                                                                  <  2:585009/150‑1 (MQ=255)
                                          tCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAatat                                                                                                             <  1:31316/150‑1 (MQ=255)
                                          tCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAatat                                                                                                             >  1:115431/1‑150 (MQ=255)
                                            gACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCTCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCa                                                                                                           >  1:370739/1‑150 (MQ=255)
                                              ccATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAAt                                                                                                         >  1:1329652/1‑150 (MQ=255)
                                              ccATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAAt                                                                                                         >  2:842153/1‑150 (MQ=255)
                                                ataattGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATac                                                                                                       <  2:579753/147‑1 (MQ=255)
                                                 ttattGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATaca                                                                                                      >  1:1410713/1‑150 (MQ=255)
                                                   attGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATacaac                                                                                                    >  1:23268/1‑150 (MQ=255)
                                                        cTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGc                                                                                               >  1:1823032/1‑150 (MQ=255)
                                                             aTTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAATAAATGATGCTATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGATTATCAATACAACAATGAATACa                                                                                          >  1:262339/1‑150 (MQ=255)
                                                             aTTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACa                                                                                          >  1:1632581/1‑150 (MQ=255)
                                                             aTTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACa                                                                                          >  2:1156121/1‑150 (MQ=255)
                                                                tAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCt                                                                                       >  2:1012511/1‑150 (MQ=255)
                                                                 aaCCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTg                                                                                      >  2:255683/1‑150 (MQ=255)
                                                                      aTCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTggcggc                                                                                 <  2:17237/150‑1 (MQ=255)
                                                                      aTCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTggcggc                                                                                 >  2:68543/1‑150 (MQ=255)
                                                                      aTCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTggcggc                                                                                 <  1:506764/150‑1 (MQ=255)
                                                                           aaGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATAc                                                                            >  1:749811/1‑150 (MQ=255)
                                                                            aGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACt                                                                           >  2:1648910/1‑150 (MQ=255)
                                                                            aGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATACGTCGTTGTGGTTCATGTGCCGAGCAAGGGAATATCAATTCAACAATGGAAACAGCTGGCGGCGATACt                                                                           >  2:319520/1‑150 (MQ=255)
                                                                             gTTTTATAATTCGCGCCACCAGTTGGTGGTGCCGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATCCtt                                                                          >  2:248982/1‑150 (MQ=255)
                                                                             gTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACtt                                                                          >  2:169557/1‑150 (MQ=255)
                                                                              ttttATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTc                                                                         >  2:694592/1‑150 (MQ=255)
                                                                               tttATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCg                                                                        >  2:1739631/1‑150 (MQ=255)
                                                                                   taATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGc                                                                    >  2:833910/1‑150 (MQ=255)
                                                                                    aaTTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGATGATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGcc                                                                   >  1:1683365/1‑150 (MQ=255)
                                                                                    aaTTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAAAACAACAATGCATACAGCTGGCGGCGATACTTCGACGcc                                                                   >  2:556942/1‑150 (MQ=255)
                                                                                          cgcCACCAGTTGGGGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGCTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCa                                                             <  2:1046895/150‑1 (MQ=255)
                                                                                            ccaccaGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATc                                                           >  2:1606947/1‑150 (MQ=255)
                                                                                            ccaccaGTTCGGGGGGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATc                                                           <  1:1629492/150‑1 (MQ=255)
                                                                                             caccaGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATcc                                                          >  2:1740825/1‑149 (MQ=255)
                                                                                             caccaGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCa                                                          >  2:16090/1‑150 (MQ=255)
                                                                                             caccaGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCa                                                          >  1:1762205/1‑150 (MQ=255)
                                                                                                caGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCAGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCAccc                                                       <  2:1306275/150‑1 (MQ=255)
                                                                                                     ggtggtGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTAcc                                                  <  1:1029733/150‑1 (MQ=255)
                                                                                                      gtggtgACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTAccc                                                 >  1:128497/1‑150 (MQ=255)
                                                                                                       tggtgACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTAcccc                                                >  2:47736/1‑149 (MQ=255)
                                                                                                       tggtgACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCca                                                >  2:154093/1‑150 (MQ=255)
                                                                                                        ggtgACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCcag                                               <  1:1149007/150‑1 (MQ=255)
                                                                                                          tgACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTTTGGTTCATTTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCTACGCCAGTGCCTCCCCCGGACCcagcg                                             >  2:1317270/1‑150 (MQ=255)
                                                                                                             cGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGCCGCCAGTGCATCACCCGTACCcagcgccg                                          >  1:922157/1‑148 (MQ=255)
                                                                                                             cGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCcagcgcag                                          >  1:56213/1‑150 (MQ=255)
                                                                                                                aaTTGAAGAAGAAATGATGCGATGCCGGATTCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATAAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGa                                       <  1:1272037/150‑1 (MQ=255)
                                                                                                                 aTTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAg                                      <  2:204476/150‑1 (MQ=255)
                                                                                                                  ttGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAgg                                     <  1:1303406/150‑1 (MQ=255)
                                                                                                                  ttGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAgg                                     <  2:1750724/150‑1 (MQ=255)
                                                                                                                  ttGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAgg                                     <  1:1371969/150‑1 (MQ=255)
                                                                                                                     aagaagaaATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGAtc                                  <  2:754381/150‑1 (MQ=255)
                                                                                                                      agaagaaATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGAtct                                 >  1:1094173/1‑150 (MQ=255)
                                                                                                                       gaagaaATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGAtctc                                >  2:1052291/1‑150 (MQ=255)
                                                                                                                        aagaaATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTcc                               >  1:332747/1‑150 (MQ=255)
                                                                                                                          caaaTGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTg                             <  1:1530083/149‑1 (MQ=255)
                                                                                                                              tgatgCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAAc                         >  2:1832382/1‑150 (MQ=255)
                                                                                                                              ggatgcgatgcCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAAAGGATTAACAATACAACAATGCATTCAGCTGGCGGCGATACTTCGACGCCAGAGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAAc                         >  1:1536141/2‑150 (MQ=255)
                                                                                                                               gatgcgatgcCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCACCGAGGATCTCCTGGAACg                        >  1:1782630/1‑150 (MQ=255)
                                                                                                                                     atgcCGGCTGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGCTCAACGGAATATCAATACAACTATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGAGCTCCTGGAACGTTTTAc                  <  2:1154129/150‑1 (MQ=255)
                                                                                                                                     acgcCGGATGCGATGCTGCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTAc                  >  1:1872705/3‑150 (MQ=255)
                                                                                                                                     acgcCGGATGCGATGCTGCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTAc                  >  1:1571955/3‑150 (MQ=255)
                                                                                                                                         cGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGt              >  2:625909/1‑150 (MQ=255)
                                                                                                                                           gATGCGCAGCATCGCATCCGGCGTTGTGGCTCATGTGCCGATCAACGGAATATCAATACAACAATGCACACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCt            <  2:949418/150‑1 (MQ=255)
                                                                                                                                            aTGCGCAGCATCGCATCCTGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTg           <  1:927995/150‑1 (MQ=255)
                                                                                                                                             tGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGa          >  1:1377832/1‑150 (MQ=255)
                                                                                                                                              gcgcAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGAc         <  2:693692/150‑1 (MQ=255)
                                                                                                                                               cgcAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGCCGCCAGTGCATCCCCCGTACCCAGCGCAGCGAGGCTCTCCTGGAACGTTTTACGCGTCTGAcc        >  2:818845/1‑150 (MQ=255)
                                                                                                                                               cgcAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGAcc        <  2:1557081/150‑1 (MQ=255)
                                                                                                                                                gcagcaTCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCg       <  2:128497/150‑1 (MQ=255)
                                                                                                                                                     atcgcatcCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCTgc  <  2:779535/150‑1 (MQ=255)
                                                                                                                                                     atcgcatcCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCTgc  >  2:583810/1‑150 (MQ=255)
                                                                                                                                                     atcgcatcCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCTgc  <  1:1636734/150‑1 (MQ=255)
                                                                                                                                                     atcgcatcCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGGTTTACGCGTCTGACCGTCTgc  >  2:725957/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCTGC  >  NC_000913/517234‑517532

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