Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,447 C→G 6.3% intergenic (+68/‑98) 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,4470CG 6.3% 617.0 / 3.7 196intergenic (+68/‑98)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (107/74);  new base G (4/11);  total (111/85)
Fisher's exact test for biased strand distribution p-value = 2.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

CGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCC  >  NC_002947/1504298‑1504596
                                                                                                                                                        |                                                                                                                                                     
cGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGcc                                                                                                                                                       <  1:4545824/150‑1 (MQ=255)
 gggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                                                      <  2:850702/150‑1 (MQ=255)
 gggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                                                      <  1:4497701/150‑1 (MQ=255)
  tgatagcgttctaagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtcgGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                                                      >  1:2690181‑M1/94‑150 (MQ=255)
  tgatagcgttctaagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtcgGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                                                      >  1:2802501‑M1/94‑150 (MQ=255)
   gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt                                                                                                                                                    <  1:6401768/150‑1 (MQ=255)
   gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt                                                                                                                                                    <  1:6050375/150‑1 (MQ=255)
    ccGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCAAAGCCGgtg                                                                                                                                                   <  1:3921628/150‑1 (MQ=255)
      gCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGtt                                                                                                                                                 <  2:4082021‑M2/150‑3 (MQ=255)
       gcgttctaagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGtt                                                                                                                                                 >  2:5453388‑M1/89‑150 (MQ=255)
       cAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttc                                                                                                                                                >  2:3384470‑M2/1‑147 (MQ=255)
        aGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCAAAGCCGgtgccat                                                                                                                                               <  1:1411942/150‑5 (MQ=255)
         gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgc                                                                                                                                              >  2:1014262‑M2/1‑145 (MQ=255)
         gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgc                                                                                                                                              >  1:5646017‑M2/1‑145 (MQ=255)
          ttctaagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTcgc                                                                                                                                              >  2:4760940‑M1/86‑150 (MQ=255)
          ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGGGGCCCCCGAAAAAGCCCCCCCGGGGCTTGCCACCGCGTTTGCggggtttgcg                                                                                                                                              >  1:5985224/1‑143 (MQ=255)
          ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcg                                                                                                                                             >  1:643921‑M2/1‑144 (MQ=255)
          ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcg                                                                                                                                             >  2:2820942‑M2/1‑144 (MQ=255)
          ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcg                                                                                                                                             >  2:3781505‑M2/1‑144 (MQ=255)
          ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcg                                                                                                                                             >  1:6371058‑M2/1‑144 (MQ=255)
          ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcg                                                                                                                                             >  1:5550884‑M2/1‑144 (MQ=255)
          ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcg                                                                                                                                             >  2:5666602‑M2/1‑144 (MQ=255)
           tctacgggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgatcggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAACACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTcgc                                                                                                                                              >  1:3129231‑M1/85‑150 (MQ=255)
           tctaagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTcgcg                                                                                                                                             >  1:6472789‑M1/85‑150 (MQ=255)
           tctaagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTcgcg                                                                                                                                             >  2:5288078‑M1/85‑150 (MQ=255)
           tctaagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTcgcg                                                                                                                                             >  2:535950‑M1/85‑150 (MQ=255)
           tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg                                                                                                                                            >  1:4085441‑M2/1‑143 (MQ=255)
           tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg                                                                                                                                            >  2:3504112‑M2/1‑143 (MQ=255)
            ctaagagagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                                            >  1:3212154‑M1/84‑150 (MQ=255)
            cGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcggg                                                                                                                                           >  2:3417824‑M2/1‑142 (MQ=255)
             taagggagtcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCggg                                                                                                                                           >  2:4837790‑M1/83‑150 (MQ=255)
             gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGGAACAGGCGCCGCGGGGCATGGCACCGGCTTTGGCGGTGttagcgggc                                                                                                                                          <  2:1266300‑M2/150‑10 (MQ=255)
             gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACACGCGCCGCCGGGCATGGCACCCGCGTTTCCGGTGttggccggc                                                                                                                                          <  2:1821478‑M2/150‑10 (MQ=255)
              aCAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggca                                                                                                                                         >  1:5076965‑M2/1‑140 (MQ=255)
               cAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTGTGCCGGTGttcgcgggcac                                                                                                                                        >  2:6378130‑M2/1‑139 (MQ=255)
                     tcatgattcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                                            >  1:1776093‑M1/75‑141 (MQ=255)
                       aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcgtgcccgcga                                                                                                                                >  2:2540630‑M2/1‑131 (MQ=255)
                          attcgttagcgaggcaacaaaaaacgggccgtcactgtagcgataggcccgcttaatggcaatcagtctGTGGGAGCGGGCGTGCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCca                                                                                                                             >  2:561076‑M1/70‑150 (MQ=255)
                             gttttggggCCTTTCTGAGCCTTAACGATGCGTGGGCACCACCGCTGGATTGTGGGGCCAAACGGGGGGGCGGGGGGGGGCCCGCTAAACTGGCGCGGGGGTGCTTGGCACGGGCTTTGCGGGTTTTGGGGGGCACGCCCGCTCCCACAg                                                                                                                          <  2:1755360/143‑1 (MQ=255)
                                       cATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCCCGGCTTTGGCGGGGTTCGCGGGCCCGCCCGCACCCACAGACTGAtttgc                                                                                                                >  2:1544815/1‑147 (MQ=37)
                                        aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCa                                                                                                               <  2:430145/150‑1 (MQ=255)
                                          tactagCCTGACCGAAGCCTTGGCACCACAGCAGGACTGTGGCTGCCAAAGGGTGGGTGCGGGCGGGCCCGCGAAACTGGCGCCGCGGTGCTTGGCACCGGCTTTGCCGTTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCAtt                                                                                                             <  2:1224445/146‑1 (MQ=37)
                                                  tGAACGACGCCTTGGCACCACAGCAGGACTGAGGCTGCAAAACGGTGGGAGCGGGCGGGCCCGCGACCCCGGCGCCGCGGGGCTTGGCACGGGCTTTGCGGGTGTGCGGGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                                                     <  2:2438841/150‑1 (MQ=255)
                                                  tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTTTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                                                     <  1:3363073/150‑1 (MQ=255)
                                                    aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCt                                                                                                   <  2:4419571/150‑1 (MQ=255)
                                                       gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGCAACAGGCGCCGCGGGGCATGGCCCCGGCGTTGGCGGGGTTCGCGGGCCCGCCCCCGCCCACCGACTGGTTGGCCTTAAGCggggcgatg                                                                                                >  1:6837017/1‑144 (MQ=255)
                                                             ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg                                                                                                                                            <  1:1467321‑M2/100‑8 (MQ=255)
                                                             ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                          <  2:4474718/150‑1 (MQ=255)
                                                             ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCGTTGCCGGTGTTCGCGGGCCCGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                          >  2:6818504/1‑150 (MQ=255)
                                                                gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCCAAAAAGGCGCCGCGGTGCATGGCACCCGCATTGCCGGTGTTGGCGGGCCCGCCCCCGCCCCCACACAGATTGCCAATAAGCGGGCCGATCGCAACCGTg                                                                                       >  2:6057917/1‑150 (MQ=255)
                                                                              cggcggCTGCCACACGGTGGGCGCGGGCGGGCCCGCGAAACCGGCGCCGCGGGGCAGGGCACCGGCTTGGCCGGGGTGCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                         <  2:379105/146‑1 (MQ=37)
                                                                               tGAGACTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCCCGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttg                                                                        >  2:4508909/1‑150 (MQ=255)
                                                                                aatggcaatcagtctGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTcgc                                                                                                                                              >  1:4836440‑M1/16‑80 (MQ=255)
                                                                                    ggtgcAAACGGGTGGGCGGGGGGGGGCCCGCCAAACGGGCGCCGCGGGGCTGGGCCCGGGCTTTGGGGGTTTGCGGGGCCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGc                                                                    <  1:1938934/146‑1 (MQ=255)
                                                                                     gccccaccGGTGGGCGCGGGGGTCCCCGCTACACCGGCGCCGGGGTGAAGGGCCCCGGCTTTGCGGGTGGCGGCGGCCCCGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCt                                                                  <  1:4335655/143‑1 (MQ=255)
                                                                                            cGGTGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTGTGCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAAc                                                           <  2:495130/150‑1 (MQ=255)
                                                                                                TGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGaa                                                        >  2:5111761‑M1/1‑150 (MQ=255)
                                                                                                      CGGGCGTGCCCGCGAAACTGGCGCCGCGGTGCATGGCACCGGCTTTGCGGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                                                  >  2:1410964‑M1/1‑150 (MQ=255)
                                                                                                      CGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                                                  >  2:1436430‑M1/1‑150 (MQ=255)
                                                                                                        GGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAc                                                 <  1:4269004‑M1/150‑1 (MQ=255)
                                                                                                         GGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACt                                                <  2:1268773‑M1/150‑1 (MQ=255)
                                                                                                             gcCCCGCGAAACCGGCGCCGCGGGGCATGGCACCGGCTTGGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                            >  2:3172492‑M1/3‑150 (MQ=255)
                                                                                                               CCCGCGAACCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                          >  2:2104406‑M1/1‑150 (MQ=255)
                                                                                                               CCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                          <  2:6756648‑M1/150‑1 (MQ=255)
                                                                                                                ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcga                                         <  1:1164553‑M2/150‑107 (MQ=255)
                                                                                                                ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcga                                         <  1:5733706‑M2/150‑107 (MQ=255)
                                                                                                                CCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                                         <  1:2200654‑M1/150‑1 (MQ=255)
                                                                                                                 CGCTAAACTGGCGCCGCGTTGCTTGGCACCGGCTTTGCCGGTGTTGGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                        >  1:3259192‑M1/1‑150 (MQ=255)
                                                                                                                 CGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                        <  1:4227797‑M1/150‑1 (MQ=255)
                                                                                                                 CGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                        <  1:5802523‑M1/150‑1 (MQ=255)
                                                                                                                 CGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                        <  1:808570‑M1/150‑1 (MQ=255)
                                                                                                                  GCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                                       <  2:2340310‑M1/150‑1 (MQ=255)
                                                                                                                  GCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                                       <  1:696798‑M1/150‑1 (MQ=255)
                                                                                                                  GCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                                       <  2:3643855‑M1/150‑1 (MQ=255)
                                                                                                                  GCGAAACAGGCGCCGCGGTGCATGGCACCGGCTGTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGAGTGCCATTAGGCGGGCCTATCGCTACAGTGACGGCCCGTGTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAgac                                       <  2:3232412‑M1/150‑2 (MQ=255)
                                                                                                                  GCGAAACAGGCGCCGCGGTGCAGGGCACCGGCTTTGCGGTTGTGCGCGGGCCCGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                                       >  2:1699138‑M1/1‑150 (MQ=255)
                                                                                                                      AACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt                                   <  1:821766‑M1/150‑1 (MQ=255)
                                                                                                                        aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcgaactgcgg                                  <  1:2216556‑M2/150‑114 (MQ=255)
                                                                                                                         cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt                                 >  2:3441551/1‑150 (MQ=255)
                                                                                                                          cggCGCGGGGGTGCCGGGCACGGGCTTTGCCGGTGTTGGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                                <  2:4269004/149‑1 (MQ=39)
                                                                                                                          aGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTCTCGCTACCGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTGAGAACGCTATc                                >  1:16873/1‑150 (MQ=255)
                                                                                                                             cgccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                             <  1:1029967/150‑1 (MQ=255)
                                                                                                                                gccgGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                          >  2:686870/3‑150 (MQ=255)
                                                                                                                                     tGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                     >  1:4682627/1‑150 (MQ=255)
                                                                                                                                     tGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                     >  2:222903/1‑150 (MQ=255)
                                                                                                                                     ccgcgaacACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCACGCAGATCTg                     >  1:3897122/8‑150 (MQ=255)
                                                                                                                                       gcgaacACCGGCAAAGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                   >  2:4148689/6‑150 (MQ=255)
                                                                                                                                       cATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                   >  1:2336289/1‑150 (MQ=255)
                                                                                                                                            cACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTCTCAAGCAGATCTGAAACGTc              >  2:6321113/1‑150 (MQ=255)
                                                                                                                                             aCCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCATGCAGATCTGAAACGTcc             >  2:103372/1‑150 (MQ=255)
                                                                                                                                             aCCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc             >  1:4257857/1‑150 (MQ=255)
                                                                                                                                             aCCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc             >  2:3921840/1‑150 (MQ=255)
                                                                                                                                              ccGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc            >  1:4614054/1‑150 (MQ=255)
                                                                                                                                               cGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                    <  1:3020867/109‑1 (MQ=255)
                                                                                                                                               cGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                    >  2:3020867/1‑109 (MQ=255)
                                                                                                                                               cGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCg           >  2:1249171/1‑150 (MQ=255)
                                                                                                                                                ggCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGa          >  1:5579897/1‑150 (MQ=255)
                                                                                                                                                   tttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActt       >  1:2552918/1‑150 (MQ=255)
                                                                                                                                                      gCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtctt    >  2:106127/1‑150 (MQ=255)
                                                                                                                                                        cGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTcc  >  2:4713340/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
CGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCC  >  NC_002947/1504298‑1504596

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: