Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,163,456 C→G 7.4% A158G (GCG→GGG)  PP_3661 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,163,4560CG 7.4% 454.0 / 3.0 177A158G (GCG→GGG) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base C (78/84);  new base G (11/4);  total (89/88)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

TATCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTG  >  NC_002947/4163307‑4163605
                                                                                                                                                     |                                                                                                                                                      
  tCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGgcgc                                                                                                                                                      >  1:1115549/1‑150 (MQ=255)
    gtgtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGgcgcgg                                                                                                                                                    >  1:4015352/1‑148 (MQ=255)
    gtgtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTg                                                                                                                                                    >  1:321527/1‑150 (MQ=255)
       tCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTgcgc                                                                                                                                                 <  2:1844748/150‑1 (MQ=255)
           tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcg                                                                                                                                             <  1:1870784/150‑1 (MQ=255)
             cTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑cc                                                                                                                                          <  1:6053875/150‑1 (MQ=255)
              tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑ccc                                                                                                                                         >  1:4328757/1‑150 (MQ=255)
               tgTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑ccga                                                                                                                                        <  1:515725/149‑3 (MQ=255)
                   gTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCTGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCAt                                                                                                                                    >  2:1092979/1‑150 (MQ=255)
                    tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCaga                                                                                                                                   >  2:1400742/1‑148 (MQ=255)
                      aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGCGGAGCG‑cccccccacaa                                                                                                                                 >  1:1949942/1‑145 (MQ=255)
                      aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATaaa                                                                                                                                 >  2:1133135/1‑150 (MQ=255)
                      aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGAAACCGCCTTGCTGACCGCCGCCGGAAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGAAGCG‑CCCCCCAaaaa                                                                                                                                 >  1:1115094/1‑150 (MQ=255)
                       caaGTACCCTGCTGATCGGCCTGTGGTGTGGAATGAAATTGTTCAAGCTCGACCGCGATACCGCATTGCTGACCGCCGCCGGCAGCGCAATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATaaaa                                                                                                                                <  1:4570280/148‑1 (MQ=255)
                         aGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAgc                                                                                                                              >  2:6557356/1‑150 (MQ=255)
                         aGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAGAAAAgc                                                                                                                              >  2:5166585/1‑150 (MQ=255)
                           tACCGTGCTGTTCGGCCTGTGGTGTGGCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAgcgc                                                                                                                            <  1:5138199/150‑1 (MQ=255)
                             cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCAAAAAGCCCCa                                                                                                                          >  1:4122624/1‑150 (MQ=255)
                              ccTGCTGATGGGCCTGTGGTGTGGCATTAAAGTGTTCAAGCTCGACCGCGATACCGCGTTGCTGACCGCCGCGGGCAGCGCCATATGCGGCGCCGCCGCCGTGCGGGAATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCAt                                                                                                                         <  2:783689/150‑1 (MQ=255)
                              ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCCTCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCCGCGCAGCC‑CCCCCCCCAAAAACGccca                                                                                                                         >  1:2942985/1‑148 (MQ=255)
                                 gctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCAAAAACCGCCAAGGc                                                                                                                      >  1:6471096/1‑150 (MQ=255)
                                  ctgATCGGCGTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCGGGAAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCACGGCa                                                                                                                     <  2:5821178/150‑1 (MQ=255)
                                  ctgATCGGCCTGTGGTGTGGCATTAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCa                                                                                                                     <  2:3867176/150‑1 (MQ=255)
                                  ctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCC‑CCCCCCAAAAAAACGCCAAGGCa                                                                                                                     >  2:2358658/1‑150 (MQ=255)
                                    gATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCCCG‑CCCCCCAAAAAAATTCCAAGGCAga                                                                                                                   >  1:5586789/1‑149 (MQ=255)
                                     ttgggCCTGTGGTGTGGCATGAATGTGTTCATGCTCGACCGCGATACCGCCTTGCTGACGGCCGCGGGCTGCGCCTTCTGGGGCGCCGCCGCCGTGGTGGCTTTCGAAGCGGGGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTc                                                                                                                  <  2:4931340/147‑1 (MQ=255)
                                     aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTTGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTc                                                                                                                  >  2:2738046/1‑150 (MQ=255)
                                     aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCAAGGCAGTc                                                                                                                  >  2:1867509/1‑150 (MQ=255)
                                        ggCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGa                                                                                                               <  2:2094504/150‑1 (MQ=255)
                                         gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAGAAAAACGCCCAGGCAGTCggg                                                                                                               >  1:3228886/1‑149 (MQ=255)
                                         gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                              <  2:4545950/150‑1 (MQ=255)
                                         gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCCAAAAAACCCCAAGGCactacgaa                                                                                                              >  1:1474689/1‑143 (MQ=255)
                                         gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCC‑CCCCCCAAAAAATCGCCCTGGCAGTCGGaa                                                                                                              >  2:4847263/1‑150 (MQ=255)
                                          ccTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCC‑CCCCCCAGAAAAGCGCCAAGGCAGTCGGGAg                                                                                                             >  1:2587132/1‑150 (MQ=255)
                                            tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCg                                                                                                           >  2:1522279/1‑150 (MQ=255)
                                              tggtgtggCATGCATGTGTTCAAGCTCCCCCGCGATACGGCATTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCAGCAGTGCTGGCATTAGAATCGGCGCTGCGCTGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                         <  2:3007826/150‑1 (MQ=255)
                                                   tggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCG‑CCCCCCCAAAAAGCGCCCGGGCCGGCGGGAGCGTAGTGCt                                                                                                    >  1:4203337/1‑150 (MQ=255)
                                                   tggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCC‑CCCCCCCAAAAAGCCCCAAGGCAGTCGGAAGCGGAGTGCt                                                                                                    >  1:2738591/1‑150 (MQ=255)
                                                   tggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCACCGCCGCCGTGCTGGCATTCGAATCGGGGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTTGTGCt                                                                                                    >  2:4126131/1‑150 (MQ=255)
                                                      cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCAAAAAACGCCCAGGCAGGAGGAAGCGTAGTGCAGtt                                                                                                 >  1:5032072/1‑150 (MQ=255)
                                                       aggAAAGTGTTCACGCTCGACCGCGATACCGCATTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTc                                                                                                <  2:4702608/148‑1 (MQ=255)
                                                       aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGGAGTGCGGTTc                                                                                                >  1:6707750/1‑150 (MQ=255)
                                                       aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCC‑CCCCCCCCAAAAGCCCCCAGGCAGTCGGAAGCGGAGTGCGGTTc                                                                                                >  1:363164/1‑150 (MQ=255)
                                                          aaaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGt                                                                                             <  2:6416571/150‑1 (MQ=255)
                                                           aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTa                                                                                            <  2:2840236/150‑1 (MQ=255)
                                                           aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCCAAAAAGCGCCAATGCAGTCGGCAGCGTAGTGCGGTTAGGtc                                                                                            >  2:5469794/1‑149 (MQ=255)
                                                            aGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAc                                                                                           >  2:3465547/1‑150 (MQ=255)
                                                            aGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCCTAAAAGCGCCATGGCAGTCGGAAGCGGAGTGCTGTTCGGtcc                                                                                           >  2:4314373/1‑148 (MQ=255)
                                                              tgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGATTCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCt                                                                                         <  2:688801/150‑1 (MQ=255)
                                                              tgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCt                                                                                         <  2:6211158/150‑1 (MQ=255)
                                                                   aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCCTTCGAATCGGCGCGGCGCAGCG‑CCCCCCCGAAAAGCGCCATGGCCGTCGGGAGCGGAGTGCTGTTCGGTACCTTGTCg                                                                                    >  2:911907/1‑150 (MQ=255)
                                                                     gcgcgcccGCGATCCCGCCTGGCCGACCGCCGCCGGCCGCGCCCGCTGCGGCGCGGCCGCGGGGCTGGCATTCGACTGGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                  <  1:172270/144‑1 (MQ=255)
                                                                     gCTCGACGGCGATCCGGCCTGGCTGCCCGCCGCGGGCTGCGCCATCGGGGGCGCCGCCGCGGTGGGGGCATTCTCCGCGGCGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                  <  1:3458525/150‑1 (MQ=255)
                                                                     gCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                  >  2:6844703/1‑150 (MQ=255)
                                                                       tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                >  2:6653835/1‑150 (MQ=255)
                                                                       gcGACCGCTATACCGCCTTGCGGACCGCCGCCGGCCGCGCCTTCTGCGGCGCCGCCGCCGTGCTGGCTTTCGAAGGGGCGCTGCGCTGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                <  2:1886436/149‑1 (MQ=255)
                                                                            cgcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCAATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCtt                                                                           <  2:173453/150‑1 (MQ=255)
                                                                              cgatcgcGCCGTGCGGACCGCCGCCGGCTGCGCCATCGGGGGCGCCGCGGCCGTGGGGGAATTCGATTGGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTa                                                                         <  1:4763876/144‑1 (MQ=255)
                                                                               gATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAc                                                                        >  2:3273947/1‑150 (MQ=255)
                                                                                     gcgtggCTGCCCGCCGCCGGCAGCGCCATGGGGGGCGCGGCCGCGGTGCTGGAATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTg                                                                  <  1:5914906/145‑1 (MQ=255)
                                                                                     gCCTTGCTGCCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTg                                                                  <  2:3310190/150‑1 (MQ=255)
                                                                                     gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATGGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTg                                                                  <  2:4884386/150‑1 (MQ=255)
                                                                                         gggggcccgccggcgGCCGCGCCTTGGGGGGCGCGGCCGGGGTGGGGGCTTTCGCCGGGGCGGGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTGGTCGATGTTCCTTTACCCACTGGCGa                                                              <  1:5476015/144‑1 (MQ=255)
                                                                                          gCTGCCCGCCGCCGGCTGCGCCATCGGGGGCGCCGCCGCGGTGGTGGCATTCGAAGCGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGAt                                                             <  1:4295955/150‑1 (MQ=255)
                                                                                           cTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCAAGGCAGTCGGAATCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATc                                                            >  2:6276916/1‑150 (MQ=255)
                                                                                            tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTTCTGGAATACGAAGCGGCGCGGCCCAGCC‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCa                                                           <  1:1079080/150‑1 (MQ=255)
                                                                                               ccgccgccgGCCGCGCCATCGGGGGCGCCGCCGCCGTGGTGGCATTCGACGCGGCGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATg                                                        <  1:5410054/150‑1 (MQ=255)
                                                                                               ccgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCCTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATg                                                        >  2:5301104/1‑150 (MQ=255)
                                                                                                cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCCGAAAAGCGCCAGGGCAGGCGGAGGCGTAGTGGTGGTCGGTACCATGTAGATGTTACTTTAACCACAGGCGATAAAAGc                                                       >  1:2946680/1‑150 (MQ=255)
                                                                                                  ccgcggGCAGCGCCTTGTGGGGCGCCGCCGCGGTGCGGGATTTCGAAGGGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCg                                                     <  1:4493635/150‑1 (MQ=255)
                                                                                                  ccgccgGCAGCGCCATGTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCg                                                     <  2:4427882/150‑1 (MQ=255)
                                                                                                    gccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggc                                                   <  2:3561553/150‑1 (MQ=255)
                                                                                                       ggCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCAAAAAGCGCCCTGGCAGTCGGAAGCGTTGTGCTGGTCGGTACCTTGTCGATGTTCCATTAACCACAGGCGATCAATTCCggcagg                                                >  1:2238739/1‑150 (MQ=255)
                                                                                                              ccATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAAAAAAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGGTCGGTACCGTGTCGATGTTACTTTAACCACAGGCGATCAATTCCGGCGGGCGGCAta                                         >  1:3952823/1‑149 (MQ=255)
                                                                                                               cATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTACTTTACCCACAGGCGATCAAATCCGGCTGGCTGCATCt                                        >  2:5556023/1‑150 (MQ=255)
                                                                                                               cATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCCTAAAAACGCCAAGGCCGTCGGGAGCGGAGTGCTGTTCGGGGCCGTGTCGGTGTTCCTTTAACCACAGGCGATCCATGCCGGCGGGCGGCATCt                                        >  1:5592300/1‑150 (MQ=255)
                                                                                                               cATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCCCAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGCTGTTGCTTTACCCACAGGCGCTCAAAGCCGGCTGGGGGCATcg                                        >  1:2126803/1‑149 (MQ=255)
                                                                                                                 tCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCCTAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGGACCGTGTCGATGTTCCTTTAACCACTGGCGATCAAAGCCGGCTGGCTGCAACTgg                                      >  2:2023856/1‑150 (MQ=255)
                                                                                                                 tCTGCGGCGCCGCCGCCGGGCTGGAATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTgg                                      <  1:3729606/150‑1 (MQ=255)
                                                                                                                      ggcgCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAAGCCGGCTGGCTGCATCTGGATAcc                                 >  2:2924411/1‑150 (MQ=255)
                                                                                                                       gcgcCGCGGCGGTGCTGGAATTCGAATCGGGGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCa                                <  1:1666988/150‑1 (MQ=255)
                                                                                                                        cgccgccgcCGTGGTGGATTTCTATGCGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAt                               <  1:5164805/150‑1 (MQ=255)
                                                                                                                         gccgccgccGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATg                              <  2:1574303/150‑1 (MQ=255)
                                                                                                                         gccgccgccGTGCTGGCATTCGAAGCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATg                              <  2:1463334/150‑1 (MQ=255)
                                                                                                                          ccgccgccgTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                             <  2:5704747/150‑1 (MQ=255)
                                                                                                                          ccgccgccgTGCTGGAATTCTAAGGGGCGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                             <  2:3046621/150‑1 (MQ=255)
                                                                                                                                    gCTGGCATTCGAAGCGGCGCTGCGCAGCG‑CCCCCCATAAACGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgct                   <  1:6202421/150‑1 (MQ=255)
                                                                                                                                           ttCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATTGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATTTTCCTTTAACCACAGGCGATCAATGCCGGCTGGCTGCATCTGGATAACATGGGCGCCGGGctgctgctg            >  2:1178064/1‑150 (MQ=255)
                                                                                                                                           ttCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCCCCAAGGCAGTCGGAAGCGTAGTTCTGTTAGGTACCTTGTCGATGTTACTTTAAACACAGGCGATCAATGCCCGCTGGCTGCATCTGGATAACAAGGGCGCCGGGctgctgctg            >  1:3103654/1‑150 (MQ=255)
                                                                                                                                             cGAAGCGGCGCTGCGCCCCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTggg          <  1:5543592/150‑1 (MQ=255)
                                                                                                                                              gAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc         <  2:5637987/150‑1 (MQ=255)
                                                                                                                                              gAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAACTTGTCGAAGTTCCTTTAACCACAGGCGATCAATTCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc         >  2:4002882/1‑150 (MQ=255)
                                                                                                                                              gAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTAGATGTTCCTTTACCCACAGGCGATCAATGCCCGCTGGCTGCATCTGGATACCATGGGCGCCGGGCAGCTGCTGggc         >  2:6155848/1‑150 (MQ=255)
                                                                                                                                               aaTCGGCGCTGCGCAGCGCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc         <  2:3228886/150‑1 (MQ=255)
                                                                                                                                                  cggcgcggcgCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGGACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa     <  2:5159499/150‑1 (MQ=255)
                                                                                                                                                  cGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCCCTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa     <  2:2065215/150‑1 (MQ=255)
                                                                                                                                                    gcgcTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt   >  2:6409533/1‑150 (MQ=255)
                                                                                                                                                    gcgcTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt   >  1:3849332/1‑150 (MQ=255)
                                                                                                                                                    gcgcTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCAATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt   >  2:1243439/1‑150 (MQ=255)
                                                                                                                                                     cgcTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTg  >  2:2150735/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TATCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTG  >  NC_002947/4163307‑4163605

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: ‑