Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,695 C→G 13.3% intergenic (+46/+56) zntA → / ← tusA zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,6950CG13.3% 1517.5 / 74.3 498intergenic (+46/+56)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (257/167);  new base G (6/64);  total (263/235)
Fisher's exact test for biased strand distribution p-value = 2.59e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 520 total aligned reads displayed.

CACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATA  >  NC_000913/3608547‑3608844
                                                                                                                                                    |                                                                                                                                                     
cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                      >  1:4170014/1‑150 (MQ=255)
 acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTg                                                                                                                                                     >  1:4931483/1‑150 (MQ=255)
  cGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                    <  2:4569980/150‑1 (MQ=255)
   gCTGTGAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCACACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   <  2:781417/150‑1 (MQ=255)
       ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               <  2:4821265/150‑1 (MQ=255)
        tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATCCGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              <  1:310177/150‑1 (MQ=255)
        tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              <  1:1339473/150‑1 (MQ=255)
        tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              >  1:3961105/1‑150 (MQ=255)
         aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                             >  1:3836990/1‑150 (MQ=255)
         aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                             >  1:4656039/1‑150 (MQ=255)
         aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                             <  1:5055145/150‑1 (MQ=255)
          gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAggg                                                                                                                                            <  2:3799276/150‑1 (MQ=255)
          gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAggg                                                                                                                                            <  2:5423851/150‑1 (MQ=255)
            gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                          >  1:6399313/1‑150 (MQ=255)
              tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                        >  1:5672505/1‑150 (MQ=255)
               gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAgg                                                                                                                                       >  1:8396804/1‑150 (MQ=255)
                 ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAgagggtgagggt                                                                                                                                     >  2:5566987/1‑150 (MQ=255)
                  cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                    >  2:6953494/1‑150 (MQ=255)
                    ggTTGTGGCTGGCAGTGCTGGCAGATCCGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                  <  2:5302991/150‑1 (MQ=255)
                    ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                  <  1:6290507/150‑1 (MQ=255)
                    ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                  <  2:5971276/150‑1 (MQ=255)
                     gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                 <  1:680315/150‑1 (MQ=255)
                     gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTACTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                 <  2:5849549/150‑1 (MQ=255)
                      ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                <  1:5537054/150‑1 (MQ=255)
                           gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                           >  2:3596696/1‑150 (MQ=255)
                             tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAAgg                                                                                                                         <  2:2426848/150‑1 (MQ=255)
                              ggcTGTGCTGGCTTTTACGGGGGCGAGGGTGGTGGTGACAGCTATTGCTTTATTTTTGTTGCGAAGTATTTATGTAAAACCGATCGCAACTTTGTGCGCGATGGGTCCCCTCGCCCCTCTGGGGTGAGGTTTAGGGTGTGGGGAAAAggc                                                                                                                        <  2:3388520/150‑1 (MQ=16)
                              ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCTGGAGATAAGGCAAACCGATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGTTGAGGGGAAAAggc                                                                                                                        <  2:3012103/150‑1 (MQ=39)
                              ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggc                                                                                                                        >  2:1774964/1‑150 (MQ=255)
                              ggcAGTGCTGGCAGATACGGGGGCGACGGGGCTGGTGACCGCGAATGCGTTAAGATGGTGGCGCAGGAGCTAAGGCAACCCGATCGCACCATTGAGGGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGTGAAAAggc                                                                                                                        <  2:1256036/150‑1 (MQ=38)
                              ggcAGGGCGGGCAGATACGGGGGCGCCGGTGATGGTGACAGCGAATGCTTTAAGATTGGTGCGCAGGAGATAAGGCAAACCGATCCCAACTTTGAGCGCGATCGGTCCCCTCGCCCCGCTGGGGAGAGGGTTCGGGTGAGGGGAAAAggc                                                                                                                        <  2:6631519/150‑1 (MQ=255)
                                  gTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCa                                                                                                                    >  2:2183360/1‑150 (MQ=255)
                                     cTGGCAGATACGGGGGCGACGGTGCTGGTTACAGCGATTGCTTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACTTTGTGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGTTTAGGGTGAGGGGAAAAGGCGGCATCg                                                                                                                 <  2:8547294/150‑1 (MQ=38)
                                        gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                              >  1:2900290/1‑150 (MQ=255)
                                         catatcggggggCGCCGGTGGTGGGGCCCGCGATTGCGTTAAGTTTGTTGCGCAGGAGATATGGCAAACCGTTCGCAACATGTTGCGCGATGGGTCCCCGCGCCCCTGTGGGTAGTGGGTTCGGGTGTGGGGAAAAGGCGGCATCGAAGc                                                                                                             <  2:9080923/143‑1 (MQ=255)
                                               cGGGGGCGAGGGTGGTGGTGACAGCTAATGCGTTAATATTGTTGCGCAGGAGATACGGCAAACCGACCGCAACTTTGAGGGCTTTCGGTCCCCTCGCCCCTGTGGGTTTTGGGTTCGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCa                                                                                                       <  2:3880131/150‑1 (MQ=255)
                                               cGGGGGCGACGGTGCTGGTGACCGCGAATGCGTTAAGTTTGTTGCGCTGGAGATAAGGCAAACCGATCGCAACATTGAGCGCTATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCa                                                                                                       <  2:4632427/150‑1 (MQ=255)
                                                 ggggggAGGGTGCTGGTGACTGCGAATGCGTTAAGATTGTTGCGCTGTAGATAAGGCAAACCGATGGAAACATTGAGGGGGTTGGGTCCCCTCGCCCCTGTGGGGAGTGGGTTCGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                     <  2:138986/150‑1 (MQ=255)
                                                 ggggCGCGGGTGGTGGTGACCGCTAATGGTTTATGTTTGTGGCGCCGGAGATCAGGCAAACCGATCGAAAATTGGCGCGCGCTGGGGCCCCTCGCCCCTGGGGGGTGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                     <  2:7363042/150‑1 (MQ=255)
                                                 ggggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCACCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                     >  2:7577562/1‑150 (MQ=255)
                                                  gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                    >  1:4475819/1‑150 (MQ=255)
                                                    gCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcccc                                                                                                  >  2:1503358/1‑150 (MQ=255)
                                                      gACGGTGCTGGGGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATCAGGCAAACCGATCGCAACAGTGAGCGCGTTCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAGGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                                >  1:2672943/1‑150 (MQ=38)
                                                       aCGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTAt                                                                                               >  2:5638688/1‑150 (MQ=255)
                                                         ggTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCa                                                                                             >  2:4986545/1‑150 (MQ=255)
                                                         ggTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCa                                                                                             >  2:6069888/1‑150 (MQ=255)
                                                           tGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                           >  1:8588569/1‑150 (MQ=255)
                                                              tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                        >  1:3719864/1‑150 (MQ=255)
                                                               ggTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGcc                                                                                       >  1:1956729/1‑150 (MQ=255)
                                                                    cAGCGAATGCTTTAAGATTGTTGCGCTGGAGATAAGGCACACCGATCGCAACATGGAGCGCGATCGGTCCCCTCGCCCCTGTGGGGTGTGGGTTAGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                  <  2:2904057/150‑1 (MQ=255)
                                                                    cAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                  >  2:8888341/1‑150 (MQ=255)
                                                                     aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                 >  2:5137796/1‑150 (MQ=255)
                                                                     aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                 >  1:5625871/1‑150 (MQ=255)
                                                                       cGAATGGTTTAAGTTTGTGGCGCAGGATATACGCCAACCCGATCGCAACTTTGAGCGCGTTCGGTCCCCCCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                               <  2:170520/150‑1 (MQ=255)
                                                                        gaaggtgTTAATTTTGTTGCGCAGTCGATAAGCCAAACCTCTCGCAAATTTGAGCGCGAACGGTCCCCTCGCCCCTCTGGGGATTGGGTTCGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGTCCCTATCAACCGCCTTTACGaaa                                                                              <  1:8901368/144‑2 (MQ=25)
                                                                          aTGCGTTAAGTTTGTTGCGCAGGAGATACGGAAAACCGAGCGCAACTTTGAGCCCCATGGGTCCCCCCGCCCCTGTGGGTAGAGGGTTAGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCa                                                                            <  2:8179476/150‑1 (MQ=255)
                                                                             cGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGACCCTATCAACCGCCTTTACGAATCAAAt                                                                         >  1:8000726/1‑150 (MQ=255)
                                                                              gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  1:1072043/1‑150 (MQ=255)
                                                                              gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  2:1909726/1‑150 (MQ=255)
                                                                                       gttgCGCAGGAGATAAGGAAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                               >  2:2289687/1‑150 (MQ=255)
                                                                                           cgcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGt                                                           <  2:9272217/150‑1 (MQ=255)
                                                                                               ggAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCAt                                                       >  1:8085324/1‑150 (MQ=255)
                                                                                                gagaTATGGCAAACCGATCGCAACTTTGAGCGCGATCGGTCCCGTCGCCCCTGTGGGTTGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                      <  2:8745716/150‑1 (MQ=255)
                                                                                                gagaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                      >  1:2606754/1‑150 (MQ=255)
                                                                                                gagaAACGGCAAACCGATCGAAAAATGGCGCGCGATCGGCCCCCTCGCCCCTCTGGGGTGAGGGTTCGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                      <  1:6089894/150‑1 (MQ=255)
                                                                                                 agaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                     >  1:3549968/1‑150 (MQ=255)
                                                                                                  gaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg                                                    >  1:2250567/1‑150 (MQ=255)
                                                                                                     acggaaaaCGGAGCGAAACATTGCGCGCGATCGGTCCCCTCGCCCCTGGGGGGAGTGGGTTAGGTTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                 <  2:8386116/145‑1 (MQ=255)
                                                                                                      aGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtctc                                                >  1:1499146/1‑150 (MQ=255)
                                                                                                      aGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTGTACGAATCAAATAACGATCAGGCAGACCATCCGtctc                                                >  1:6473195/1‑150 (MQ=255)
                                                                                                       ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                               >  2:7201174/1‑150 (MQ=255)
                                                                                                        gaaaaCCGTTCGCAAATTTGTGCGCGATCGGTCCCCCCGCCCCTCGGGGGTGTGGTTTCGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              <  1:7746147/148‑1 (MQ=255)
                                                                                                        gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              >  1:2495881/1‑150 (MQ=255)
                                                                                                        gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              >  1:760451/1‑150 (MQ=255)
                                                                                                        gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              >  1:6630881/1‑150 (MQ=255)
                                                                                                        gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              >  1:5401949/1‑150 (MQ=255)
                                                                                                          aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            >  2:2680797/1‑150 (MQ=255)
                                                                                                          aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            >  1:7506232/1‑150 (MQ=255)
                                                                                                                 tCGCAACTTTGTGCGCGTTGGGTCCCGGCGCCCCTGTGGTGAGCGGTTGTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCATCCCCTATCAACCGCCTTTACGAATCAAAGAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                     <  2:4317679/150‑1 (MQ=255)
                                                                                                                 tCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                     <  2:6263220/150‑1 (MQ=255)
                                                                                                                   gaacctttGCGCGCGATCGGCCCCCTCGCCCCTGTGGGTAGTGGTTTTGGGTGTGGGCACAAGGCGGTATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                                   <  1:2338186/144‑1 (MQ=255)
                                                                                                                        aTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                              >  2:7081072/1‑150 (MQ=255)
                                                                                                                         ttGCGCGCGCTCGGGCCCCTCGCCCCTCGGGGGTGCGGGTTCGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                             <  2:8299941/150‑1 (MQ=255)
                                                                                                                             gcgcGATCGGACCCCTCGCCCCTCTGGGGAGTGGTTTCGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCa                         <  1:5273165/150‑1 (MQ=255)
                                                                                                                              cgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                        >  1:2749313/1‑150 (MQ=255)
                                                                                                                                 gctCGCTCCCCCCGCCCCTGTGGGTTGTGGGTTTGGGTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                     <  1:3260053/148‑1 (MQ=255)
                                                                                                                                  aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                    >  1:7096303/1‑150 (MQ=255)
                                                                                                                                      gcccccTCGCCCCTGTGGTTTGTGGTTTTGGGTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                <  2:4188254/148‑1 (MQ=255)
                                                                                                                                      gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                >  2:1025276/1‑150 (MQ=255)
                                                                                                                                       aacccTCTCCCCAGAGGGGCGCGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCCATCAGCCCCTATCAACCGCCTGTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCACCCAGTTCGTGTTCCATAAAGGTACa               >  1:1720558/3‑150 (MQ=255)
                                                                                                                                          ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa            >  2:3728930/1‑150 (MQ=255)
                                                                                                                                           cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc           >  2:6345924/1‑150 (MQ=255)
                                                                                                                                            tCGCCCCTCGGGGGAGAGGGGGAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATACAGGTACAAAAcc          >  2:318665/1‑150 (MQ=255)
                                                                                                                                              gCCCCTGTGGGGAGTGGTTTTGGTTGCGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        <  1:7331298/150‑1 (MQ=255)
                                                                                                                                              gCCCCTCTGGGGAGAGGGTTTGGGGGTGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        <  1:3400452/150‑1 (MQ=255)
                                                                                                                                              gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        >  1:8720320/1‑150 (MQ=255)
                                                                                                                                              gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        >  1:6013502/1‑150 (MQ=255)
                                                                                                                                               ccccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg       >  1:9084998/1‑150 (MQ=255)
                                                                                                                                                cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAgg      >  2:8917671/1‑150 (MQ=255)
                                                                                                                                                  ctctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa    >  1:7948096/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATA  >  NC_000913/3608547‑3608844

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: