Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 517,380 A→T 5.9% intergenic (+18/+45) 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,3800AT 5.9% 344.6 / 7.8 118intergenic (+18/+45)fetB/ybbNiron export ABC transporter permease; peroxide resistance protein/DnaK co‑chaperone, thioredoxin‑like protein
Reads supporting (aligned to +/- strand):  ref base A (66/45);  new base T (3/4);  total (69/49)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 120 total aligned reads displayed.

CAGATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCG  >  NC_000913/517234‑517527
                                                                                                                                                  |                                                                                                                                                   
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)
  gATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcatc                                                                                                                                                >  2:1734966/1‑150 (MQ=255)
      aTGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAAATCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcatcgcat                                                                                                                                            <  2:894747/150‑1 (MQ=255)
           gACCTTTATGCTGCTCTCAACCGCCAGCTTGGCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGGGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGATGCTGCGCATCCGGc                                                                                                                                       <  2:365697/150‑1 (MQ=255)
             ccTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGt                                                                                                                                     >  1:570307/1‑150 (MQ=255)
             ccTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGt                                                                                                                                     <  1:206385/150‑1 (MQ=255)
              cTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTTTCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGtt                                                                                                                                    <  2:512323/150‑1 (MQ=255)
              cTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGtt                                                                                                                                    <  1:1648053/150‑1 (MQ=255)
              cTTTATGCTGCTCTCAACCGCCAGCTTCGCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGtt                                                                                                                                    <  2:238415/150‑1 (MQ=255)
               tttATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGACGCCGGATGCGATGCTGCGCATCCGGcatcg                                                                                                                                   <  2:1544942/150‑5 (MQ=255)
               tttATGCTGCGCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTGTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGACGCCGGATGCGATGCTGCGCATCCGGcatcg                                                                                                                                   <  2:1537461/150‑5 (MQ=255)
                   tgctgcTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGt                                                                                                                               >  2:1303406/1‑150 (MQ=255)
                   tgctgcTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGt                                                                                                                               <  1:1819170/150‑1 (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:1371969/1‑150 (MQ=255)
                            cAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGcc                                                                                                                      <  2:1447821/150‑1 (MQ=255)
                              aCCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCTATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGa                                                                                                                    >  2:832909/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)
                                       ttGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGaa                                                                                                           >  1:1155800/1‑150 (MQ=255)
                                          tCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAatat                                                                                                        <  1:31316/150‑1 (MQ=255)
                                          tCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAatat                                                                                                        >  1:115431/1‑150 (MQ=255)
                                            gACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCTCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCa                                                                                                      >  1:370739/1‑150 (MQ=255)
                                              ccATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCGACGGAATATCAAt                                                                                                    <  2:437555/150‑1 (MQ=255)
                                              ccATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAAt                                                                                                    >  2:842153/1‑150 (MQ=255)
                                              ccATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAAt                                                                                                    >  1:1329652/1‑150 (MQ=255)
                                                ataattGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATac                                                                                                  <  2:579753/147‑1 (MQ=255)
                                                 ttattGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATaca                                                                                                 >  1:1410713/1‑150 (MQ=255)
                                                   attGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATacaac                                                                                               <  2:780827/150‑1 (MQ=255)
                                                   attGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATacaac                                                                                               >  1:23268/1‑150 (MQ=255)
                                                       ccTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATg                                                                                           <  2:606698/150‑1 (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)
                                                                  aCCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTgg                                                                                >  2:190742/1‑150 (MQ=255)
                                                                   ccTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGCTCAACGGAATATCAATACAACAATGCATACAGCTggc                                                                               >  1:1578155/1‑150 (MQ=255)
                                                                      aTCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTggcggc                                                                            >  2:68543/1‑150 (MQ=255)
                                                                      aTCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTggcggc                                                                            <  1:506764/150‑1 (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)
                                                                                    aaTTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGATGATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGcc                                                              >  1:1683365/1‑150 (MQ=255)
                                                                                    aaTTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAAAACAACAATGCATACAGCTGGCGGCGATACTTCGACGcc                                                              >  2:556942/1‑150 (MQ=255)
                                                                                        cgcgCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTg                                                          <  2:520913/150‑1 (MQ=255)
                                                                                        cgcgCCACCAGTTGGAGGTGACGCAATTGAAGAAGAACAGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTg                                                          <  1:702095/150‑1 (MQ=255)
                                                                                        cacgcCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCATAACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTg                                                          >  2:1694817/3‑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)
                                                                                                       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)
                                                                                                                 aTTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAg                                 <  1:385744/150‑1 (MQ=255)
                                                                                                                  ttGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAgg                                <  2:1750724/150‑1 (MQ=255)
                                                                                                                  ttGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAgg                                <  1:1371969/150‑1 (MQ=255)
                                                                                                                  ttGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAgg                                <  1:1303406/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)
                                                                                                                          gaaATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTg                        <  2:1738056/150‑1 (MQ=255)
                                                                                                                          caaaTGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTg                        <  1:1530083/149‑1 (MQ=255)
                                                                                                                              ggatgcgatgcCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAAAGGATTAACAATACAACAATGCATTCAGCTGGCGGCGATACTTCGACGCCAGAGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAAc                    >  1:1536141/2‑150 (MQ=255)
                                                                                                                               gatgcgatgcCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCACCGAGGATCTCCTGGAACg                   >  1:1782630/1‑150 (MQ=255)
                                                                                                                                atgcgatgcCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATAAATACAACACTGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGt                  >  2:1129159/1‑150 (MQ=255)
                                                                                                                                     atgcCGGCTGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGCTCAACGGAATATCAATACAACTATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGAGCTCCTGGAACGTTTTAc             <  2:1154129/150‑1 (MQ=255)
                                                                                                                                        ccGGATGCGATGCTGCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTAcgcg          >  1:964259/1‑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)
                                                                                                                                               cgcAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGAcc   >  1:131518/1‑150 (MQ=255)
                                                                                                                                                gcagcaTCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCg  <  2:128497/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAGATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCG  >  NC_000913/517234‑517527

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: