Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,452,283 T→G 13.0% A19A (GCT→GCG PP_5556 → bacteriophage lysis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,452,2830TG13.0% 704.1 / 19.7 316A19A (GCT→GCGPP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base T (149/114);  new base G (16/34);  total (167/149)
Fisher's exact test for biased strand distribution p-value = 1.82e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 323 total aligned reads displayed.

CGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCA  >  NC_002947/3452134‑3452432
                                                                                                                                                     |                                                                                                                                                     
cGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCt                                                                                                                                                       <  1:6051226/150‑1 (MQ=255)
   cATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggc                                                                                                                                                    <  1:3343175/150‑1 (MQ=255)
      ccTGGCGATCACCACACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggc                                                                                                                                                 <  2:927801/150‑1 (MQ=255)
       cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggcc                                                                                                                                                <  1:2184019/150‑2 (MQ=255)
       cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggcc                                                                                                                                                <  1:5706083/150‑2 (MQ=255)
        tGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggcca                                                                                                                                               <  1:866822/150‑3 (MQ=255)
        tGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggcca                                                                                                                                               <  1:8258440/150‑3 (MQ=255)
        tGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggcca                                                                                                                                               <  1:5421021/150‑3 (MQ=255)
           cGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGTGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCg                                                                                                                                            <  1:654594/150‑1 (MQ=255)
            gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCgg                                                                                                                                           >  1:6128037/1‑150 (MQ=255)
            gATCACCAAACGTATAATTGGTGGTAGAAAGGGTCTGGTGGATCGCAAGGCCCTCTACACGCGAGCGCTGGCGGTGCTGCTGTGCGGGCTTGGGCCACCCGCCTGGTCGCGGTGGCGCTGCTGCGGGCCGCCAGCGCTGGCGGCGCCCgg                                                                                                                                           <  2:2188911/150‑1 (MQ=255)
             aTCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCggg                                                                                                                                          >  2:1826053/1‑150 (MQ=255)
             aTCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCggg                                                                                                                                          >  2:5981365/1‑150 (MQ=255)
             aTCACAAAACGTATCAATGGTGGTAGCAAGGGTCTGGTGGATCGCATGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGTTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGGTGGCGGCGCCCggg                                                                                                                                          <  2:5148667/150‑1 (MQ=255)
                 ccAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGc                                                                                                                                      >  1:4235835/1‑150 (MQ=255)
                 ccAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGc                                                                                                                                      >  2:8236049/1‑150 (MQ=255)
                 ccAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGc                                                                                                                                      >  2:4745691/1‑150 (MQ=255)
                   aaaCGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCt                                                                                                                                    >  2:6908811/1‑150 (MQ=255)
                   aaaCGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCt                                                                                                                                    >  2:2953128/1‑150 (MQ=255)
                    aaCGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTg                                                                                                                                   >  2:7091202/1‑150 (MQ=255)
                    aaCGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGGCCGGGCAGccgg                                                                                                                                   >  1:2927353/1‑148 (MQ=255)
                     ccGTATCATTGGTGGTCGCAACGGTCTGGTGGATCGCACGGCCCTCTATAAGCGAGCGCTCGCGGTGCTGCTGTGAGCGCTGGGGCCACCCGCCTGGGCGCCGGGGCGCTGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTgg                                                                                                                                  <  2:3153987/149‑1 (MQ=255)
                       gTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGgt                                                                                                                                >  2:7847108/1‑150 (MQ=255)
                       gTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCCCTGGCGGCGCCCGGGCAGCCTgggg                                                                                                                                >  1:2533017/1‑149 (MQ=255)
                         aTCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGgtgt                                                                                                                              <  2:7335181/150‑1 (MQ=255)
                           cAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCCCCAGGGTGTgg                                                                                                                            >  1:3123867/1‑150 (MQ=255)
                             ttggtggTAGCACGGGTCTGGAGGATCGCACGGCCCTCTACCAGCGAGCGCTCGCGGGGCTGCTGTGAGCGCTTGGGCCACCCGCCTGGTCGCGGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCa                                                                                                                          <  1:1324087/149‑1 (MQ=255)
                             aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTATACAAGCGAGCGCTCGTGGTGCTGCAGTGTGCGCTTGGGCCACCCGCCTGGTCGCCTTGGCGCTGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCa                                                                                                                          <  1:3390777/150‑1 (MQ=255)
                              tggtggTTGCAAGGGTCTGGAGGCTCGCACGGCCCTCTACACGCGCGCGCTCGCGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCGGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAg                                                                                                                         <  2:3678378/150‑1 (MQ=255)
                               ggtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGGGTGGCAgg                                                                                                                        >  2:8493772/1‑150 (MQ=255)
                                gtggtAGCACCGGTCTGGAGGATCGCACGGCCCTCTCCACGCGCGCGCTCGAGGTGCTGCCGTGAGCGCTTGGGCCACCCGCCTGGTGGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGc                                                                                                                       <  2:6658754/150‑1 (MQ=255)
                                    tAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAAc                                                                                                                   >  2:1531814/1‑150 (MQ=255)
                                      gCACCGGTCTGGAGTATCGCAAGGCCCTCTACAAGCGCGCGCTCGCGGTGCTGCTGTGAGCGCTTGGGCCACCCGCGTGGTCGCCGTGGCGCTGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGc                                                                                                                 <  2:767357/150‑1 (MQ=255)
                                        aaCGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCt                                                                                                               >  2:3357303/1‑150 (MQ=255)
                                          gggTCTGGCGGCTCGCCCGGCCCCCTACCGGCGCGCGCTGGTGGGGCTGCGGGGCGGGGTGGGGCCCCCCGGGGGGGCGCGGGGGCGGGGGGGCGGGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCtat                                                                                                             <  1:5956674/149‑1 (MQ=255)
                                           ggTCTGGAGGATCGCACGGCCCTCTACAAGCGCGCGCTCGTGGTGCTGCAGTGCGCGCTTGGGCCACCGGCCGGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCtata                                                                                                            <  2:879651/150‑1 (MQ=255)
                                           ggTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCtata                                                                                                            >  2:4086080/1‑150 (MQ=255)
                                                 gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgcgc                                                                                                      >  2:6243822/1‑150 (MQ=255)
                                                 gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGCCAACGCCGATAAGgcgc                                                                                                      >  1:5881210/1‑150 (MQ=255)
                                                 gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCAATAAAgcgc                                                                                                      >  1:6169095/1‑150 (MQ=255)
                                                 gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCTCTAAGgcgc                                                                                                      >  1:3205449/1‑150 (MQ=255)
                                                 gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGCGAAGgcgc                                                                                                      >  2:6431553/1‑150 (MQ=255)
                                                    gaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGc                                                                                                   >  2:7596140/1‑150 (MQ=255)
                                                      tCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTc                                                                                                 >  2:6586071/1‑150 (MQ=255)
                                                       cGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCGATAAGGAGCCGCACg                                                                                                >  1:5141016/1‑150 (MQ=255)
                                                        gCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGc                                                                                               >  2:3468106/1‑150 (MQ=255)
                                                         cAAGGCCCTCTACACGCGAGCGCTCGTGGTGCGGCGGGGAGCGCTGGGGCCCCCCGGGTGGTGGCCGTGGGGGTGCTGGGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCt                                                                                              <  1:390698/150‑1 (MQ=255)
                                                         cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCt                                                                                              >  1:4458177/1‑150 (MQ=255)
                                                         cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCt                                                                                              >  2:7517683/1‑150 (MQ=255)
                                                          aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCAATAAGGCGCAGCTCGCTg                                                                                             >  1:7411993/1‑150 (MQ=255)
                                                                 tctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGc                                                                                      >  2:3552982/1‑150 (MQ=255)
                                                                   tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCGGGGGGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCg                                                                                    >  1:4934647/1‑150 (MQ=255)
                                                                   tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCg                                                                                    >  1:1242223/1‑150 (MQ=255)
                                                                     cacgcgagcgCTCGAGGGGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGAt                                                                                  <  2:5256466/147‑1 (MQ=255)
                                                                     cAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGAt                                                                                  >  2:7148583/1‑150 (MQ=255)
                                                                      aaGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATg                                                                                 >  2:4931985/1‑150 (MQ=255)
                                                                      aaGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCACCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACCCCGATAAGGCGCAGCTCGCGGTTCAGGCCGagg                                                                                 >  1:145844/1‑148 (MQ=255)
                                                                        gcgagcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGAc                                                                               >  2:611746/1‑150 (MQ=255)
                                                                           agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTc                                                                            <  2:568631/150‑1 (MQ=255)
                                                                                tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCAGCCGGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCGCGGATCAGGCCGATGACTTCGGCaa                                                                       >  1:6675957/1‑150 (MQ=255)
                                                                                tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCaa                                                                       >  2:1350513/1‑150 (MQ=255)
                                                                                tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAAAGCCAATAAGGAGCAGCTCGCTGATCAGGCCGAGGACTTCGGcca                                                                       >  1:7355971/1‑148 (MQ=255)
                                                                                  gcggTGCTGCAGTGAGCGCTGGGGCCACCCGCCGGGGCGCCGGGGCGCTGCGGGTGGCCCCCCGGGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGCTGACTTCGGCAAGc                                                                     <  1:7682176/148‑1 (MQ=255)
                                                                                      tgctgcTGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCTGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGtt                                                                 <  2:2250693/150‑1 (MQ=255)
                                                                                        ctgcCGTGCGCGCGTGGGCCCCCCGCGGGGTCGCGGTGGCGCTGCTGCTGGCCGCCAGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTgg                                                               <  1:900654/150‑1 (MQ=255)
                                                                                         tgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGc                                                              >  2:3915271/1‑150 (MQ=255)
                                                                                          gccgTGCGCGCGTGGGCCCCCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCg                                                             <  1:6370867/147‑1 (MQ=255)
                                                                                           cAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCgg                                                            >  1:3559967/1‑150 (MQ=255)
                                                                                            tgTGCGCGCTTGGGCCACCCGCCTGGTCGCCGGGGCGCTGCTGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGga                                                           <  1:1999972/149‑1 (MQ=255)
                                                                                              tGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGgaga                                                         <  2:2780869/150‑1 (MQ=255)
                                                                                              tGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGGGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCGGATCAGGCCGAGGACTTAGGCCAGCAGTTGGCGgaga                                                         >  1:411643/1‑150 (MQ=255)
                                                                                               gcgcgcTTGGGCCCCCCGCCTGGTCGCGTTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGaa                                                        <  1:8322026/148‑1 (MQ=255)
                                                                                               gAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGGGGTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCTAGCAGTTGGCGGAGaa                                                        <  2:1104356/150‑1 (MQ=255)
                                                                                               gAGCGCTTGGGCCACCCGCCTGGGCGCCGTGGCGCGGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGaa                                                        <  1:384909/150‑1 (MQ=255)
                                                                                                     gggggCCACCCGCCGGGTCGCGGGGGCGCTGCTGGGGGCCGCCAGCGCGGCCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcg                                                  <  1:6837194/148‑1 (MQ=255)
                                                                                                       gggCCACCCGCCTGGTCGCCGGGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcg                                                <  1:7243209/150‑1 (MQ=255)
                                                                                                        ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                                               <  2:4586968/150‑1 (MQ=255)
                                                                                                                 ccTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCt                                      >  2:2498019/1‑150 (MQ=255)
                                                                                                                      gcGCCGTGGCGCTGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGc                                 <  2:3138634/149‑1 (MQ=255)
                                                                                                                       cGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCg                                <  2:8043013/150‑1 (MQ=255)
                                                                                                                       cGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCg                                <  2:1706415/150‑1 (MQ=255)
                                                                                                                        gCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg                               >  2:7017674/1‑150 (MQ=255)
                                                                                                                         ccGTGGCGCTGCTGCTGGCCGCCTGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCCGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCggg                              <  2:1242223/150‑1 (MQ=255)
                                                                                                                         ccGTGGCGCTGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCggg                              <  2:2801189/150‑1 (MQ=255)
                                                                                                                          cGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGga                             <  2:3997061/150‑1 (MQ=255)
                                                                                                                             ggCGCTGCTGGTGGCCGCCTGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa                          <  1:1067376/150‑1 (MQ=255)
                                                                                                                             ggCGCTGCTGCTGGCCGCCTGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa                          <  1:8036239/150‑1 (MQ=255)
                                                                                                                               cgcTGCTGCTGGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAgg                        <  1:8493772/150‑1 (MQ=255)
                                                                                                                                gctgctgctGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGc                       >  2:4249201/1‑150 (MQ=255)
                                                                                                                                  ggctgggggCCGCCTGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTg                     <  1:5253790/143‑1 (MQ=255)
                                                                                                                                   gctgcggGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgc                    <  1:4627546/150‑1 (MQ=255)
                                                                                                                                   cccgggcgccgccAGCGCTGGCGGCCCCCCGGCACCCTGGGTGTGGCAGGCCAACCCCGAAAACGCGCCGCGCGCGGATCAGGCCGATGACATCCGCAAGCAGTTGGCGGCGAACGATCGCGCCCACCGGCTGGAGCGGGAGAAGGCAgc                    >  2:5597753/8‑150 (MQ=255)
                                                                                                                                      accagtggccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagc                 >  2:6050014/8‑150 (MQ=255)
                                                                                                                                        tgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcag               >  2:1062932/1‑150 (MQ=255)
                                                                                                                                        aggCCGCCAGAGCTCGCGGCCCCCGCGCAGCCTGGGTGAGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCCATGACTTCCCCAAGCAGTTCGCGGAGAAGGATCGCGCGCACGGTCTGGCGCGGGAGAAGGCTgcagcac               >  2:2594893/2‑149 (MQ=255)
                                                                                                                                         cgccgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcagc              >  1:7821955/2‑150 (MQ=255)
                                                                                                                                            cgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCagcagcagc           >  2:1453541/1‑150 (MQ=255)
                                                                                                                                               ccgcgcTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCTATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCt        <  1:278680/148‑1 (MQ=255)
                                                                                                                                               cAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCt        >  1:1650333/1‑150 (MQ=255)
                                                                                                                                                     tGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCa  >  2:7424680/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCA  >  NC_002947/3452134‑3452432

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: