Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,160,437 A→C 20.6% 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,160,4370AC20.6% 307.8 / 25.8 190S161R (AGC→CGC) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base A (70/71);  new base C (20/22);  total (92/98)
Fisher's exact test for biased strand distribution p-value = 8.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 193 total aligned reads displayed.

CGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAA  >  NC_002947/4160288‑4160586
                                                                                                                                                     |                                                                                                                                                      
  gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                                      <  2:2087414/150‑1 (MQ=255)
  gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCGGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                                      <  2:1833472/150‑1 (MQ=255)
          aGTCCCCAGTACCCTGCTGTGCGGCCTGTGGTGTGGCATGAAAGTTTTCATGCTCGCCCGCTATACCGCGTGGCTGACGGCCGCGGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCTATTCGGCGCTgcgccgcg‑ccgattc                                                                                                                                             <  2:1270538/150‑6 (MQ=255)
           gTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CCCCCCAt                                                                                                                                            >  2:4796574/1‑150 (MQ=255)
            tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑cccccccac                                                                                                                                           >  1:7787564/1‑147 (MQ=255)
             cACCAGTACCCTGCTGATGGGCCTGTGGTGTGGCATTAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCTTTCGATTCGGCGCTGCGCTGCG‑CCCCCCATaa                                                                                                                                          <  2:661550/150‑1 (MQ=255)
              aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATaaa                                                                                                                                         <  2:2775216/150‑1 (MQ=255)
              aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAaaaa                                                                                                                                         >  1:8515202/1‑150 (MQ=255)
               ccAGTACCGTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGTTACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTCGAATCGGCGCTGCGCTGCG‑CCCCCCATaaaa                                                                                                                                        <  2:8918253/150‑1 (MQ=255)
               ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCaaaaaa                                                                                                                                        >  2:8557876/1‑150 (MQ=255)
               ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATaaaa                                                                                                                                        >  2:26091/1‑150 (MQ=255)
               ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATaaaa                                                                                                                                        >  2:25836/1‑150 (MQ=255)
                cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGTGGGCAGCG‑CGCCCCATAAAAg                                                                                                                                       >  1:3714161/1‑150 (MQ=255)
                cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCCCAAAAg                                                                                                                                       >  1:4458515/1‑150 (MQ=255)
                 aGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAgc                                                                                                                                      <  2:5836456/150‑1 (MQ=255)
                     cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAAGCGCCa                                                                                                                                  >  2:3934612/1‑150 (MQ=255)
                      ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CCCCCCCTAAAAGCGCCAt                                                                                                                                 >  2:9056078/1‑150 (MQ=255)
                      ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCAAAAAGCGCCAt                                                                                                                                 >  2:2484626/1‑150 (MQ=255)
                      ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCCAAAAACCGccca                                                                                                                                 >  1:4650955/1‑148 (MQ=255)
                      ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAAAAAAGCGCCaa                                                                                                                                 >  2:7948709/1‑149 (MQ=255)
                        tgctgTTGGGCTTGTGGTGTGGCATGAAAGTGTTCATGCTCGACGGCGATACCGCGTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTGGAATCGGCGCTGCGCTGCG‑CCCCCCATAAAAGCGCCATgg                                                                                                                               <  2:8104995/150‑1 (MQ=255)
                        tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCC‑CCCCCCCAAAAATCGCCCTgg                                                                                                                               >  1:624952/1‑150 (MQ=255)
                         gctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCTAAAAGCGcccaggc                                                                                                                              >  1:9107547/1‑145 (MQ=255)
                         gctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCCAAAAAGCGcccaggc                                                                                                                              >  1:2298556/1‑145 (MQ=255)
                                 gcgtgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                                      <  2:6170827/147‑1 (MQ=255)
                                 gCCTGTGGTGTGGCTTGAAAGTGTTCAAGCTCGCCCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                                      <  2:9015805/150‑1 (MQ=255)
                                 gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                                      <  2:7082765/150‑1 (MQ=255)
                                 gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                                      >  2:3688902/1‑150 (MQ=255)
                                  ccTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CCCCCCAAAAAAAACCCAAGGCAGTCGGAAg                                                                                                                     >  1:4060853/1‑150 (MQ=255)
                                  ccTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAAGCGCCAAGGCAGACGGaaa                                                                                                                     >  1:7616753/1‑149 (MQ=255)
                                    tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCg                                                                                                                   <  2:2121039/150‑1 (MQ=255)
                                    tgtggtgtggATTGAATGTTTTCAAGCCCGACGGCGATACGGCGTTGCTGACCCCGGCGGGCAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGCTTTCTATGCGGGGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCg                                                                                                                   <  1:2796172/150‑1 (MQ=255)
                                    tgtggtggggCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATGTGCGGCGCCGCCGCCGTGCTGGCTTTCGACTGGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCg                                                                                                                   <  1:2795606/150‑1 (MQ=255)
                                      tggtgtggCTTGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCCGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                                 <  1:8388276/150‑1 (MQ=255)
                                         tgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCTTGTGGGGCGCCGCCGCCGTGGTGGCATTCTAATCGGCGCTGCGCAGCG‑CTCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTg                                                                                                              <  1:4086036/150‑1 (MQ=255)
                                         tgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CTCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTg                                                                                                              <  1:4085946/150‑1 (MQ=255)
                                           gggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATCAAAGCGCCATGGCAGTCGGAAGCGTAGTGCt                                                                                                            <  2:3171868/149‑1 (MQ=255)
                                             gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGGATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCCTGGCAGTCGGAAGCGTAGTGCTGt                                                                                                          >  2:1850498/1‑150 (MQ=255)
                                              cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCCCGCCCCCCCAAAAAACCCCCAGGGC‑GTCGGAAGCGGAGTGCTGtt                                                                                                         >  1:3121188/1‑150 (MQ=255)
                                               aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTTGTGCTGTTc                                                                                                        >  2:2836051/1‑150 (MQ=255)
                                               aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTc                                                                                                        >  2:7846012/1‑150 (MQ=255)
                                               aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCCCAAAAACGCCATGGCAGTCGGAAGCCGAGTGCGGTTc                                                                                                        >  1:7352018/1‑150 (MQ=255)
                                                ttattgtgtTCAAGCTCGACCGCGATACCGCCGTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAACCGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCg                                                                                                       <  1:1495443/145‑1 (MQ=255)
                                                tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCG‑CCCCCCCAAAAAGCGCCCTGGCCGTCGGAAGCGTAGTGGTGTTCg                                                                                                       >  2:3343990/1‑150 (MQ=255)
                                                 gataGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCTTCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCgg                                                                                                      <  1:6864294/147‑1 (MQ=255)
                                                 gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGGAGTGCTGTTCgg                                                                                                      >  2:5376424/1‑150 (MQ=255)
                                                    aGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAc                                                                                                   <  2:5640295/150‑1 (MQ=255)
                                                    aGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAc                                                                                                   <  2:5640638/150‑1 (MQ=255)
                                                     gtgtTCACGCTCGCCCGCGCGCCCGCCTTGCTGCCCGCCGCCGGCCGCGCCATGGGCGGCGCCGCCGCGGTGCGGGCTTTCGACGCGGGGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAcc                                                                                                  <  1:7969214/150‑1 (MQ=255)
                                                        ttCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAGAAAAGCGCCAGGGCAGTCGGGAGCGTAGTGCTGTTCGGTAccaag                                                                                               >  1:8201450/1‑147 (MQ=255)
                                                         tCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGt                                                                                              <  2:1742608/150‑1 (MQ=255)
                                                          cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                             >  2:3811238/1‑150 (MQ=255)
                                                          cAAGCTCCACCGCCATACCGCCTTGCTGACCGCCGCCGGCAGCGCAATCTGCGGCGCCGCCGCCGTGCTGCAATTCAAATCGGCGCTCCGCACCC‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                             <  1:319156/150‑1 (MQ=255)
                                                            aGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCACCG‑CCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGGAGTGGAGTTCGGTACCGTGTcca                                                                                           >  1:8309155/1‑148 (MQ=255)
                                                             gCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                          >  2:7923211/1‑150 (MQ=255)
                                                               tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTGGTCGATGt                                                                                        >  2:2509218/1‑150 (MQ=255)
                                                               ccgcccGCGATCCCGCGTTGCGGCCCGCGGCGGGCTGCGCCTTGGGGGGCGCCGCCGCGGGGGGGGATTTCGACGCGGGGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                        <  1:4620497/146‑1 (MQ=255)
                                                                  cccGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTcc                                                                                     <  1:5964590/149‑1 (MQ=255)
                                                                   ccGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCt                                                                                    <  2:2882462/150‑1 (MQ=255)
                                                                    cgcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGATTCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCtt                                                                                   <  2:7044384/150‑1 (MQ=255)
                                                                       gATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCG‑CCCCCCCAAAAATCGCCAAGGCCGGCGGGAGAGTAGTTGTGTTGGGGACAATGTCGATGTTGCGTTaa                                                                                >  1:4846818/1‑149 (MQ=255)
                                                                        aTACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCG‑CCCCCCAAAAAAGCGCCCAGGCCGTCGGAAGCGGAGTGCTGTTGGGGACCATGTCGATGTTGCGTTaaa                                                                               >  1:3432730/1‑148 (MQ=255)
                                                                           ccGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGACTCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGGAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCAc                                                                            <  2:5762557/150‑1 (MQ=255)
                                                                           ccGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCAc                                                                            <  2:2818361/150‑1 (MQ=255)
                                                                            cGCCTTGCTGACCGCCGCCGGCAGCGCCTTCTGCGGCGCCGCCGCCGTGCGGGCATTCGAAGCGGCGCGGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACt                                                                           <  2:3271148/150‑1 (MQ=255)
                                                                            cGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACt                                                                           >  2:9167805/1‑150 (MQ=255)
                                                                             gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTg                                                                          <  2:8414258/150‑1 (MQ=255)
                                                                                   cTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCG‑CCCCCCATAAAATCGCCAAGGCAGGCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCATTACCCACAGGCGATc                                                                    >  1:7470467/1‑150 (MQ=255)
                                                                                      aCCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAt                                                                 <  2:1562287/150‑1 (MQ=255)
                                                                                      aCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCTGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAt                                                                 <  2:1272777/150‑1 (MQ=255)
                                                                                        gcacgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCCAAAAAGCGCCCAGGCCGCCGGCAGCGGAGTGCTGTTCGGTACCTTGTCGAGGTTGCTGTTCCCACAGGCGAACacagc                                                               >  1:6356429/4‑146 (MQ=255)
                                                                                           cgccgGCCGCGCCATCTGCGGCGCCGCCGCCTTGGTGGCATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCgg                                                            <  2:13515/150‑1 (MQ=255)
                                                                                               ggCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctgg                                                        <  2:7738364/150‑1 (MQ=255)
                                                                                                       cATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                                                <  2:6736071/150‑1 (MQ=255)
                                                                                                       cATCGGGGGCGCCGCCGCCGTGCTGGCATTCGTCGCGGGGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                                                <  1:1309121/150‑1 (MQ=255)
                                                                                                          cggcggcgCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                             <  2:7193036/148‑1 (MQ=255)
                                                                                                                 gccgccgccGTTCTGGAATTCGAATGGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATg                                      <  1:1896685/150‑1 (MQ=255)
                                                                                                                  ccgccgcggTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                                     <  2:7017015/150‑1 (MQ=255)
                                                                                                                         cctGCTGGAATTCGAATCGGCGCTCCCCACCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCgg                              <  2:8955342/148‑1 (MQ=255)
                                                                                                                          gTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCggg                             <  2:591252/150‑1 (MQ=255)
                                                                                                                          gTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCggg                             <  1:6811277/150‑1 (MQ=255)
                                                                                                                           ggCTGGCATCCGACTCGGCGCTGCGCAGCG‑CCCCCCATACCAGCGCCATGGCCGTCGGAAGCGTCGTGCGGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgc                            <  2:3668191/149‑1 (MQ=255)
                                                                                                                                gCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAAGCCCGCTGGCTGCATCTGGAAACCAAGGGCGCCGGgctgcag                       >  2:6835926/1‑148 (MQ=255)
                                                                                                                                   ttCGAATCGGCGCTGCGCAGCG‑CCCCCCAAAAAATCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAACCACAGGCGATAAAATCCGGCTGGCTGCATCTGGATAACAAGGGCGCCGGGctgctgctg                    >  2:4227513/1‑150 (MQ=255)
                                                                                                                                    tatggggggCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTACTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg                   >  1:6121921/8‑150 (MQ=255)
                                                                                                                                    tCGATTGGGGGCGGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg                   <  1:1337946/150‑1 (MQ=255)
                                                                                                                                    tCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg                   <  2:4132738/150‑1 (MQ=255)
                                                                                                                                        gggggCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAAGGGCGCCGGGCTGCTGCTGggcgg               >  1:1001851/4‑150 (MQ=255)
                                                                                                                                        cgcggcgcggcgCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcgg               <  1:4978001/148‑1 (MQ=255)
                                                                                                                                        agcGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcgg               <  2:624982/148‑1 (MQ=255)
                                                                                                                                        aTCGGCGCTGCGCAGCG‑CCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCAACTGGATACCAAGGGCGCCGGGCTGCTGCTGggcgg               >  1:8629689/1‑150 (MQ=255)
                                                                                                                                          gggCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa             >  1:435756/2‑150 (MQ=255)
                                                                                                                                          gggCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAAGGGCGCCGGGCTGCTGCTGGGCGGCa             >  1:435461/2‑150 (MQ=255)
                                                                                                                                              gcTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt         >  2:7273526/1‑150 (MQ=255)
                                                                                                                                              gcTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt         >  2:4437646/1‑150 (MQ=255)
                                                                                                                                               cTGCGCCGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTc        <  1:4922752/150‑1 (MQ=255)
                                                                                                                                                tGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTcc       >  1:3752437/1‑150 (MQ=255)
                                                                                                                                                ggcgcAGCG‑CCCCCCATAAAAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTcc       <  2:7372600/149‑1 (MQ=255)
                                                                                                                                                 gcgcAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCa      >  2:905981/1‑150 (MQ=255)
                                                                                                                                                     aGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGaa  >  1:7160726/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCG‑CCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGAA  >  NC_002947/4160288‑4160586

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