Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,163,464 A→C 20.8% S161R (AGC→CGC)  PP_3661 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,163,4640AC20.8% 282.3 / 21.8 179S161R (AGC→CGC) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base A (67/64);  new base C (25/17);  total (93/86)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

CGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACG  >  NC_002947/4163315‑4163611
                                                                                                                                                     |                                                                                                                                                    
  gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                                    <  1:379219/150‑1 (MQ=255)
   tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcg                                                                                                                                                   <  1:1870784/150‑1 (MQ=255)
     cTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑cc                                                                                                                                                <  1:6053875/150‑1 (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)
              aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGAATCGGCGCGGCGCCGCG‑cccccccgcaa                                                                                                                                       >  2:4201963/1‑145 (MQ=255)
              aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGAAACCGCCTTGCTGACCGCCGCCGGAAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGAAGCG‑CCCCCCAaaaa                                                                                                                                       >  1:1115094/1‑150 (MQ=255)
               caaGTACCCTGCTGATCGGCCTGTGGTGTGGAATGAAATTGTTCAAGCTCGACCGCGATACCGCATTGCTGACCGCCGCCGGCAGCGCAATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATaaaa                                                                                                                                      <  1:4570280/148‑1 (MQ=255)
                cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAGAAAAg                                                                                                                                     >  1:2757823/1‑150 (MQ=255)
                 aGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAGAAAAgc                                                                                                                                    >  2:5166585/1‑150 (MQ=255)
                   tACCGTGCTGTTCGGCCTGTGGTGTGGCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAgcgc                                                                                                                                  <  1:5138199/150‑1 (MQ=255)
                     cccTGCTTATCGGCCTGTGGTGTGTCATGAAAGTGTTCAAGCTCGACGGCTATACCGCCTTGCTGACCGCCGCCGGAAGCGCAATCTGCGGCGCCGCCGCGGTGCTGGAATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCa                                                                                                                                <  2:3096800/150‑1 (MQ=255)
                     cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCAAAAAGCCCCa                                                                                                                                >  1:4122624/1‑150 (MQ=255)
                      ccTGCTGATGGGCCTGTGGTGTGGCATTAAAGTGTTCAAGCTCGACCGCGATACCGCGTTGCTGACCGCCGCGGGCAGCGCCATATGCGGCGCCGCCGCCGTGCGGGAATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCAt                                                                                                                               <  2:783689/150‑1 (MQ=255)
                         gctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCAAAAACCGCCAAGGc                                                                                                                            >  1:6471096/1‑150 (MQ=255)
                          ctgATCGGCCTGTGGTGTGGCATTAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCa                                                                                                                           <  2:3867176/150‑1 (MQ=255)
                          ctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCC‑CCCCCCAAAAAAACGCCAAGGCa                                                                                                                           >  2:2358658/1‑150 (MQ=255)
                            gATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCAAGGCCGt                                                                                                                         >  1:1676693/1‑150 (MQ=255)
                             agcGGCCTGTGGTGTGGCATGAATGTGTTCACGCTCGACCGCGAGACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTc                                                                                                                        <  2:4931322/148‑1 (MQ=255)
                             aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTTGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTc                                                                                                                        >  2:2738046/1‑150 (MQ=255)
                              tCGGCCTGTGGTGTGGCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGAAGCGCCATGTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCg                                                                                                                       <  2:3304139/150‑1 (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)
                                 gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCC‑CCCCCCAAAAAATCGCCCTGGCAGTCGGaa                                                                                                                    >  2:4847263/1‑150 (MQ=255)
                                  cgtgtggtgtgtCATGAAATTGTTCAAGCTCGACCGCGATCCCGCCTTGCTGCCCGCCGCGGGCCGCGCCTTCGGCGGCGCCGCCGCCGTGCTGGCTTTCGAATGGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAg                                                                                                                   <  1:3409996/148‑1 (MQ=255)
                                  ccTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCC‑CCCCCCAGAAAAGCGCCAAGGCAGTCGGGAg                                                                                                                   >  1:2587132/1‑150 (MQ=255)
                                    tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCCAAAAGCGCCCAGGCAGTCGGAAGCg                                                                                                                 >  1:4380776/1‑150 (MQ=255)
                                      tggtgtggCTTGAATGTGTTCAAGCTCGACCGCTATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                               <  2:6487138/150‑1 (MQ=255)
                                      tggtgtggCATGCATGTGTTCAAGCTCCCCCGCGATACGGCATTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCAGCAGTGCTGGCATTAGAATCGGCGCTGCGCTGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                               <  2:3007826/150‑1 (MQ=255)
                                      tggtgtggCATGAAAGTGTTCACGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCCGTGGTGGCTTTCGATTCGGGGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                               <  2:1874474/150‑1 (MQ=255)
                                      tggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCC‑CCCCCCCTAAAACCGCCCAGGCCGTCGGAcgcgga                                                                                                               >  2:5476015/1‑148 (MQ=255)
                                        gtgtggCATTAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATATGCGGCGCCGCCGCGGTGCTGGCATTCGAATCGGCGGTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGt                                                                                                             <  2:2067976/150‑1 (MQ=255)
                                           tggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCG‑CCCCCCCAAAAAGCGCCCGGGCCGGCGGGAGCGTAGTGCt                                                                                                          >  1:4203337/1‑150 (MQ=255)
                                           tggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCACCGCCGCCGTGCTGGCATTCGAATCGGGGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTTGTGCt                                                                                                          >  2:4126131/1‑150 (MQ=255)
                                              cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCAAAAAACGCCCAGGCAGGAGGAAGCGTAGTGCAGtt                                                                                                       >  1:5032072/1‑150 (MQ=255)
                                               aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGGAGTGCGGTTc                                                                                                      >  1:6707750/1‑150 (MQ=255)
                                                 gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCC‑CCCCCCCGAAAAGCGCCCAGGCACTCGGAAGCGGAGTGCTGGTCgg                                                                                                    >  2:6602113/1‑150 (MQ=255)
                                                   atgttttCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCAATCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTa                                                                                                  <  2:3040254/145‑1 (MQ=255)
                                                   aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTa                                                                                                  <  2:2840236/150‑1 (MQ=255)
                                                    aGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAc                                                                                                 >  2:3465547/1‑150 (MQ=255)
                                                      tgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGATTCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCt                                                                                               <  2:688801/150‑1 (MQ=255)
                                                      tgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCt                                                                                               <  2:6211158/150‑1 (MQ=255)
                                                          aacgCTCGACCGCGATACCGCTTTGCTGCCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                           <  1:3626794/147‑1 (MQ=255)
                                                           aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCCTTCGAATCGGCGCGGCGCAGCG‑CCCCCCCGAAAAGCGCCATGGCCGTCGGGAGCGGAGTGCTGTTCGGTACCTTGTCg                                                                                          >  2:911907/1‑150 (MQ=255)
                                                            aGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCG‑CCCCCCCAAAAAACCCCCAGGGAGTCGGGAGCCGAGTGCTGTTGGGTGCCgtgtggt                                                                                         >  1:3697895/1‑147 (MQ=255)
                                                             gCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                        >  2:6844703/1‑150 (MQ=255)
                                                               tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                      >  2:6653835/1‑150 (MQ=255)
                                                                   ccGCTATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCt                                                                                  <  2:2017986/150‑1 (MQ=255)
                                                                    cgcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCAATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCtt                                                                                 <  2:173453/150‑1 (MQ=255)
                                                                     gcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCttt                                                                                >  2:1864782/1‑150 (MQ=255)
                                                                      cgatcgcGCCGTGCGGACCGCCGCCGGCTGCGCCATCGGGGGCGCCGCGGCCGTGGGGGAATTCGATTGGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTa                                                                               <  1:4763876/144‑1 (MQ=255)
                                                                       gATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CCCCCCCTAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCATTaa                                                                              >  2:3273890/1‑149 (MQ=255)
                                                                          aCCGCCTTGGTGCCCGCCGCCGGCCGCGCCCCCGGGGGCGCCGCCGCGGGGCTGGAATTCTACTCGGCGGGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCa                                                                           <  1:1855145/150‑1 (MQ=255)
                                                                             gcgtggCTGCCCGCCGCCGGCAGCGCCATGGGGGGCGCGGCCGCGGTGCTGGAATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTg                                                                        <  1:5914906/145‑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)
                                                                                    tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGGGCCGCG‑CCCCCCAAAAAAAAGCCAAGGCAGTAGGAAGCGTAGTGCAGTTCGGTACCATTTAGATGTTGCTTTAACCACAGGCGagaa                                                                 >  1:5558273/1‑147 (MQ=255)
                                                                                       ccgccgccgGCCGCGCCATCGGGGGCGCCGCCGCCGTGGTGGCATTCGACGCGGCGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATg                                                              <  1:5410054/150‑1 (MQ=255)
                                                                                       ccgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCCTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATg                                                              >  2:5301104/1‑150 (MQ=255)
                                                                                        cgccgcggGCCGCGCCTTGTGGGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGc                                                             <  1:3967593/150‑1 (MQ=255)
                                                                                        cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCCGAAAAGCGCCAGGGCAGGCGGAGGCGTAGTGGTGGTCGGTACCATGTAGATGTTACTTTAACCACAGGCGATAAAAGc                                                             >  1:2946680/1‑150 (MQ=255)
                                                                                          ccgcggGCAGCGCCTTGTGGGGCGCCGCCGCGGTGCGGGATTTCGAAGGGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCg                                                           <  1:4493635/150‑1 (MQ=255)
                                                                                          ccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCg                                                           >  2:6193092/1‑150 (MQ=255)
                                                                                            gccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggc                                                         <  2:3561553/150‑1 (MQ=255)
                                                                                               ggCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCAAAAAGCGCCCTGGCAGTCGGAAGCGTTGTGCTGGTCGGTACCTTGTCGATGTTCCATTAACCACAGGCGATCAATTCCggcagg                                                      >  1:2238739/1‑150 (MQ=255)
                                                                                                gCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGACGCGGCGCTGCGCAGCG‑CCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctggc                                                     <  2:2038975/150‑1 (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)
                                                                                                                cgccgccgcCGTGGTGGATTTCTATGCGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAt                                     <  1:5164805/150‑1 (MQ=255)
                                                                                                                 gccgccgccGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCagg                                    >  2:1479665/1‑148 (MQ=255)
                                                                                                                 gccgccgccGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATg                                    <  2:1574303/150‑1 (MQ=255)
                                                                                                                  ccgccgccgTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                                   <  2:5704747/150‑1 (MQ=255)
                                                                                                                       gccgTGCTGGCATGCGACTCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGcc                              <  2:5981770/150‑1 (MQ=255)
                                                                                                                         ggTGCTGGCATTCGAAGCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCgg                            <  2:1211486/149‑1 (MQ=255)
                                                                                                                            gCTGGCATTCGAAGCGGCGCTGCGCAGCG‑CCCCCCATAAACGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgct                         <  1:6202421/150‑1 (MQ=255)
                                                                                                                            gCTGGCATGCGACGCGGCGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgct                         <  1:272812/150‑1 (MQ=255)
                                                                                                                              ttgatttCGACGCGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgc                       <  1:1523432/145‑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)
                                                                                                                                      gAAGCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc               <  2:3410895/150‑1 (MQ=255)
                                                                                                                                      gAAGCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc               <  2:4328757/150‑1 (MQ=255)
                                                                                                                                       aaTCGGCGCTGCGCAGCG‑CCCCCCAAAAAAACCCCAAGGCAGTCGGAAGCGTAGTGCTGTTAGGTACCATGTCGATGTTACATTACACAAACGCGAAAAAAGACGGCTGGCTGCATCAGGATAACAATGGCGCCGGGCTGCTGCTGggcg              >  1:1379457/1‑150 (MQ=255)
                                                                                                                                         tCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcggc            >  1:513964/1‑150 (MQ=255)
                                                                                                                                          cggcgcggcgCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGGACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa           <  2:5159499/150‑1 (MQ=255)
                                                                                                                                            gcgcTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt         >  2:6409642/1‑150 (MQ=255)
                                                                                                                                             cgcTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTg        >  2:2150735/1‑150 (MQ=255)
                                                                                                                                              gcTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt       >  1:6552324/1‑150 (MQ=255)
                                                                                                                                               cTGCGCAGCG‑CCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTc      <  2:5032072/150‑1 (MQ=255)
                                                                                                                                                ggcgccgcg‑ccCCCCATACCCGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTcc     <  1:385331/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACG  >  NC_002947/4163315‑4163611

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