Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,163,464 A→C 16.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,4640AC16.8% 344.0 / 21.3 187S161R (AGC→CGC) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base A (78/71);  new base C (19/16);  total (97/90)
Fisher's exact test for biased strand distribution p-value = 8.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

GGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGA  >  NC_002947/4163316‑4163612
                                                                                                                                                    |                                                                                                                                                     
ggTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAg                                                                                                                                                      <  1:3348649/150‑1 (MQ=255)
  tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcg                                                                                                                                                    >  2:3629161/1‑150 (MQ=255)
  tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcc                                                                                                                                                    >  2:4774286/1‑149 (MQ=255)
     tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGTGCGCAGCGccc                                                                                                                                                 >  1:267916/1‑150 (MQ=255)
     tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccc                                                                                                                                                 <  1:523342/150‑1 (MQ=255)
     tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccc                                                                                                                                                 <  2:6497709/150‑1 (MQ=255)
     tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccc                                                                                                                                                 <  2:1603436/150‑1 (MQ=255)
     tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGGTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGACTCGGCGCGGCGCCGCGccc                                                                                                                                                 >  2:6840431/1‑150 (MQ=255)
       gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGGTGCGCAGCGccccc                                                                                                                                               <  2:3505479/150‑1 (MQ=255)
       gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccccc                                                                                                                                               >  2:850635/1‑150 (MQ=255)
        tAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccc                                                                                                                                              <  2:4944778/150‑1 (MQ=255)
           tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTGGGGAGCGcccccccaa                                                                                                                                           >  1:2318050/1‑147 (MQ=255)
             aCCAGTACCCTGCTGTTGGGCCTGTGGGGTGGCATGAATGTGTTCAAGCTCGACGGCGATACCGCGTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTGGATTGGGGGCTGGGCAGCGCCCCCCATaaa                                                                                                                                         <  2:2901717/150‑1 (MQ=255)
             aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTGCGCACCGcccccccgaaa                                                                                                                                         >  1:4498268/1‑145 (MQ=255)
             aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccccaaaa                                                                                                                                         >  1:3814296/1‑145 (MQ=255)
             aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAaaaa                                                                                                                                         >  1:3635985/1‑150 (MQ=255)
              ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATTAAAGTGTTCACGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATGTGCGGCGCCGCGGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATaaaa                                                                                                                                        <  2:7376860/150‑1 (MQ=255)
              ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCTGCGCCCCCCATaaaa                                                                                                                                        <  2:5731604/150‑1 (MQ=255)
                cgTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCTATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCTGCGCCCCCCATAAAAgc                                                                                                                                      <  2:2326689/149‑1 (MQ=255)
                aGTACCCTGCTGTTCGGCCTGTGGTGTGGAATTAATTTTTTCATGCTCTACCGCTATACCGCCTTGCTGACGGCCGCCGGCAGCGCCATCTGCGGCGCGGCCGCCGTGCTGGCATTCGAATCGGCGGTGCGCAGCGCCCCCCATAAAAgc                                                                                                                                      <  2:1380873/150‑1 (MQ=255)
                aGTACCCTGCTGATCGGCCTGTGGTGTGGAATGAAAGTGTTCAAGCTCGACCGCTATACCGCCTTGCTGACGGCCGCCGGCAGCGCATTCTGCGGCGCCGCCGCCGTGCTGGATTTCGATTCGGCGCTGCGCAGCGCCCCCCATAAAAgc                                                                                                                                      <  2:6995330/150‑1 (MQ=255)
                     ccTGCTGATCGGCCTGTGGTGTGGAATGAAAGTGTTCAAGCTCGCCCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCAt                                                                                                                                <  2:5581521/150‑1 (MQ=255)
                      ctgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCATGCTCGACCGCGATACGGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATg                                                                                                                               <  2:3737344/150‑1 (MQ=255)
                      ctgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCCGCGCAGACCCCCCCAAAAAAAAG‑CCaag                                                                                                                               >  1:2924944/1‑148 (MQ=255)
                       tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATgg                                                                                                                              >  2:5372563/1‑150 (MQ=255)
                          tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCAGGGccg                                                                                                                           >  2:6179152/1‑148 (MQ=255)
                           gATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAAGCG‑CCCAGGCAGt                                                                                                                          >  1:4298301/1‑150 (MQ=255)
                            aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCACAAAAACG‑CCAAGGCAGTc                                                                                                                         >  1:456809/1‑150 (MQ=255)
                                gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGaa                                                                                                                     <  2:459008/150‑1 (MQ=255)
                                gCCTGGGGTGTGGCATGAATGTGTTCATGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGATTCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGaa                                                                                                                     <  2:2460703/150‑1 (MQ=255)
                                  cTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGc                                                                                                                   <  2:267916/150‑1 (MQ=255)
                                     tggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCC‑CCCGGGGCGTCGGCAGCGtc                                                                                                                >  1:4592167/1‑149 (MQ=255)
                                       gtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCAAGGCAGTCGGGAGCGTAGt                                                                                                              >  2:2177358/1‑150 (MQ=255)
                                        tgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTg                                                                                                             <  2:4925751/150‑1 (MQ=255)
                                          tggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCCCCCCCCCCCAAAAAGCC‑CCCAGGCCGACGGAAGCGGAGTggg                                                                                                           >  1:5739799/1‑148 (MQ=255)
                                           ggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCACAAAGAG‑CCAAGGCAGGCCGGAGCGTAGttgtg                                                                                                          >  1:318600/1‑146 (MQ=255)
                                            gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGt                                                                                                         >  2:5114140/1‑150 (MQ=255)
                                            gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCCCCCCCCCAAAAAACC‑CCAGGGCCCTCGGGAGCCGAGGGCTGt                                                                                                         >  1:3714218/1‑150 (MQ=255)
                                             cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAAAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGtt                                                                                                        >  2:5652358/1‑150 (MQ=255)
                                              aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAAGCG‑CCCTGGCAGTCGGCAGCGGAGTGCTGTTc                                                                                                       >  1:6825209/1‑150 (MQ=255)
                                              aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTc                                                                                                       >  2:4236727/1‑150 (MQ=255)
                                               tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCg                                                                                                      >  2:5814297/1‑150 (MQ=255)
                                                taatgttttCAAGCTCGACCGCTATACCGCTTTGCTGACCGCCGCCGGCAGCGCCATCGGGGGCGCCGCCGCGGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCgg                                                                                                     <  1:5510084/143‑1 (MQ=255)
                                                 aatttgtTCACGCCCGCCCGCGACCCGGCCTTGCGGCCCGCCGCGGGCCGCGCCAGGGGCGGCGCCGCCGCGGGGGGGGATTTCTCCGGGGCGCGGCGCCGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGt                                                                                                    <  1:2190671/146‑1 (MQ=255)
                                                    gtgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGCGCCGCGCCCCCCATAAAATCG‑CCCTGGCAGTCGGGAGCGTAGTGCTGGTAGGGAcc                                                                                                 >  1:498929/1‑150 (MQ=255)
                                                    gtgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCAAAAATCG‑CCAAGGCAGTCGGCAGCGTAGTGGGGTTGGGTAcc                                                                                                 >  1:983311/1‑150 (MQ=255)
                                                      gtTCACGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTCCTCACCGGCGCTGCGCAGCGCCCCCCATACACGCG‑CCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCtt                                                                                               <  2:1464239/150‑1 (MQ=255)
                                                        tCACGCTCGACCGCGTTACCGCCTTGCTGACCGCCGCGGGCAGCGCCATGTGGGGCGCCGCCGCCGTGCTGGCATTCGACGCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGt                                                                                             <  1:4237080/150‑1 (MQ=255)
                                                        tCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCAGCGCCCCCCCCAAAAGCG‑CCAAGGCACTCGGAAGCCTAGTGGTGGTGGGGAccgagt                                                                                             >  1:4381733/1‑146 (MQ=255)
                                                         cacgCTCGACCGCTATCCCGCGTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTCGATTCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                            <  2:2045234/147‑1 (MQ=255)
                                                          aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAATCGTAGTGCTGTTCGGTACCTTGTCg                                                                                           >  2:5948265/1‑150 (MQ=255)
                                                          aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGGGCCGCGCCCCCCCTAAAAGCG‑CCCAGGCAGGCGGAAGCGTAGTGCTGTTGGGTACCTGGTCg                                                                                           >  1:1139982/1‑150 (MQ=255)
                                                              tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGTTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                       <  2:3709672/150‑1 (MQ=255)
                                                              gcgcccGCGATCCGGCGTTGCTGCCCGCCGCGGGCAGCGCCTTGGGGGGCGCCGCCGCGGTGCGGGCATTCGAAGCGGCGCGGCGCCGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                       <  1:6404586/146‑1 (MQ=255)
                                                                gACCGCTATCCCGCTTTGCTGCCCGCGGCGGGCAGCGCCATCTGGGGCGCCGCCGCGGGGGGGGAATTCGACGCGGGGCTGCGCCGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTc                                                                                     <  1:85661/150‑1 (MQ=255)
                                                                gACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTc                                                                                     <  2:5407304/150‑1 (MQ=255)
                                                                         aCCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGACTCGGCGCGGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCa                                                                            <  2:2545310/150‑1 (MQ=255)
                                                                          ccGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGTGCGCCGCGCCCCCCCGGAAATCGCCCATGGC‑GTAGGAAGGGTTGTGCGGGTGGGTACGGTGTGGGTGGTGCTTTAGCCAc                                                                           >  1:4579514/1‑150 (MQ=255)
                                                                           cGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACt                                                                          <  2:663643/150‑1 (MQ=255)
                                                                            gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCGCCCCCCCCAAAAACG‑CCCAGGCAGTCGGGAGCGTACTGGTGGTGGGGACCATGTGGATGTTGCTTTTACCAcag                                                                         >  1:3150584/1‑148 (MQ=255)
                                                                             ccTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCCTGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTgg                                                                        >  2:6121592/1‑150 (MQ=255)
                                                                              cgtGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCGGGCATTCGACTCGGCGGTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGc                                                                       <  2:4647703/148‑1 (MQ=255)
                                                                                tGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCTTTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTGGTCGATGTTCCTTTACCCACTGGCGa                                                                     <  2:4460377/150‑1 (MQ=255)
                                                                                 gCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCCGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGTCGAt                                                                    >  2:6703587/1‑150 (MQ=255)
                                                                                  cTGACCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCGGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATc                                                                   <  2:1511820/150‑1 (MQ=255)
                                                                                   tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCCTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGGAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCa                                                                  >  2:6864225/1‑150 (MQ=255)
                                                                                   tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGGTCGGTACCATGTCGATGTTCCTTTAACCACTGGCGATCa                                                                  >  2:1289363/1‑150 (MQ=255)
                                                                                       cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGc                                                              >  2:7336254/1‑150 (MQ=255)
                                                                                       cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAttc                                                              >  2:1221463/1‑148 (MQ=255)
                                                                                            ccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct                                                         <  2:550890/150‑1 (MQ=255)
                                                                                                ccgcgcCATCGGCGGCGCCGCCGCCGTGCTGGCATTCGACTCGGCGCGGCGCCGCGCCCCCCATAACAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctggct                                                     <  2:3113739/148‑1 (MQ=255)
                                                                                                 aGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCAAGGCATTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTACTTTACCCACAGGCGATCAATGCCGgctggctg                                                    >  2:1203185/1‑150 (MQ=255)
                                                                                                       cccTGCGGCGCCGCCGCCGTGCTGGCATTCGACTGGGCGGTGCGCCGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTg                                              <  1:778863/148‑1 (MQ=255)
                                                                                                       cacggcggcggcgCCGCAGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGCTGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTg                                              >  2:5809631/5‑150 (MQ=255)
                                                                                                       aTCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCAAAAAAGCG‑CCATGGCCGTCGGAAGCGTAGTGCTGTTGGGGTCCATGTCGATGTTACTTTAACCACTGGCGATAAAATCCGGCTGGCTGCAACTg                                              >  1:1087008/1‑150 (MQ=255)
                                                                                                             ggcgCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCCCCCCCCATCAAAGCG‑CCCAGGCAGTCGGGCGCGGAGAGCAGTTGGGGGCCTTGTCGATGTTGCTATAACCACAGGCGCACAAAACCCGCAGGCTTCATCAGGATaaa                                        >  1:7294957/1‑148 (MQ=255)
                                                                                                                 ccgccgcggTGCTGGAATTCGAATGGGGGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                                    <  1:4917008/150‑1 (MQ=255)
                                                                                                                   gccgccgTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACAGGCGATCAATGCCGGCTGGCTGCAACTGGATACCATGGgc                                  >  1:6077171/1‑150 (MQ=255)
                                                                                                                     cgccgTGCTGGCAGTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCGGTTCGGTACCTTGTCGATGTTCCGTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGgcgc                                >  2:3074416/1‑150 (MQ=255)
                                                                                                                       cgcagatGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATTGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCACAGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCg                              >  2:5259294/7‑150 (MQ=255)
                                                                                                                       ccgTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCg                              <  2:2492150/150‑1 (MQ=255)
                                                                                                                        ggTGCTGGTATTCGAGGCGGGGCTGCGCCGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCgg                             <  1:689843/149‑1 (MQ=255)
                                                                                                                                cATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgc                     <  2:5059128/150‑1 (MQ=255)
                                                                                                                                 aTTCGACTCGGCGCGGCGCAGCGCCCCCCCGAAAAGCG‑CCATGGCAGTCGGAAGCGTAGGGCGGGTCGGTACCGTGTCGATGTTCCGTTACCCACAGGCGATCAAAGCCGGCGGGCGGCATCTGGATACCATGGGCGCCGGgcggctgct                    >  2:937915/1‑150 (MQ=255)
                                                                                                                                 aTTCGACCCGGCGCTGCGCCGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgct                    <  1:613360/150‑1 (MQ=255)
                                                                                                                                  ttCGAATCGGCGCTGCGCAGCGCCCCCCATAACAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg                   <  2:456482/150‑1 (MQ=255)
                                                                                                                                   actaaTCGGCGCTGCGCACCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg                  <  1:3924510/147‑1 (MQ=255)
                                                                                                                                     gatgcGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc                <  1:6511783/146‑1 (MQ=255)
                                                                                                                                     gAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc                <  2:1447237/150‑1 (MQ=255)
                                                                                                                                     gAATCGGCGCGGCGCAGCGCCCCCCAAAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCATGTCGATGTTACATTAACCACATGCGATCAAATCCGGCTGGCAGCATAATGATACCATGGGCGCCCGGCTGCTGCTGggc                >  2:710990/1‑150 (MQ=255)
                                                                                                                                      aagcGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg               <  2:3914610/147‑1 (MQ=255)
                                                                                                                                      aaTCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg               <  2:2318050/150‑1 (MQ=255)
                                                                                                                                         cggcgcggcgCCGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa            <  2:3677890/150‑1 (MQ=255)
                                                                                                                                             gcTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt        >  2:3404605/1‑150 (MQ=255)
                                                                                                                                             gcTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt        >  1:6559318/1‑150 (MQ=255)
                                                                                                                                              cggcgcAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTc       <  1:1773981/148‑1 (MQ=255)
                                                                                                                                                gcgccgcgccCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCa     <  1:6965756/150‑1 (MQ=255)
                                                                                                                                                gcgcAGCGCCCCCCATCAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCCCTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCa     <  2:3503227/150‑1 (MQ=255)
                                                                                                                                                gcgcAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCa     >  1:5915832/1‑150 (MQ=255)
                                                                                                                                                   cAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGa  >  1:2065852/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGA  >  NC_002947/4163316‑4163612

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