Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,906,229 T→G 8.3% I181S (ATC→AGC) ‡ mntP → putative Mn(2+) efflux pump, mntR‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,906,2290TG 8.3% 259.4 / 14.1 97I181S (ATC→AGC) ‡mntPputative Mn(2+) efflux pump, mntR‑regulated
Reads supporting (aligned to +/- strand):  ref base T (48/40);  new base G (5/3);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

ACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATG  >  NC_000913/1906082‑1906375
                                                                                                                                                   |                                                                                                                                                  
aCATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATcc                                                                                                                                                  >  2:615842/1‑150 (MQ=255)
 cATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCct                                                                                                                                                 >  2:217138/1‑150 (MQ=255)
      tCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGa                                                                                                                                            <  2:551802/150‑1 (MQ=255)
         cgACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGc                                                                                                                                         >  1:1657696/1‑150 (MQ=255)
          gACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAGGATCTGGACGcc                                                                                                                                        <  1:996586/150‑2 (MQ=255)
           aCCGCATTGGCCATCGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCAc                                                                                                                                       <  1:94245/150‑1 (MQ=255)
            ccGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGTGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAGGATCTGGACGccga                                                                                                                                      <  2:1473317/150‑4 (MQ=255)
            ccGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACt                                                                                                                                      >  2:459107/1‑150 (MQ=255)
               cATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTcc                                                                                                                                   >  2:1400618/1‑150 (MQ=255)
               cATTGGCCATTGGTTGTGCAACCTTGATTATGCCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTcc                                                                                                                                   >  2:1476094/1‑150 (MQ=255)
                aTTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCTCTTCCa                                                                                                                                  <  2:5809/150‑1 (MQ=255)
                aTTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGAGGATCTGGACGccgatgcc                                                                                                                                  <  2:52279/150‑8 (MQ=255)
                   ggCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACgg                                                                                                                               >  1:234485/1‑150 (MQ=255)
                     ccATTGGTTGTGCAACCTTTATTATGTATACATTAGGGATGATGGTTTGTCGCTTTATCGGCTCTATTATTTTTTATATAGCTTAATTTTTCGTCGGGCTGGTGCTGATCGGGATATGCTTCAATATCATCTGGACGCACTTCCACGGtt                                                                                                                             <  2:1146509/150‑1 (MQ=255)
                     ccATTGGTGGTGCACCCTTGATTATGTCAACATTAGGGATGATGGGTGGGCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGtt                                                                                                                             <  1:769171/150‑1 (MQ=255)
                         tGGTTGTGCAACCTTGATTATGTCAAGATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCTTCCAGATCCTCTGGACGCACTTCCACGGTtaata                                                                                                                         <  1:673575/150‑1 (MQ=255)
                         tGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaata                                                                                                                         >  1:1110343/1‑150 (MQ=255)
                          ggTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaataa                                                                                                                        >  1:800295/1‑150 (MQ=255)
                          ggTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaataa                                                                                                                        >  2:1069697/1‑150 (MQ=255)
                          ggTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaataa                                                                                                                        >  2:620423/1‑150 (MQ=255)
                          ggTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTtaataa                                                                                                                        >  2:1339896/1‑150 (MQ=255)
                           gTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAg                                                                                                                       >  1:734075/1‑150 (MQ=255)
                            ttGTGCAACCTTGATTATGTAAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGTAAAAATGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGa                                                                                                                      <  2:664779/150‑1 (MQ=255)
                               tgCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACg                                                                                                                   <  2:1202061/150‑1 (MQ=255)
                                 cAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGGGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCt                                                                                                                 <  1:1046505/150‑1 (MQ=255)
                                 cAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATAGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCt                                                                                                                 >  2:1594118/1‑150 (MQ=255)
                                  aaCCTTGATTATGTCAACATTCGGGATGTTGGTGGGTCGCTTAAGCGGCTCAATTATTGGTAAAAACGCGGAAATTCTCGGCGGGCTGGGGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTg                                                                                                                <  2:589802/150‑1 (MQ=255)
                                  aaCCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTg                                                                                                                >  2:401454/1‑150 (MQ=255)
                                    ccTTTATTTTGTCAAAATTAGTTTTGATGGTTGGGCTTTTTTTCGTATAATTTTTTGGGAAAACAGCTTAAATTCTCGCGGGGCGGGTGCGTATCGGCATCGGGGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGcc                                                                                                              <  2:1292680/150‑1 (MQ=255)
                                    ccTTGTTTATGGCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGGGCGGATCGGCAGCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGcc                                                                                                              <  2:357729/150‑1 (MQ=255)
                                    ccTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGcc                                                                                                              >  1:354456/1‑150 (MQ=255)
                                     cTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTTTTTGGAAAAATGCGGAAATTCTCGGCGGGCTGGTGATGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATATTAACGCTGCCa                                                                                                             >  2:1642657/1‑150 (MQ=255)
                                        gATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCagag                                                                                                          <  1:218728/150‑1 (MQ=255)
                                         aTTATGTCTACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGTGAAAAACGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATGGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAgg                                                                                                         <  2:278189/150‑1 (MQ=255)
                                          ttttGTAAACTTTAGGGTTTAGGGTTGGGCGCTTATTCGGCAAAATTTTGGGTAAAAAAGCGGAAATTCTCGGCGGGCGGGTGCTGATCGGCATCGGCGCCCAGATCCTCTGTACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGt                                                                                                        <  2:353433/147‑1 (MQ=255)
                                           tATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTg                                                                                                       >  2:1043811/1‑150 (MQ=255)
                                           tATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTg                                                                                                       >  2:1313962/1‑150 (MQ=255)
                                           tATGTCAACATTAGGGATGAGGGTGGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTg                                                                                                       <  2:824721/150‑1 (MQ=255)
                                             tGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGaa                                                                                                     >  2:1637119/1‑150 (MQ=255)
                                               tCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAata                                                                                                   >  2:1570811/1‑150 (MQ=255)
                                                cAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAatat                                                                                                  >  1:1694955/1‑150 (MQ=255)
                                                aaaCATTCGGGATTTTGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCTGGGTGGTGCTGATCGGCATCGGCGTCCATATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAatat                                                                                                  <  2:1687776/149‑1 (MQ=255)
                                                     ttcgggATGATGGTTGGGCGATTTATCGGCTCATTTCTTGGGTAAAAAGGGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTGCTGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAg                                                                                             <  1:339051/147‑1 (MQ=255)
                                                     ttAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAg                                                                                             <  2:783345/150‑1 (MQ=255)
                                                         ggATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCg                                                                                         >  2:408805/1‑150 (MQ=255)
                                                           atgatgGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTc                                                                                       <  2:1020325/150‑1 (MQ=255)
                                                           atgatgGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTc                                                                                       >  1:819116/1‑150 (MQ=255)
                                                            tgatgGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCt                                                                                      <  2:1110343/150‑1 (MQ=255)
                                                            tgatgGTTGGGCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGGGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCt                                                                                      <  1:1025066/150‑1 (MQ=255)
                                                               tggttggtGGCTTTATCGGCTCAATTTTTGGGAAAAAAGCGGAAATTCGCGGCGGGCTGGTGCTGATCGGCAGCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGc                                                                                   <  1:671143/150‑1 (MQ=255)
                                                                 gttggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGGCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAg                                                                                 <  2:245694/150‑1 (MQ=255)
                                                                    ggtCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGAGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCg                                                                              >  2:1013855/1‑150 (MQ=255)
                                                                          tttAGCGGCTCAATTATGGGGAAAAACGCGGAAATTCTCGGCGGGCTGGTGCGGATCGGCAGCGGCGGCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACAct                                                                        <  1:1567528/150‑1 (MQ=255)
                                                                               cGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGAGCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCttt                                                                   <  1:110897/150‑1 (MQ=255)
                                                                                ggCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGCCGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCGGTCGACCACTTCtttt                                                                  >  2:1714234/1‑150 (MQ=255)
                                                                                ggCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCtttt                                                                  <  2:234485/150‑1 (MQ=255)
                                                                                  cTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTtgc                                                                <  1:217138/150‑1 (MQ=255)
                                                                                   tCAATTATTGGGAAAAAAGCGCAAATGCTCGGCGGGCGGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTtgct                                                               <  1:640187/150‑1 (MQ=255)
                                                                                    cAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGTCATCGGCGTCCAGATCCTCTTGACGCACTTCAACGGTTAATAAGAACCCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTttctg                                                              >  1:329611/1‑150 (MQ=255)
                                                                                      attattGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACATCTTTTGCTGcc                                                            >  1:1249062/1‑150 (MQ=255)
                                                                                       ttattGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCACATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCTGTCGAACACTTCTTTTGCTGCCa                                                           >  1:119258/1‑150 (MQ=255)
                                                                                        tattGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCGACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAg                                                          >  1:19676/1‑150 (MQ=255)
                                                                                           tGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAgtgt                                                       >  1:1003579/1‑150 (MQ=255)
                                                                                            gggAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGtt                                                      >  2:839478/1‑150 (MQ=255)
                                                                                            gggAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGATGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCATTTGCTGCAAGTGtt                                                      >  2:110230/1‑150 (MQ=255)
                                                                                              gAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTg                                                    >  2:1357833/1‑150 (MQ=255)
                                                                                               aaaaaaGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGGGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGc                                                   <  2:1134563/150‑1 (MQ=255)
                                                                                                aaaaaGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGcc                                                  <  2:1361827/150‑1 (MQ=255)
                                                                                                   aaGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGGCGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGCATATTAAAGTCCGTCTGGCAGTCGAGCACTTCTTTTGCTGCCAGTGTTTGGCAgc                                               >  1:1470474/1‑149 (MQ=255)
                                                                                                    aGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGAc                                              >  2:1241401/1‑150 (MQ=255)
                                                                                                    aGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGAc                                              <  2:34581/150‑1 (MQ=255)
                                                                                                             tctcGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCtt                                     <  1:615842/150‑1 (MQ=255)
                                                                                                              ctcGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTc                                    >  2:1527148/1‑150 (MQ=255)
                                                                                                                cggcggGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGc                                  >  1:654730/1‑150 (MQ=255)
                                                                                                                 ggcggGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCAATTCCACGGTTAATAAGAACGCTGCCAGAGGTTAATATTAAAGTCCGTCTGGAAGTCGAACACTTATTTAGCTTCCACTGTTTGACATACATCTGGCTTCGCa                                 >  1:1516133/1‑150 (MQ=255)
                                                                                                                     ggCTGGTGCTGATAGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGcc                             <  1:1594118/150‑1 (MQ=255)
                                                                                                                        tGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAgg                          >  1:1638479/1‑150 (MQ=255)
                                                                                                                         ggTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGAAAGGc                         >  2:988203/1‑150 (MQ=255)
                                                                                                                         ggTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAAAGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGc                         >  1:439435/1‑150 (MQ=255)
                                                                                                                            gCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCaaa                      >  1:191229/1‑150 (MQ=255)
                                                                                                                              tGATCGGCATCGGCGTCCAGATCCTCTGGACGCCCTTCCACGGATAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTCGCAGTCGAACACGGCGGTTGCTGCCAGTGTTTGCCCGACGTCGGGCTTCGCCCGCCAGGCACACg                    >  1:311582/1‑150 (MQ=255)
                                                                                                                               gctcggcatcggcGGCCAGATCCTCTGGCCGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACgg                   <  1:1515832/148‑1 (MQ=255)
                                                                                                                                atcggcagcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGc                  <  1:1625155/150‑1 (MQ=255)
                                                                                                                                 tcggcatcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCg                 <  1:1013855/150‑1 (MQ=255)
                                                                                                                                 tcggcatcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCg                 >  2:1247154/1‑150 (MQ=255)
                                                                                                                                    gcatcggcGTCCAGAACCTCTGGACGCACTTCCACGGTTAATAAGAAGGCTGCCAGAGGTGAATATTAACGTCCGTCTGGCAGTCGAACACTACTTTTGCTGCCAGTGATTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCa              <  1:1583254/150‑1 (MQ=255)
                                                                                                                                     catcggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCAt             >  2:494573/1‑150 (MQ=255)
                                                                                                                                      atcggcGTCCATATCCTTTGTAGGCACTTAGACGGTTAATAAGAACGCTGCCAGAGGTTAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTTCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATc            <  2:1516133/150‑1 (MQ=255)
                                                                                                                                         tgcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc         >  1:242320/2‑150 (MQ=255)
                                                                                                                                         tgcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc         >  1:190990/2‑150 (MQ=255)
                                                                                                                                         tgcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc         >  1:108872/2‑150 (MQ=255)
                                                                                                                                         ggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTTAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc         >  1:1557376/1‑150 (MQ=255)
                                                                                                                                         ggcGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGc         >  1:1369286/1‑150 (MQ=255)
                                                                                                                                          gcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCTAACGGCGTCATCTGCa        >  2:1050024/1‑150 (MQ=255)
                                                                                                                                          gcGTCCAGAGGATCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCa        >  2:1558805/1‑150 (MQ=255)
                                                                                                                                                aGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATg  <  2:819116/150‑1 (MQ=255)
                                                                                                                                                aGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATg  <  2:1657696/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACATTATCGCGACCGCATTGGCCATTGGTTGTGCAACCTTGATTATGTCAACATTAGGGATGATGGTTGGTCGCTTTATCGGCTCAATTATTGGGAAAAAAGCGGAAATTCTCGGCGGGCTGGTGCTGATCGGCATCGGCGTCCAGATCCTCTGGACGCACTTCCACGGTTAATAAGAACGCTGCCAGAGGTGAATATTAAAGTCCGTCTGGCAGTCGAACACTTCTTTTGCTGCCAGTGTTTGCCAGACTTCTGGCTTCGCACGCCAGGCAAACGGCGTCATCTGCAATAATG  >  NC_000913/1906082‑1906375

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