Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 622,455 G→C 5.2% W52C (TGG→TGC entS → enterobactin exporter EntS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913622,4550GC 5.2% 398.4 / 8.8 116W52C (TGG→TGCentSenterobactin exporter EntS
Reads supporting (aligned to +/- strand):  ref base G (53/57);  new base C (2/4);  total (55/61)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

Only 100 of 119 total aligned reads displayed.

AAACAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCTTA  >  NC_000913/622306‑622603
                                                                                                                                                     |                                                                                                                                                    
aaaCAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTgg                                                                                                                                                      >  2:268473/1‑150 (MQ=255)
aaaCAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTgg                                                                                                                                                      >  1:187654/1‑150 (MQ=255)
 aaCAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGc                                                                                                                                                     <  1:268473/150‑1 (MQ=255)
  aCAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCa                                                                                                                                                    >  1:741470/1‑150 (MQ=255)
   cAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAg                                                                                                                                                   >  1:1115595/1‑150 (MQ=255)
   cAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAg                                                                                                                                                   >  2:616746/1‑150 (MQ=255)
       ccTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTgg                                                                                                                                               <  1:610512/150‑1 (MQ=255)
       ccTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTgg                                                                                                                                               <  1:272199/150‑1 (MQ=255)
       ccTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTgg                                                                                                                                               <  2:422425/150‑1 (MQ=255)
         tGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTgggg                                                                                                                                             <  1:356038/150‑1 (MQ=255)
         tGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTgggg                                                                                                                                             >  2:1065370/1‑150 (MQ=255)
          ggCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaat                                                                                                                                            <  1:990027/150‑4 (MQ=255)
           gctgctTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaatg                                                                                                                                           <  1:862875/150‑5 (MQ=255)
            ctgctTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGCCAGGTggaatgt                                                                                                                                          <  1:676611/150‑6 (MQ=255)
            ctgctTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCAACTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCtt                                                                                                                                          >  1:692478/1‑150 (MQ=255)
                ttAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCgg                                                                                                                                      <  1:1398398/150‑1 (MQ=255)
                 tAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGt                                                                                                                                     >  1:1448945/1‑150 (MQ=255)
                    ccTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTg                                                                                                                                    <  2:751337/148‑1 (MQ=255)
                    ccTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGAc                                                                                                                                  <  1:1554772/150‑1 (MQ=255)
                    ccTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTACTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGAc                                                                                                                                  >  2:35440/1‑150 (MQ=255)
                     cTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACg                                                                                                                                 <  2:1448945/150‑1 (MQ=255)
                        aGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTg                                                                                                                              <  1:1246856/150‑1 (MQ=255)
                             gTTGAAAACGCACCCGGGGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGCCcgg                                                                                                                         >  1:862322/1‑150 (MQ=255)
                              ttGAAAACGCACCCGGCGTTTCGCGCATTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACcggc                                                                                                                        >  1:299990/1‑150 (MQ=255)
                               tGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACcggcg                                                                                                                       >  1:1441798/1‑150 (MQ=255)
                                     cGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGa                                                                                                                 <  2:1372200/150‑1 (MQ=255)
                                      gCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGGGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAt                                                                                                                <  1:1527972/150‑1 (MQ=255)
                                       cACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtg                                                                                                               <  1:1018819/150‑1 (MQ=255)
                                           cGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttg                                                                                                           <  1:870272/150‑1 (MQ=255)
                                            ggCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGAtgtttgt                                                                                                          >  2:627625/1‑150 (MQ=255)
                                              cGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTttgttg                                                                                                        <  2:1441798/150‑1 (MQ=255)
                                               gTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTgg                                                                                                       >  2:1076607/1‑150 (MQ=255)
                                               gTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTgg                                                                                                       >  1:504308/1‑150 (MQ=255)
                                               gTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTgg                                                                                                       >  2:1481628/1‑150 (MQ=255)
                                               gTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTgg                                                                                                       <  2:187654/150‑1 (MQ=255)
                                               gTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTgg                                                                                                       <  2:229797/150‑1 (MQ=255)
                                                 ttCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTACACCTGGGAGGTGGGGCTGTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGcc                                                                                                     >  2:940306/1‑150 (MQ=255)
                                                     cgcAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGAt                                                                                                 <  1:1511183/150‑1 (MQ=255)
                                                     cgcAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGGGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGAt                                                                                                 <  1:1505547/150‑1 (MQ=255)
                                                         gTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTc                                                                                             >  1:814227/1‑150 (MQ=255)
                                                              ccTCGGTCGTTTCATCTCAATTGTGGCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTcggcgg                                                                                        <  1:305385/150‑1 (MQ=255)
                                                                tcgctcgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGgtg                                                                                      <  2:370895/150‑1 (MQ=255)
                                                                 cgctcgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGCTCCAGATGATGACACATTCCACCTGGCAGGTGGCGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGgtgt                                                                                     >  2:1364280/1‑150 (MQ=255)
                                                                  gctcgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGtgtg                                                                                    <  1:71782/150‑1 (MQ=255)
                                                                  gctcgTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGtgtg                                                                                    <  2:1302601/150‑1 (MQ=255)
                                                                      gTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTgg                                                                                >  1:828187/1‑150 (MQ=255)
                                                                       tttCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGc                                                                               <  2:757915/150‑1 (MQ=255)
                                                                          cATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGa                                                                            <  2:1174416/150‑1 (MQ=255)
                                                                            tctcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATc                                                                          >  2:496608/1‑150 (MQ=255)
                                                                              tcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGc                                                                        >  1:774197/1‑150 (MQ=255)
                                                                              tcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGc                                                                        <  1:156710/150‑1 (MQ=255)
                                                                              tcAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGc                                                                        <  1:1190700/150‑1 (MQ=255)
                                                                               cAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCt                                                                       >  1:91445/1‑150 (MQ=255)
                                                                               cAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCt                                                                       >  2:1598298/1‑150 (MQ=255)
                                                                               cAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCt                                                                       >  2:31811/1‑150 (MQ=255)
                                                                                aaTTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTa                                                                      <  2:791553/150‑1 (MQ=255)
                                                                                  ttGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATg                                                                    <  1:527389/150‑1 (MQ=255)
                                                                                  ttGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATg                                                                    <  2:1389098/150‑1 (MQ=255)
                                                                                    gtgtCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAg                                                                  <  2:1115595/150‑1 (MQ=255)
                                                                                       tctcTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAgcgc                                                               >  1:34311/1‑150 (MQ=255)
                                                                                         tctGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCaa                                                             <  1:1466635/150‑1 (MQ=255)
                                                                                              gTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAg                                                        <  1:996041/150‑1 (MQ=255)
                                                                                               tttGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGt                                                       <  1:1314511/150‑1 (MQ=255)
                                                                                                  gCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAt                                                    >  1:1306923/1‑150 (MQ=255)
                                                                                                   cTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAg                                                        >  2:856554/1‑145 (MQ=255)
                                                                                                   cTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGAtt                                                   >  1:358104/1‑150 (MQ=255)
                                                                                                      ggCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTg                                                <  2:943654/150‑1 (MQ=255)
                                                                                                         gTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGTGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTg                                             <  1:1334229/150‑1 (MQ=255)
                                                                                                         gTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTg                                             >  1:346900/1‑150 (MQ=255)
                                                                                                           cgcgGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGGGGGCGGATCGCTATGAGCGCAAAAAAGTGATGTTGCTGgc                                           >  2:689635/1‑150 (MQ=255)
                                                                                                            gcgGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGgcg                                          <  1:951775/150‑1 (MQ=255)
                                                                                                                   cggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCa                                   <  1:1517169/150‑1 (MQ=255)
                                                                                                                    ggtgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCAc                                  >  1:555093/1‑150 (MQ=255)
                                                                                                                      tgcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCt                                >  1:484767/1‑150 (MQ=255)
                                                                                                                       gcAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCtg                               <  1:689635/150‑1 (MQ=255)
                                                                                                                        cAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCtgt                              >  2:1468117/1‑150 (MQ=255)
                                                                                                                          gATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTgg                            >  2:1149834/1‑150 (MQ=255)
                                                                                                                          gATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTgg                            >  1:400966/1‑150 (MQ=255)
                                                                                                                          gATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTgg                            <  1:998512/150‑1 (MQ=255)
                                                                                                                          gATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTgg                            >  1:1414016/1‑150 (MQ=255)
                                                                                                                            tCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCa                          >  1:768775/1‑150 (MQ=255)
                                                                                                                               aGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTg                       <  2:821717/150‑1 (MQ=255)
                                                                                                                               aGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTg                       <  1:1149834/150‑1 (MQ=255)
                                                                                                                                gatgatGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCAGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTgg                      >  1:1101385/1‑150 (MQ=255)
                                                                                                                                 atgatgACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGc                     >  1:1291691/1‑150 (MQ=255)
                                                                                                                                 atgatgACACATTCCACCTGGCAGGTGGGGCTTTCGGGGACGCGGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGGGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGc                     <  1:747721/150‑1 (MQ=255)
                                                                                                                                  tgatgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCt                    <  2:1312327/150‑1 (MQ=255)
                                                                                                                                  tgatgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCt                    >  1:88882/1‑150 (MQ=255)
                                                                                                                                  tgatgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTCTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTCGCt                    >  2:639749/1‑150 (MQ=255)
                                                                                                                                    atgaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTc                  >  2:1445699/1‑150 (MQ=255)
                                                                                                                                      gaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCAt                >  2:551578/1‑150 (MQ=255)
                                                                                                                                      gaCACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCAt                >  1:815913/1‑150 (MQ=255)
                                                                                                                                         acaTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTgg             <  2:814227/150‑1 (MQ=255)
                                                                                                                                          caTTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGa            >  2:775436/1‑150 (MQ=255)
                                                                                                                                          agccccACCTGCCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGa            >  1:1004716/5‑150 (MQ=255)
                                                                                                                                           gccccACCTGCCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGAc           >  2:755315/4‑150 (MQ=255)
                                                                                                                                              ccACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACtgt        >  2:433923/1‑150 (MQ=255)
                                                                                                                                               cACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACtgtg       >  1:995492/1‑150 (MQ=255)
                                                                                                                                                 ccTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCTGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGcc     <  2:1101385/150‑1 (MQ=255)
                                                                                                                                                  cTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCt    <  2:1617773/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAACAATCCTGGCTGCTTAACCTCAGCCTGTTGAAAACGCACCCGGCGTTTCGCGCAGTATTCCTCGCTCGTTTCATCTCAATTGTGTCTCTGGGTTTGCTCGGCGTCGCGGTGCCGGTGCAGATCCAGATGATGACACATTCCACCTGGCAGGTGGGGCTTTCGGTGACGCTGACCGGCGGCGCGATGTTTGTTGGCCTGATGGTCGGCGGTGTGCTGGCGGATCGCTATGAGCGCAAAAAAGTGATTTTGCTGGCGCGCGGCACCTGTGGCATTGGCTTCATTGGACTGTGCCTTA  >  NC_000913/622306‑622603

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: