Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,443 T→A 5.3% intergenic (+64/‑102) trpS → / → zapE tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,504,4430TA 5.3% 437.6 / 10.6 148intergenic (+64/‑102)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (71/60);  new base A (3/5);  total (78/70)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 151 total aligned reads displayed.

CGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAC  >  NC_002947/1504295‑1504589
                                                                                                                                                      |                                                                                                                                                  
tgcttgatagcgttctaagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCttt                                                                                                                                                   >  1:1360790‑M1/98‑150 (MQ=255)
cGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCttt                                                                                                                                                   <  1:4030708/150‑1 (MQ=255)
  cttgatagcgttctaagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGc                                                                                                                                                 >  2:3757631‑M1/96‑150 (MQ=255)
   cGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGcc                                                                                                                                                <  1:5955307/150‑1 (MQ=255)
   cGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGcc                                                                                                                                                <  1:5955424/150‑1 (MQ=255)
     ggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCgg                                                                                                                                              <  1:1939246/150‑1 (MQ=255)
      gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGACCGAAGCCATGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt                                                                                                                                             <  2:3481194/150‑1 (MQ=255)
      gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCAAAGCCGgt                                                                                                                                             <  2:5199447/150‑1 (MQ=255)
      gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCGTTGCCGgt                                                                                                                                             >  2:7391745/1‑150 (MQ=255)
       ccGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCGTTGCCGgtg                                                                                                                                            >  2:513697/1‑150 (MQ=255)
         gCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGGGGCCCCGCAAAAGGCGCCGCCGGGGCATGGACCCCGCTGTGCCGgggt                                                                                                                                           >  1:3362480/1‑150 (MQ=255)
          cAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttc                                                                                                                                         >  1:467941‑M2/1‑147 (MQ=255)
          cAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCAAAGCCGgtgcca                                                                                                                                         <  2:6772262/150‑4 (MQ=255)
          cAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCAAAGCCGgtgcca                                                                                                                                         <  1:722654/150‑4 (MQ=255)
           gttctaagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGACCACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTc                                                                                                                                         >  1:1805647‑M1/87‑150 (MQ=255)
             ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcg                                                                                                                                      >  2:5124223‑M2/1‑144 (MQ=255)
             tctaaggtagtcatgattctttagcgaggcaacaaaaaacgggacgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCTGTGCATGGCACCGGCTTTGCCGGTTTtcgct                                                                                                                                      >  2:3079801‑M1/85‑149 (MQ=255)
              tcgagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                                     >  2:5402418‑M1/84‑150 (MQ=255)
              ctaagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                                     >  1:1503452‑M1/84‑150 (MQ=255)
              ctaagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                                     >  1:7448051‑M1/84‑150 (MQ=255)
               taagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCggg                                                                                                                                    >  2:6113627‑M1/83‑150 (MQ=255)
                gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggc                                                                                                                                   >  2:1148968‑M2/1‑141 (MQ=255)
                   ggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTcgc                                                                                                                                       <  1:625903‑M1/65‑1 (MQ=255)
                   ggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTcgc                                                                                                                                       >  2:625903‑M1/79‑143 (MQ=255)
                   aaCGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGCGGTGCATGGCACCCGCTTTGCCGGTGgtcgcgggcacg                                                                                                                                <  2:2743719‑M2/150‑13 (MQ=255)
                      gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCCCGGTGCATGGCACCCGCGTTGCCGGTGttcgcgggcacgccc                                                                                                                             <  2:5840797‑M2/150‑16 (MQ=255)
                       tcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTc                                                                                                                                         >  2:1762243‑M1/75‑137 (MQ=255)
                         gAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAAAAGGCGCCCCGGGGCCAGGCACCGGCGTTTCCGGTGttcgccggcccgcccccg                                                                                                                          <  2:62643‑M2/150‑19 (MQ=255)
                           aGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCAAGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgcgcc                                                                                                                        <  2:4573410‑M2/150‑21 (MQ=255)
                               cGTGTTGGGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACCGCAGGACTGAGGCTGCCAACGGGGGGGCGGGGGGGTGCCCGCTAAACCGGCGCCGGGGTGCATGGCACGGGCTTTGCGGGTGTTCGGGGGCACGCCCGCTCCcaca                                                                                                                    <  2:2216152/150‑1 (MQ=25)
                                   aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGCGGTGCAAGGCACCGGCTTTGCCGGTGttagcgggcacgcccgctcccacacaca                                                                                                                <  2:6989740‑M2/150‑29 (MQ=255)
                                     gCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGAGTGAGGCTGCAAAACGGTGGGAGGGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGa                                                                                                              <  1:4862231/150‑1 (MQ=255)
                                       ggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGGCCGCGGAAAAAGGGGCGCGGTGGATTTGACCCGGTTTTGCGGTGTTGGCGGGGACCCCCGCTCACAAAGACttttt                                                                                                            >  2:2511013/1‑146 (MQ=255)
                                       ggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg                                                                                                                                     <  1:4608904‑M2/125‑8 (MQ=255)
                                       ggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg                                                                                                                                     >  2:4608904‑M2/1‑118 (MQ=255)
                                       ggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGCGGTGCATGGCACCCGCGTTGGCGGTGTTAGCGGGCACCCCCGCTCCCCCAGACTGCtt                                                                                                            >  2:7090118/1‑150 (MQ=37)
                                          cATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGACACCGGCAAAGCCGgtgcc                                                                                                                                          >  1:1176066/1‑115 (MQ=37)
                                           aTTCCGCGCCTGAACGACGCCTTGGCACCACCGCTGGAGTGCGGGTGCCAATGGGTGGGCGCGGGGGTGCCCGCGAAACTGGCGCCGCGGTGCTTGGCACCGGCTTGGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCa                                                                                                        <  1:5494521/150‑1 (MQ=21)
                                                 agCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGGCCGCGAAAAAGGCGGCGCGGTGGGTGGCACCGGGGTTTGCGGTGTTGGCGGGCGCGCCCGCTCCCACAGAATGTTTTGCATTAAGc                                                                                                  >  2:7107064/1‑150 (MQ=25)
                                                     tGAACGAAGCCTTGGCACCACCGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCGGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                                              <  2:3858187/150‑1 (MQ=255)
                                                         cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGCAACAGGCGCCGCGGGGCCTGGCCCCCGCGGTGCCGGGGTTGGCGGGCCCGCCCCCGCCCACACCCTGCGTGCCCATAAGCGGGccgaa                                                                                          >  1:6765010/1‑147 (MQ=255)
                                                           acgCCTTGGCACCACCGCTGGACTGCGGCTGCAAAAGGGGGGGAGCGGGCGGGCCCGCGAAACCGGCGCCGCGGTGCTTGGCACGGGCTTTGCCGGTGTTGGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                        <  1:7258307/148‑1 (MQ=25)
                                                             gCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAAAAGGCGCCGCGGTGCATGGCACCCGCGTTTCCGGTGTTAGCGGGCACGCCCCCGCCCACAGACAGATTGCCCTTAAGCGGGCCCAACGCt                                                                                      >  1:6749171/1‑150 (MQ=255)
                                                                ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGCATGGCACCCGCGTTGCCGGTGTTCGCGGGGCCGCCCGCTCCCACAGCCGGATTGCCCATAAGCGGGCCGATAGCTACa                                                                                   >  2:7145940/1‑150 (MQ=255)
                                                                ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCGTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                   >  2:2409354/1‑150 (MQ=255)
                                                                        acaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGc                                                                           >  2:7390518/1‑150 (MQ=255)
                                                                          agccgGACTGTGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                         <  2:3709237/150‑1 (MQ=255)
                                                                          agcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTTAGCGGGCCTATCGCTACAGTGTCGGccc                                                                         >  2:5555275/1‑150 (MQ=255)
                                                                               gccgggggCTCCCAACGGGGGGGCGGGGGCGGGCCCGCGAACCCGGCGCCGGGGTGCTGGGCACCGGCTTTGCGGGGGTCCGGGGGCCCCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttt                                                                    <  1:3508689/144‑1 (MQ=255)
                                                                                 cTGCGGCTGCCAACGGGGGGGTGCGGGCGTGCCCGCGACACCGGCGCCGCGGTGCAGGGCACCGGCTTTGCCGGTGTCCGCGGCCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                  <  1:545341/150‑1 (MQ=37)
                                                                                 cTGAGGCTGCCAACGGGTGGGAGCGGGCGGGCCCGCGAAACTGGCGCCGCGGTGCTTGGCACCGGCTTTGCGGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                  <  2:1236821/150‑1 (MQ=37)
                                                                                 cTGAGGCTGCAAACGGGTGGGCGCGGGCGTGCCCGCTAAACAGGCGCGGCGGTGCATGGCACCGGCTTTGCGGGTGTGCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                  <  1:4198797/150‑1 (MQ=37)
                                                                                     ggCTGCCAAACGGTGGGAGCGGGGGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACGGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                                              <  2:211777/150‑1 (MQ=255)
                                                                                      gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCCGCGTTGCCGGTGTTCGCGGGCACGCCCGCGCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCGACAGTGACGGCCCGTTTTTTGTTGc                                                             >  2:3702020/1‑150 (MQ=255)
                                                                                                 GTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGa                                                  >  1:4957171‑M1/1‑150 (MQ=255)
                                                                                                    GGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                                               >  2:1498712‑M1/1‑150 (MQ=255)
                                                                                                     gcGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCTTGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                              >  2:1611837‑M1/3‑150 (MQ=255)
                                                                                                      cGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTGTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                             >  2:6659242‑M1/2‑150 (MQ=255)
                                                                                                           GCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                        <  2:3032961‑M1/150‑1 (MQ=255)
                                                                                                           GCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTGTGTGTTGCCTCGCTCACGAATCATGACTc                                        <  2:4385956‑M1/150‑1 (MQ=255)
                                                                                                               GCCCGCGAATCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCTCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCtt                                    >  2:6066631‑M1/1‑150 (MQ=255)
                                                                                                                cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcg                                   <  2:1138357‑M2/150‑106 (MQ=255)
                                                                                                                 ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcga                                  <  1:3313607‑M2/150‑107 (MQ=255)
                                                                                                                 ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcga                                  <  2:3774453‑M2/150‑107 (MQ=255)
                                                                                                                 CCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                                  <  1:4173650‑M1/150‑1 (MQ=255)
                                                                                                                 CCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                                  <  2:454507‑M1/150‑1 (MQ=255)
                                                                                                                  cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcgaa                                 <  2:2909715‑M2/150‑108 (MQ=255)
                                                                                                                  cgcgAAACAGGCGCCGCGGTGCATGGCACCGGC‑‑‑GCCTGTtcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcgaactgc                               <  2:3286645‑M2/150‑111 (MQ=255)
                                                                                                                  cgcgAAACAGGCGCCGCGGTGCATGGCACCGGC‑‑‑GCCTGTtcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcgaactgc                               <  2:1036845‑M2/150‑111 (MQ=255)
                                                                                                                  CGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                 <  2:842705‑M1/150‑1 (MQ=255)
                                                                                                                  CGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                 <  2:4158238‑M1/150‑1 (MQ=255)
                                                                                                                   gctaaaCAGGCGCCGCGGTGCATGGCACGGGCTTTGCCGGTTTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                                <  1:7320485/147‑1 (MQ=255)
                                                                                                                   GCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                                <  1:6356131‑M1/150‑1 (MQ=255)
                                                                                                                    CGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAAc                               <  2:6516476‑M1/150‑1 (MQ=255)
                                                                                                                     gAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcgaactg                              <  1:6075313‑M2/150‑111 (MQ=255)
                                                                                                                       AACCGGGGCCGGGTTGCTTGGCACCGCCTTGGCGGGTTTGCGCGGGCCCGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt                            >  1:1500328‑M1/1‑150 (MQ=255)
                                                                                                                         aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa                           >  2:4895011/1‑150 (MQ=255)
                                                                                                                         aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa                           >  2:1057233/1‑150 (MQ=255)
                                                                                                                         aCAGGCGCCGCGGTGCATGGCACCGGC‑‑‑GCCTGTtcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgt                                       <  2:534147‑M2/136‑103 (MQ=255)
                                                                                                                           aGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                         >  1:3786951/1‑150 (MQ=255)
                                                                                                                            ggCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcgtgcccgcga                                                                                                                         <  2:676193‑M2/53‑20 (MQ=255)
                                                                                                                             gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa                           <  2:4777109/146‑1 (MQ=255)
                                                                                                                             gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa                           >  1:4777109/1‑146 (MQ=255)
                                                                                                                             gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCaa                       <  2:6918906/150‑1 (MQ=255)
                                                                                                                             gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCaa                       >  1:6558485/1‑150 (MQ=255)
                                                                                                                             gcgcCGCGGTGCATGGCAACGGCTATGCCGGTGTACGCGCGCACGCCCGATCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGACTCGCTACCGAATCATGCCTCCCTACGAACGCTATCaa                       >  2:5933578/1‑150 (MQ=255)
                                                                                                                               gccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc                     >  1:2387640/1‑150 (MQ=255)
                                                                                                                                  gcgGTGCATGGCATCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGa                  >  1:5193707/1‑150 (MQ=255)
                                                                                                                                      gTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCt               >  1:7309344/1‑150 (MQ=255)
                                                                                                                                             gCACCGGCTTTGCCGGTGttcgcgggcgagcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcgaactgcggcccgtcggcaacctgctcg        <  2:2673891‑M2/150‑133 (MQ=255)
                                                                                                                                             gCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGt        >  1:6143896/1‑150 (MQ=255)
                                                                                                                                              cACGGGCTTTGCCGGTGTTGGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTc       <  1:782589/150‑1 (MQ=255)
                                                                                                                                               aCCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGAt                                                                                                             >  2:7138659/1‑47 (MQ=255)
                                                                                                                                               aCCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGAt                                                                                                             <  1:7138659/47‑1 (MQ=255)
                                                                                                                                               aCCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc      >  1:1099764/1‑150 (MQ=255)
                                                                                                                                               aCCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc      >  2:3468769/1‑150 (MQ=255)
                                                                                                                                               aCCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc      >  1:1099758/1‑150 (MQ=255)
                                                                                                                                               aCCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCGCGCAAACGAAACCAGACTCCCGTAGCACCCGATCAAGCAGAACAGAAACCTcc      >  2:1876669/1‑150 (MQ=37)
                                                                                                                                                   gCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc  >  1:4015601/1‑150 (MQ=255)
                                                                                                                                                   gCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc  >  1:2123754/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAC  >  NC_002947/1504295‑1504589

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