Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,163,474 A→C 11.1% H164P (CAT→CCT)  PP_3661 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,163,4740AC11.1% 484.6 / 15.4 189H164P (CAT→CCT) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base A (78/86);  new base C (23/2);  total (101/88)
Fisher's exact test for biased strand distribution p-value = 1.88e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

AGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGGCGCAGG  >  NC_002947/4163325‑4163623
                                                                                                                                                     |                                                                                                                                                     
aGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTCGAATCGGCGCTGCGCAGCGCCCCCCa                                                                                                                                                       <  2:4143925/150‑1 (MQ=255)
  tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGGGGCGCGGCCCCCaga                                                                                                                                                     >  1:2073079/1‑148 (MQ=255)
  tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCaaa                                                                                                                                                     >  2:348193/1‑148 (MQ=255)
     cctgTACCCTGCTGATCGGCCTGTGGTGTGGCTTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGGTGACCGCCGCCGGCTGCGCATTCTGCGGCGCCGCGGCCTTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATaaaa                                                                                                                                                  <  2:70394/147‑1 (MQ=255)
     ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATaaaa                                                                                                                                                  <  2:2094784/150‑1 (MQ=255)
      cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGTCGCAGCGCCCCCCagcaaac                                                                                                                                                 >  1:4482228/1‑144 (MQ=255)
       aGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCAAAAAAgc                                                                                                                                                >  1:3710544/1‑150 (MQ=255)
         tACCCTGCTGATCGGCCTGGGGTGTGGCATGAACGGGTTCACGCTGGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCGGGCTTTGGAATCGGGGCTGCGCAGCGCCCCCCATAAAAgcgc                                                                                                                                              <  2:4056283/150‑1 (MQ=255)
            ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAt                                                                                                                                           <  2:4614799/150‑1 (MQ=255)
            ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCAAAAAAGCGCCAt                                                                                                                                           >  2:2514484/1‑150 (MQ=255)
            ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCTAAAAGCGCCAt                                                                                                                                           >  2:5337593/1‑150 (MQ=255)
            ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCCCCCCCCCAAAAAACCCCaa                                                                                                                                           >  1:201077/1‑149 (MQ=255)
              tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCAAgg                                                                                                                                         >  1:4600386/1‑150 (MQ=255)
              tgctgATCGGCCTGTGGTGTGGCATGAAAGTGATCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCCTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCAtcg                                                                                                                                         >  2:5558509/1‑148 (MQ=255)
              tgctgATAGGCCTGTCGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCATGCTGACCGCCACAGGCAGAGCCATCTGCGCCGCCGACGCATAGCTGCCAAACGAATCGTCCCAGCGCAGAGCCCCCCAAAAAAGCGCCaatg                                                                                                                                         >  2:1874367/1‑147 (MQ=255)
                cctaTCTGCCTGTGGTGTGGCATGAATGTGTTAAAGCTCGACCGCGTTACGGCGTTGCTGACCGCCGCGGGCTGCGCAATGTGGGGCGCGGCCGCGGTGGTGGCATTCTATTCGGGGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCa                                                                                                                                       <  2:1136256/147‑1 (MQ=255)
                 tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCACCGCCCCCCCCAAAAACCCCAACGCAg                                                                                                                                      >  1:2458629/1‑150 (MQ=255)
                  gATCGGCCTGTGGTGTGGCTTGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCTTATGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGt                                                                                                                                     <  2:4205296/150‑1 (MQ=255)
                    tCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCAGCGCCCCCCCCAAAAGCGCCCAGGCAGGCg                                                                                                                                   >  1:5196808/1‑150 (MQ=255)
                       gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGCCCCGCGCCCCCCAAAAAAAACCCCAGGCAGTCGGaa                                                                                                                                >  1:3189551/1‑150 (MQ=255)
                            tggtgtggCATGAATGTTTTCAAGCTCGACCGCGATACCGCCTTGCTGCCCGCCGCCGGCCGCGCCATGTGCGGCGCCGCCGCCGTGCTGGCATTCGATTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                                           <  2:150088/150‑1 (MQ=255)
                            tggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCgga                                                                                                                           >  2:3545183/1‑148 (MQ=255)
                             ggtgtggCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAg                                                                                                                          <  2:632627/150‑1 (MQ=255)
                              gtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGt                                                                                                                         >  2:4633069/1‑150 (MQ=255)
                                 tggCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCCTTGTTGGCATTCGAATCGGCGCTGCGAAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCt                                                                                                                      <  2:439563/150‑1 (MQ=255)
                                  ggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCCCCCCCCAAAAAACCGCCAAGGCAGTCGGAAGCCTAGTGCTg                                                                                                                     >  1:3361671/1‑150 (MQ=255)
                                     aTGAACGTGTTCAAGCTCGACCGCGATACCGCCTGGCTGCCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTc                                                                                                                  <  2:2850422/150‑1 (MQ=255)
                                     aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCGGTTc                                                                                                                  >  2:2377334/1‑150 (MQ=255)
                                         aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCACCGCCCCCCAAAAAAACGCCAAGGCACTCGGAACCGGAGTGCTGTTAGGTa                                                                                                              >  1:4486714/1‑150 (MQ=255)
                                           gtgtTCAAGCTCGACCGCGATACCGCCTTCCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCACCCCTAAAAGCGCCAAGGCACTCGGAAGCGAAGTTCTGTGAGGTAcc                                                                                                            >  2:2940242/1‑150 (MQ=255)
                                             gtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCCCCCCCCAAAAAAGCCCCCAGGCACTCGGAACCGGAGTGGTGGTCGGTAccgt                                                                                                          >  1:2886428/1‑148 (MQ=255)
                                               tCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGGAGTGCTGTTCGGTACCATGt                                                                                                        >  2:3206678/1‑150 (MQ=255)
                                                cAAGCTCGACCGCTATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                                       <  2:1821776/150‑1 (MQ=255)
                                                cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCCCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                                       <  2:2412519/150‑1 (MQ=255)
                                                cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCAAAAAGCGCCATGGCAGTCGGAAGCGGAGTGGTGGTCGGTACCgtgtg                                                                                                       >  1:5014195/1‑149 (MQ=255)
                                                 aaGCTCGACCGCTTTACCGCCTTGCTGCCCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCACCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCg                                                                                                      <  1:1258546/150‑1 (MQ=255)
                                                 aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCGCCCCCCCCAAAAGCGCCCCGGCACGGGGGAGCGTAGTGCGGTTGGGGACCGTGGCg                                                                                                      >  1:702479/1‑150 (MQ=255)
                                                 aaGCACGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTAGGCAGCGTAGTGCGGTTCGGTACCTTGTCg                                                                                                      >  2:2395696/1‑150 (MQ=255)
                                                   gCTCGCCCGCTTTCCCGCCTTGCTGCCCGCGGCCGGCAGCGCCCTGGGGGGCGCCGCCGCCGGGGTGTCTTTCGACTCGGGGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                                    <  1:781766/150‑1 (MQ=255)
                                                     tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTGGTTGt                                                                                                  >  1:454232/1‑150 (MQ=255)
                                                     tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCGAAAATCGCCCGGGCCGTCGGAAGCGTAGTGGTGTTCGGTACCGTGTCGgtgt                                                                                                  >  1:2800317/1‑150 (MQ=255)
                                                      cGACCGCGATACCGCCTTGCTGACCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGtt                                                                                                 <  2:3133225/150‑1 (MQ=255)
                                                         ccGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCCCCGCGCCCCCCAGAAAAGCGCCCTGGCAGTACGAAGCGGAGTGCGGTTAGGTACCTGGGCGATGTTCCt                                                                                              >  2:5495443/1‑150 (MQ=255)
                                                          cgcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCAGTTCGGTACCTTGGCGATGTTccat                                                                                             >  2:5176060/1‑148 (MQ=255)
                                                            cgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGCTTTCGACGGGGCGCTGCCCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTa                                                                                           <  1:1287356/150‑1 (MQ=255)
                                                                cccGCTTTGCGGCCCGCCGCCGGCAGCGCCATGGGGGGCGCCGCCGCGGTGCTGGCATTCGATGCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCa                                                                                       <  1:315020/149‑1 (MQ=255)
                                                                  cGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGGGGCGCAGCGCCCCCCAGAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGGTGGTGGGTACCTTGTCGATGTTGCGTTAACCAca                                                                                     >  1:360662/1‑149 (MQ=255)
                                                                  cGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCACCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACt                                                                                     >  1:30279/1‑150 (MQ=255)
                                                                   gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTg                                                                                    <  2:2318515/150‑1 (MQ=255)
                                                                    ccTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAGAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGGTGTTCGGTACCATGTCGAGGTTCCGTTAACCACAgg                                                                                   >  1:1365106/1‑150 (MQ=255)
                                                                       tggtGACGGCGGCCGGCTGCGCGATGTGCGGCGCGGCCGCGGTGCGGGATTTCGATTGGGGGCTGCGCGGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGa                                                                                <  1:5495443/147‑1 (MQ=255)
                                                                        gCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGAt                                                                               <  2:3139849/150‑1 (MQ=255)
                                                                            cccgccgccgGCAGCGCCCTCGGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAt                                                                           <  1:4205643/149‑1 (MQ=255)
                                                                              ctccgccgGCTGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGc                                                                         <  2:22053/148‑1 (MQ=255)
                                                                                ccgccgGCATCGCCATGTGGGGGGCCGCCGCCGTGCTGGCATTCTAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCagcaggg                                                                       <  2:2008046/150‑7 (MQ=255)
                                                                                ccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGGCGTGCTGTTCGGTACCTTGTCGAGGGTCCTTTACCCACTGGCGCTCAATGCCg                                                                       >  2:4786625/1‑150 (MQ=255)
                                                                                ccgccgGCAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCg                                                                       <  2:3861487/150‑1 (MQ=255)
                                                                                     ggCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAATCGCCCTGGGAGTCGGAATCGTTGTGGTGTTGGGTACCTTGTGGATGTTGCTTTTACCACTGGCGATCAATTCCggctgg                                                                  >  2:2487831/1‑150 (MQ=255)
                                                                                     ggCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctgg                                                                  >  2:4231492/1‑150 (MQ=255)
                                                                                      gCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctggc                                                                 <  2:2267841/150‑1 (MQ=255)
                                                                                       cAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGGTGTTCGGTACCTTGTCGATGTTACTTTACCCACAGGCGATAAAATCCggctggct                                                                >  1:2783085/1‑150 (MQ=255)
                                                                                       cAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctggct                                                                <  2:4496017/150‑1 (MQ=255)
                                                                                             cttCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                                                          <  2:1158475/148‑1 (MQ=255)
                                                                                                gggcggcgCCGCCGCGGTGCGGGCATTCTTCGCGGCGCGGCGCCGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                                       <  1:1895245/148‑1 (MQ=255)
                                                                                                cTGCGGCGCCGCCGCCGTGGTGGATTTCGAGTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                                       <  1:3677117/150‑1 (MQ=255)
                                                                                                cTGCGGCGCCGCCGCCGTGCTGGAATTCTAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                                       <  2:1400573/150‑1 (MQ=255)
                                                                                                  gcggcgCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATa                                                     >  2:1767117/1‑150 (MQ=255)
                                                                                                    ggcgCCGCCGCCGTGCTGGCATTGGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATAcc                                                   <  2:1616227/150‑1 (MQ=255)
                                                                                                    ggcgCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCCTGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCGTTACCCACGGGCGATCAAAGCCGGCTGGCTGCATCTGGATAcc                                                   >  2:4967241/1‑150 (MQ=255)
                                                                                                       gccgccgccGTGCTGGCATTCGAACCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATg                                                <  2:1917676/150‑1 (MQ=255)
                                                                                                        ccgccgccgTGCTGGCATTCGACTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                                               <  1:4625504/150‑1 (MQ=255)
                                                                                                               cgTGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCgg                                        <  2:2601258/150‑1 (MQ=255)
                                                                                                                  gCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATACAAGCGCCATGGCAGTCGGACGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACGGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgct                                     <  2:4119311/150‑1 (MQ=255)
                                                                                                                  gCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgct                                     <  2:1484620/150‑1 (MQ=255)
                                                                                                                    tGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgc                                   >  2:4591021/1‑150 (MQ=255)
                                                                                                                    tGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACAGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGggcggc                                   >  2:455988/1‑150 (MQ=255)
                                                                                                                         ttcgccgcggcgcggcgCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg                              <  1:348193/143‑1 (MQ=255)
                                                                                                                         ttCTAAGCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg                              <  1:5337593/150‑1 (MQ=255)
                                                                                                                         tacgccgcGGCGCTGCGCCCCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg                              <  1:2514484/143‑1 (MQ=255)
                                                                                                                          tCGACTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg                             <  1:5627443/150‑1 (MQ=255)
                                                                                                                          tCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAACCGTAGTGCTGTTCGGTACCATGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCggg                             >  2:3218374/1‑150 (MQ=255)
                                                                                                                            gattgggcgcggcgcTGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc                           <  2:702102/145‑1 (MQ=255)
                                                                                                                            gAATCGGCGCTGCGCAGCGCCCCCCAAAAAATCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATTTTCAATTAACCACAGGCGATCAATGCCCGCTGGCTGCATCTGGATAACAATGGCGCCGGGCTGCTGCTGggc                           >  1:827577/1‑150 (MQ=255)
                                                                                                                            gAATCGGCGCGGCGCCGCGCCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc                           <  2:4516961/150‑1 (MQ=255)
                                                                                                                             aagcGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg                          <  2:224205/147‑1 (MQ=255)
                                                                                                                             aaTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTAGGTACCTTGTCGATGTTCCTTTACCCACAGGCGATCAATTCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg                          >  1:382451/1‑150 (MQ=255)
                                                                                                                              aTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcgg                         >  2:869591/1‑150 (MQ=255)
                                                                                                                                 ggCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCAc                      >  1:1633747/1‑150 (MQ=255)
                                                                                                                                  gcgcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt                     >  2:3618466/1‑150 (MQ=255)
                                                                                                                                  gcgcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt                     >  2:1509978/1‑150 (MQ=255)
                                                                                                                                    gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt                   >  1:3869242/1‑150 (MQ=255)
                                                                                                                                     cTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTc                  >  2:271138/1‑150 (MQ=255)
                                                                                                                                         gcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACg              >  1:3673475/1‑150 (MQ=255)
                                                                                                                                           cccgcCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGaa            <  1:4252242/148‑1 (MQ=255)
                                                                                                                                              gcCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTg         <  2:5399519/150‑1 (MQ=255)
                                                                                                                                               ccccccATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTgg        >  1:5416014/1‑150 (MQ=255)
                                                                                                                                                 ccccATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGgcg      >  2:1660906/1‑150 (MQ=255)
                                                                                                                                                  cccATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGgcgc     >  1:1186428/1‑150 (MQ=255)
                                                                                                                                                   ccATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGGCGCa    >  2:4028500/1‑150 (MQ=255)
                                                                                                                                                     aTAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGGCGCAgg  >  1:1414876/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAAGTGGCGCAGG  >  NC_002947/4163325‑4163623

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