Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,695 C→G 9.0% 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,6950CG 9.0% 873.7 / 18.3 268intergenic (+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 (151/88);  new base G (2/27);  total (153/115)
Fisher's exact test for biased strand distribution p-value = 2.87e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

CACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATA  >  NC_000913/3608547‑3608844
                                                                                                                                                    |                                                                                                                                                     
  cGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATAGGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                    <  1:463156/150‑1 (MQ=255)
   gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   >  2:3450251/1‑150 (MQ=255)
   gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   <  1:1279832/150‑1 (MQ=255)
   gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   <  2:1505415/150‑1 (MQ=255)
     tgtgAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                 <  2:444640/150‑1 (MQ=255)
      gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                <  1:3857908/150‑1 (MQ=255)
      gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                <  1:3794739/150‑1 (MQ=255)
       ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               <  1:4545960/150‑1 (MQ=255)
       ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               <  1:3621515/150‑1 (MQ=255)
       ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               <  1:2828542/150‑1 (MQ=255)
         aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                             <  1:4388054/150‑1 (MQ=255)
          gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGTCGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                            >  2:1915114/1‑150 (MQ=255)
            gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAgggga                                                                                                                                          <  1:4456920/150‑3 (MQ=255)
              tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                        <  2:1341320/150‑1 (MQ=255)
               gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                       >  2:3532605/1‑150 (MQ=255)
               gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                       <  1:3393791/150‑1 (MQ=255)
               gACCGGGTGGGGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                       <  1:3287901/150‑1 (MQ=255)
                aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                      <  1:1195611/150‑1 (MQ=255)
                  cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                    >  2:4320677/1‑150 (MQ=255)
                  cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                    >  1:2921538/1‑150 (MQ=255)
                   gggTTGTGGCTGGCAGTGCTGGCAGTTACGGGGGCGACGGTGCTGGTGACAGCGATTGCGTTAAGATTGTGGCGCAGGAGATAAGGCAAACCGATCGCAACATTGTGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                   <  1:2107097/150‑1 (MQ=255)
                   gggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCTAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                   <  2:2416005/150‑1 (MQ=255)
                    ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                  <  1:2015811/150‑1 (MQ=255)
                         tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTTGTGACAGCGAATGCGTTAAGATTGTTGAGCAGGAGATAATGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCATCTGGGGAGAGGGTTGGGGTGAGGGGaa                                                                                                                             >  2:1741844/1‑150 (MQ=255)
                         tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                             >  1:3548353/1‑150 (MQ=255)
                           gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                           >  1:3145560/1‑150 (MQ=255)
                           gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                           >  1:511984/1‑150 (MQ=255)
                           gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                           >  1:594945/1‑150 (MQ=255)
                           gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                           >  2:3381315/1‑150 (MQ=255)
                            ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                                          <  2:3971272/150‑1 (MQ=255)
                            ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                                          >  1:1728945/1‑150 (MQ=255)
                             tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                         <  2:1527343/150‑1 (MQ=255)
                               gcATTGCTGGCAGATACGGGGGCGACGGTGCTGGTTAAAGCGAATGCTTTAAGATTTTTGCGAAGGATATAAGGAAAACCGATCGCAACATTGAGCGCTATCGGTCCCCTCGCCCCTCTGGGTATAGGGTTAGGGTGAGGGGAAAAggcg                                                                                                                       <  2:1041623/150‑1 (MQ=37)
                                     cTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCg                                                                                                                 >  2:3126519/1‑150 (MQ=255)
                                      tGGCTGATACGGGGGCGACGGTGCGGGTGACAGCGATTGGGTTATTATTGTTGCTCAGGAGATAAGGCAAACCGCTCGCAACATTGTGCGCGTTCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGAGGGTAAAAGGCGGCATCGa                                                                                                                <  2:3908963/150‑1 (MQ=37)
                                      tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGGGATAAGGCAAACCGATCGCAACAGTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                >  2:422552/1‑150 (MQ=255)
                                      tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                >  1:1407653/1‑150 (MQ=255)
                                               cGGGGGCGCGGGTGCTGGTTACTGCTAATGGTTTAAGTTTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATGGGTCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCa                                                                                                       <  2:552757/150‑1 (MQ=37)
                                                         ggTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCa                                                                                             >  1:706801/1‑150 (MQ=255)
                                                           tgctgcTGACAGCGAATGCGTTAAGGTTGTTGTGCAGGAGATAAGGCAAACCGATCGCAACAGTGAGCGCGATCGGTCCCCTCGCCCCACTGGGGAGAGGGATAGGGTGAGGGGACCAGGCGGCATCGAAGCCAATCCCCCCCGATCAAc                                                                                           >  2:3655365/1‑150 (MQ=255)
                                                             ggggTGACTGCGATTGGTTTATGTTTGTGGCGGTGGTTATCTGGCAACCCGATCGCAACTTGGCGCGCGTGGGGGCCCCTCGCCCGTGTGGGGTGTGGGTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCg                                                                                         <  1:3573698/148‑1 (MQ=255)
                                                               ggTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGcc                                                                                       >  2:2252488/1‑150 (MQ=255)
                                                                 tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                     >  2:4614718/1‑150 (MQ=255)
                                                                  gACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                    >  1:2379137/1‑150 (MQ=255)
                                                                  gACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                    >  2:3687910/1‑150 (MQ=255)
                                                                  gACAGCGAATGCGTTAAGATTGTTGCACAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                    >  1:3108576/1‑150 (MQ=255)
                                                                    cAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTGAGGGTGAGGGGACGAGGCGGCCTCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                  >  2:3740269/1‑150 (MQ=255)
                                                                      gCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                                >  1:149948/1‑150 (MQ=255)
                                                                         tttGCTTTATGTTTGTTGCGCTGTATATAAGGCAAACCGATCGAAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGTTGTGGGTTCGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATc                                                                             <  2:1065097/148‑1 (MQ=255)
                                                                              gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  2:2321043/1‑150 (MQ=255)
                                                                              gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  1:696990/1‑150 (MQ=255)
                                                                               ttatttttgttgCGCTGTATATCCGGAAACCCGATCGCAACTTTGCGCGGTTTGGGTCCCCGGGCCCCTGGGGGGTGTGTTTTGGGGGGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                       <  2:2844696/144‑1 (MQ=255)
                                                                                  tgATTTTTGCTCTGGAGATAAGGCAACCCGATCGCAACTTTGAGGGCGATGGGTCCCCTCGCCCCTCTGGGGTGTGGTTAAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                    <  2:1111721/149‑1 (MQ=255)
                                                                                    tttgttgCGCTGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATGGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                                  <  2:592375/149‑1 (MQ=255)
                                                                                     ttgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                 >  2:3331050/1‑150 (MQ=255)
                                                                                         tgCGCTGGAGATAAGGAAAACCGATCGAAACATTGAGCGCTATCGGTCCCCTCGCCCCTCTGGGGAGTGGTTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                             <  2:289944/150‑1 (MQ=39)
                                                                                          gcgcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAg                                                            >  2:839700/1‑150 (MQ=255)
                                                                                             cAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTcc                                                         >  2:3206584/1‑150 (MQ=255)
                                                                                                gagaTAAGGCAACCCGATCGCAACATTGCGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGTGGGTTAGGTTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                      <  1:3534988/150‑1 (MQ=255)
                                                                                                 agaTAAGGCAACCCGATCGCAACTTTGAGCGCGATGGGCCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                     <  2:2522210/150‑1 (MQ=255)
                                                                                                   aTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                   >  2:4696974/1‑150 (MQ=255)
                                                                                                   aTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                   >  2:4780723/1‑150 (MQ=255)
                                                                                                    tAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtc                                                  >  1:1908862/1‑150 (MQ=255)
                                                                                                       ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                               >  2:1421027/1‑150 (MQ=255)
                                                                                                        gCAAACCTATGGCAACTTTGCGCGCGATCGGTCCCCGCGCCCCTCTGGGGAGTGGGTTCGGGTGAGGCGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              <  2:4188782/150‑1 (MQ=255)
                                                                                                        gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              >  2:4609648/1‑150 (MQ=255)
                                                                                                        gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              >  2:369675/1‑150 (MQ=255)
                                                                                                        ccAACCCGACCGCAACATGGCGCGCGAGCGGCCCCCTCGCCCCCGGGGGTAGCGGGTTCGGGGGAGGGCACACGGCGGCATCGCCGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              <  1:228343/149‑1 (MQ=255)
                                                                                                         cAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                             <  2:2422504/150‑1 (MQ=255)
                                                                                                          aacccGAGCGCAAATTTGCGCGCGATCGGTCCCCCCGCCCCTGTGGGGTGGGGGTTTGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            <  2:1891687/147‑1 (MQ=255)
                                                                                                          aaaCCGATCTCACCATTGAGCGCGAGCGGTCCCCTCGCCCCGGTGGGGAGTGGGTTGGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            <  2:2110591/150‑1 (MQ=255)
                                                                                                          aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            >  1:4018811/1‑150 (MQ=255)
                                                                                                          aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            >  2:2946514/1‑150 (MQ=255)
                                                                                                           aaCCTATCGCAACATTGAGCGCGATCGGTCCCCGCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAg                                           <  2:3533200/150‑1 (MQ=255)
                                                                                                           aaCCGTGCGCAACTTTGCGCGCGTTGGGGCCCCTCGCCCCTCTGGGTTGTGGGTTTGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAg                                           <  1:2978270/150‑1 (MQ=255)
                                                                                                                aTCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcc                                      >  1:2661771/1‑150 (MQ=255)
                                                                                                                aTCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcc                                      >  2:1645300/1‑150 (MQ=255)
                                                                                                                   gCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                                   >  2:1087957/1‑150 (MQ=255)
                                                                                                                    cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                  >  1:2580807/1‑150 (MQ=255)
                                                                                                                    cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                  >  1:4124960/1‑150 (MQ=255)
                                                                                                                    cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                  >  1:863530/1‑150 (MQ=255)
                                                                                                                     aataTTGCGCGCGAGGGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGTGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                 <  2:4532517/147‑1 (MQ=255)
                                                                                                                     aaCATGGCGCGCGATGGGTCCCCCCGCCCCTCTGGGGTGTGGGTTCGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                 <  2:3472669/150‑1 (MQ=255)
                                                                                                                        aTTGCGCGCGATCGGTCCCCTCGCCCCTGTGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                              <  1:726238/150‑1 (MQ=255)
                                                                                                                            aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                          <  1:1971787/150‑1 (MQ=255)
                                                                                                                               gcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATa                       >  1:2517165/1‑150 (MQ=255)
                                                                                                                                 gATGGGTCCCGCCGCCCCTCGGGGTAGAGGGTTCGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                     <  1:1520450/150‑1 (MQ=255)
                                                                                                                                 gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                     >  2:1999044/1‑150 (MQ=255)
                                                                                                                                  aTCGGTCCCCCCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                    <  1:2634462/150‑1 (MQ=255)
                                                                                                                                     ggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa                 >  1:4229049/1‑150 (MQ=255)
                                                                                                                                      taacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                >  2:1145252/4‑150 (MQ=255)
                                                                                                                                         cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa             >  1:2404329/1‑150 (MQ=255)
                                                                                                                                           cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc           >  1:1035912/1‑150 (MQ=255)
                                                                                                                                            tCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc          <  2:2993793/150‑1 (MQ=255)
                                                                                                                                            tCGCCCCGCTGGGTTGTGTGTTTGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGAAACCAGTTCGTGTTCCATAAAGGTACAAAAcc          <  2:2695545/150‑1 (MQ=255)
                                                                                                                                              gCCCCTCAGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACTATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTACGTGTTCCATAAAGGTACAAAACCCa        >  1:2837575/1‑150 (MQ=255)
                                                                                                                                               ccccTCTGGGGAGAGGGTTAGGGTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg       <  1:289418/150‑1 (MQ=255)
                                                                                                                                                   tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat   >  2:923070/1‑150 (MQ=255)
                                                                                                                                                    ctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAata  >  1:2548547/1‑150 (MQ=255)
                                                                                                                                                    ctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAata  >  1:2775591/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: