Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,160,437 A→C 16.0% 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,4370AC16.0% 243.1 / 10.7 138S161R (AGC→CGC) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base A (54/55);  new base C (15/10);  total (69/69)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

GTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGA  >  NC_002947/4160290‑4160585
                                                                                                                                                   |                                                                                                                                                     
 tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCcgcg                                                                                                                                                    >  2:2907320/1‑150 (MQ=255)
    tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGccc                                                                                                                                                 >  2:1795169/1‑150 (MQ=255)
     ggTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccc                                                                                                                                                <  2:326190/150‑1 (MQ=255)
      gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccccc                                                                                                                                               >  2:4132391/1‑150 (MQ=255)
      gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGccccc                                                                                                                                               >  2:2991811/1‑150 (MQ=255)
         gTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGccccccct                                                                                                                                            >  2:826301/1‑148 (MQ=255)
          tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGGGCAGCGCCCCCCATa                                                                                                                                           >  2:832753/1‑150 (MQ=255)
              cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCAAAAAg                                                                                                                                       >  1:531476/1‑150 (MQ=255)
                 tACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCTCAAAgcgc                                                                                                                                    >  2:4355077/1‑150 (MQ=255)
                 tACCCTGCTGATCGGCCTGGGGTGTGGCATGAATGTGTTCAAGCTCGACCGCGTTACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCTGCGCCCCCCATAAAAgcgc                                                                                                                                    <  2:1302897/150‑1 (MQ=255)
                   cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTTACCGCCGCCGGCAGCGCCATCTGCGGCGCGGCCGCGGTGCTGGATTTCGAATGGGCGGTGCGCAGCGCCCCCCATAAAAGCGCCa                                                                                                                                  <  2:717185/150‑1 (MQ=255)
                   cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATGTGCGGCGCCGCCGCCGTGCTGGCATTCGATTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCa                                                                                                                                  <  2:3837753/150‑1 (MQ=255)
                    ccTGCTGTTCGGCCTGTGGTGTGGCATGAATGTGTTAAAGCTGGACCGCGTTACCGCCTTGATGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGGTGGCTTTCGAATGGGCGATGCGCTGCGCCCCCCATAAAAGCGCCAt                                                                                                                                 <  2:2362740/150‑1 (MQ=255)
                    ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGGGCCGCGCCCCCCCGAAAATCGcccg                                                                                                                                 >  1:1117648/1‑148 (MQ=255)
                    ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCCCCCCCCATAAAAGCGCCaa                                                                                                                                 >  2:77767/1‑149 (MQ=255)
                      tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCAGCGCCCCCCATAAAAGCGCCAAgg                                                                                                                               >  2:2550927/1‑150 (MQ=255)
                      tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCCTgg                                                                                                                               >  2:1851081/1‑150 (MQ=255)
                      tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCACAAAGCGCCATgg                                                                                                                               >  2:236035/1‑150 (MQ=255)
                      tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATgg                                                                                                                               >  2:5089817/1‑150 (MQ=255)
                      tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAAgg                                                                                                                               >  2:677897/1‑150 (MQ=255)
                      tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCATgg                                                                                                                               >  2:2520246/1‑150 (MQ=255)
                      tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCAAAAAACGCCAAgg                                                                                                                               >  1:2337843/1‑150 (MQ=255)
                            tCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCCAAAAACCGCCAAGGCAGtag                                                                                                                         >  1:2675952/1‑148 (MQ=255)
                                ccTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTGGGCATCGGCCCCCATAAAATAGGCATGGGAGTGGGTAg                                                                                                                     >  1:4638282/1‑150 (MQ=255)
                                 cTGTGGTGTGGCATGTATGTGTTCAAGCTCGCCCGCGATACCGCGTTGCTGACCGCCGCGGGCTGCGCCATCGGCGGGGCCGCCGCCGTGCTGGCATTCGAATGGGCGCTGCGAAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGc                                                                                                                    <  2:201209/150‑1 (MQ=255)
                                 cTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGc                                                                                                                    >  2:3907833/1‑150 (MQ=255)
                                    tggtgtggCATGAATGTTTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCAATCTGCGGCGCCGCCGCGGTGCTGTCATTCTATTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                                 <  2:1300608/150‑1 (MQ=255)
                                    tggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTGGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                                 <  2:3503027/150‑1 (MQ=255)
                                         tggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCt                                                                                                            <  2:5012146/150‑1 (MQ=255)
                                         gggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCTTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCt                                                                                                            <  2:198580/149‑1 (MQ=255)
                                           gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGGAGTGCTGt                                                                                                          >  1:426526/1‑150 (MQ=255)
                                           gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAATCCCCATGGCAGTCGGAAGCGTAGTGCGGt                                                                                                          >  1:4580079/1‑150 (MQ=255)
                                           gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCGAAAAGCCCCATGGCAGTCGGAAGCGGAGTGgtgt                                                                                                          >  1:1585359/1‑150 (MQ=255)
                                           gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTCCTGGCATTCGAATCGGCGCTGCCCAGCGCCCCCCATAAAACCGCCAAGGCAGTCGGAAGCCGAGGGCTGt                                                                                                          >  2:359182/1‑150 (MQ=255)
                                            cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGtt                                                                                                         >  2:2039858/1‑150 (MQ=255)
                                             attaaaGTTTTCAAGCTCGACCGCGATACCGCGTTGCTGCCCGCCGCCGGCAGCGCCTTCTGGGGCGCCGCCGCGGTGCTGGAATTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTc                                                                                                        <  1:810793/147‑1 (MQ=255)
                                             aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCAAAAAACGCCCGGGCAGTCGGAAGCGGAGGTCTGTTc                                                                                                        >  1:1498000/1‑150 (MQ=255)
                                                    tgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCt                                                                                                 <  2:1127620/150‑1 (MQ=255)
                                                    tgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCt                                                                                                 <  2:1532369/150‑1 (MQ=255)
                                                        cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCCTTCGAATCGGCGCTGCGCAGCGCCCCCCCTAAAAGCGCCATGGCCGTCGGAAGCGTAGTGCTGTTCGGTACCATGTc                                                                                             >  2:1534808/1‑150 (MQ=255)
                                                          cgCTCGACCGCTATCCCGCGTTGCTGCCCGCGGCCGGCTGCGCCATGTGGGGGGCCGCCGGGGTGTTTGCTTTCTACTGGGGGCGGGGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGa                                                                                           <  1:305793/149‑1 (MQ=255)
                                                          aGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGa                                                                                           >  2:3834563/1‑150 (MQ=255)
                                                           gCTCGCCCGCTATACCGCCTTGGGGCCCGCCGCCGGCCGCGCCTTGTGCGGCGCCGCCGCGTTGGTGGAATTCTAAGGGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                          <  1:1273778/150‑1 (MQ=255)
                                                           gCTCGACCGCGATCCCGCGTTGCGGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                          <  2:3954112/150‑1 (MQ=255)
                                                           gCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                          <  2:2266305/150‑1 (MQ=255)
                                                           gCTCGACCGCGATACCGCCTGGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCACCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGAt                                                                                          <  2:662000/150‑1 (MQ=255)
                                                            cTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGGGCTGGCATTCGACTCGGCGCGGCGCAGCGCCCCCCCGACACGCGCCATGGCAGTCGGACGCGTAGGGCTGGTCGGTACCTGGTCGATg                                                                                         >  2:4106409/1‑150 (MQ=255)
                                                             tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCCGAAAAACCCCCTGGCCCGCGGGAGCGGTGTGGTGTTGGGGGCCTTGGGGATGt                                                                                        >  1:3911912/1‑150 (MQ=255)
                                                             tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGGAGTGCGGTTCGGTACCTGGTCGCTGt                                                                                        >  2:2403147/1‑150 (MQ=255)
                                                               gcccGCGATCCCGCCTGGCTGCCCGCCGCCGGCTGCGCCCTCGGGGGCGCCGCCGCGGGGCGGGACTTGGCCGGGGCGGGGCGCCGCGCCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTc                                                                                      <  1:3139710/148‑1 (MQ=255)
                                                               gcccGCGATACCGCGTTGCTCCCCGCCGCCGGCAGCGCACTCGGCGGCGCCGCCGCCGGGCTGGCATTCGAAGCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTc                                                                                      <  1:2877800/148‑1 (MQ=255)
                                                                cccGCTAGACCGCGTTGCGGCCCGCCGCCGGCCGCGCCCGGGGGGGCGCCGCGGCGGTGGTGGATTTCGACGGGGGGGGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTcc                                                                                     <  1:2253512/149‑1 (MQ=255)
                                                                aCCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTcc                                                                                     <  2:3852885/150‑1 (MQ=255)
                                                                 ccGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCCCGCCCCCCCAAAAACCGCCAAGGCAGGCGGAAGCGGGCTGCGGTTCGGTACCGTGTAGATGttaca                                                                                    >  1:298011/1‑147 (MQ=255)
                                                                  cgctATCCCGCGTTGCTGCCCGCGGCGGGCTGCGCCATGGGCGGGGCCGCCGCCGTGCTGGCATTCGTATCGGCGCGGCGCCGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCt‑t                                                                                  <  1:557131/150‑1 (MQ=255)
                                                                  cgcgATACCGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATACCAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCt‑t                                                                                  <  2:3408186/150‑1 (MQ=255)
                                                                  cgcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCTATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCt‑t                                                                                  <  2:2926167/150‑1 (MQ=255)
                                                                   gcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCGAAAAGCGCCCAGGCCCGCGGAAGAGGAGTGGGGTTGGGGTCCGTGTCGGTGTTCCg‑tt                                                                                 >  1:2192465/1‑150 (MQ=255)
                                                                    cgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCA‑TTa                                                                                >  2:4051532/1‑150 (MQ=255)
                                                                     gATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAGAAAAGCGCCAAGGCCGTCGGAAGCGTAGTGCTGTTCGGTACCATGTCGATGGTCCg‑tttc                                                                               >  1:768634/1‑148 (MQ=255)
                                                                        aCCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCa                                                                            <  2:2623611/150‑1 (MQ=255)
                                                                            ccTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTgg                                                                        <  2:2768872/150‑1 (MQ=255)
                                                                               tGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGGTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGa                                                                     <  2:4843859/150‑1 (MQ=255)
                                                                                gCTCACCCCCGCCGGCAGCCCCAACTGCGGCGCCGCCGCCGTGCTGCCATCCTAATCGGCGCTGCGCACCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGAt                                                                    <  2:200442/150‑1 (MQ=255)
                                                                                 gtGACCGCCGCCGGCAGCGCCTTGTGGGGCGCCGCCGCCGTGGTGGAATTCGACTGGGGGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATc                                                                   <  1:2413265/149‑1 (MQ=255)
                                                                                  tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATAAAATCGCCATTGCAGTCGGAAGCGTATTGCTTTTGGGTAACATGTAGATGTTCCT‑TTAAACACATGCGAtaa                                                                  >  1:2033555/1‑148 (MQ=255)
                                                                                    aCCGCCGCCGGCTGCGCCTTCTGCGGCGCCGCCGCGGTGCTGGATTTCTATTGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAAt                                                                <  2:640330/150‑1 (MQ=255)
                                                                                       gccgccgGCAGCGCCATCTGCGGCGCCGCCGCGGTGGTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGcc                                                             <  2:3276511/150‑1 (MQ=255)
                                                                                        ccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTGGGAATCGTAGTGCTGTTGGGTTCCTTGTGGATTTTACT‑TTTCACAATGGCGATAAAttacg                                                            >  1:5074141/1‑146 (MQ=255)
                                                                                        ccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCAGCGCCCCCCCGCAAAACGCCAGGGCAGTCGGGAGCGGATTGGTGTTCGGTAACTTGTCGATGTTGCT‑TTAAACACAAGCGATCAAAGCCg                                                            >  1:987412/1‑150 (MQ=255)
                                                                                             ggCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCCGAAAAGCGCCAAGGCAGTGGGAACCGTAGTGCTGGTCGGTCCCGGGTCGATGTTGCTGTTACCC‑CAGGCGGGCAAAGCCggcggg                                                       >  1:3065467/1‑150 (MQ=255)
                                                                                              gccgcgcAACCTGCGGCGCCGCCGCCTTGCGGGCTTGCGAAGCGGCGCTGCGCAGCGCCCCCCATAAATGCTCAATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACGGGCGATAAATGCCgtctgga                                                      <  2:4716164/147‑2 (MQ=255)
                                                                                              gCAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCggctggc                                                      <  2:4716341/150‑1 (MQ=255)
                                                                                               cAGCGCCATGTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCggctggct                                                     <  2:867757/150‑1 (MQ=255)
                                                                                               cAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCCATGTTCCG‑TTACCCACAGGCGATCAAAGCCggctggct                                                     >  2:2602314/1‑150 (MQ=255)
                                                                                               cAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCTCCATGGCAGTCGGAAGCGTAGTGCGGGTCGCTACCTCGTCGATTTTCCT‑TTACCCACGGGCGATCAATGCCggctggct                                                     <  2:867765/150‑1 (MQ=255)
                                                                                                  cgcCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGGTAGGTACCTTGTCGAATTTACT‑TTAACCACATGCGATCAAAGCCGGCAGGCTGCa                                                  >  2:4415634/1‑150 (MQ=255)
                                                                                                     cATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                                               >  2:3698542/1‑150 (MQ=255)
                                                                                                             gcgcCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCa                                       >  2:759330/1‑150 (MQ=255)
                                                                                                                ccgccgccgTGCTGGATTTCGACTGGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                                    <  1:1004825/150‑1 (MQ=255)
                                                                                                                                aTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAATCCCCATGGCATTCGGAAGCGTAGTGCTGTTCGGTACCTTTTCGATGTACCA‑TTACCAAATTTCGATAAATGCCGGCTTGCTGAAACTTTATACCATGGGCGCCGGgctgctgct                    >  2:3399006/1‑150 (MQ=255)
                                                                                                                                 ttCGAATCGGCGCTGCCCACCGCCCCCCAAAAAATCGCCAAGGCAGTAGGAAGCGTAGTGCTGTTAGGTACCATGTCGATTTTACT‑TTAACCACTGGCGATCAATTCCCGCTGGCTGCATCATGATAACATGGGCGCCGGGctgctgctg                   >  1:303737/1‑150 (MQ=255)
                                                                                                                                  ccgacgcGGCGCCGCGCCCCGCCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg                  <  1:4355077/144‑1 (MQ=255)
                                                                                                                                   cGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCT‑TTACCCACAGGCGATCTATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTggg                 >  2:4865837/1‑150 (MQ=255)
                                                                                                                                    gAATCGGCGCTGCGCAGCGCCCCCCAAAAAACCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCATGTCGATTTTCCT‑TTAAACACAGGCGATCAAAGCCGGCTGGCTGCATCAGGATACCATAGGCGCCGGGCGGCTTCTGggc                >  1:3996848/1‑150 (MQ=255)
                                                                                                                                    gAATCGGCGCTGCGCAGCCCCCCCCATAAAAGCGCCATGGCAGGCCGAAGCGTAGTGCTGGTCGGTACCGTGGCGGTGTGCCG‑TTACCCACGGGCGGTCATTGCCGGCTGGCTGCCGCTGGCTACCATGGGCGCCGGGCGGCGGCGGggc                >  1:2109237/1‑150 (MQ=255)
                                                                                                                                        cggcgcggcgAAGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa            <  1:1947082/150‑1 (MQ=255)
                                                                                                                                        cGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa            <  1:233239/150‑1 (MQ=255)
                                                                                                                                           cgcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTg         >  2:1426763/1‑150 (MQ=255)
                                                                                                                                           cgcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTg         <  2:1102297/150‑1 (MQ=255)
                                                                                                                                           cgcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTg         >  2:4005717/1‑150 (MQ=255)
                                                                                                                                            gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCTTGGGCGCCGGGCTGCTGCGGGGCGGCACTGt        >  2:4803626/1‑150 (MQ=255)
                                                                                                                                            gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt        >  2:4803524/1‑150 (MQ=255)
                                                                                                                                            gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt        >  2:2214726/1‑150 (MQ=255)
                                                                                                                                             cTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTc       <  2:225444/150‑1 (MQ=255)
                                                                                                                                              tGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTcc      >  2:3141317/1‑150 (MQ=255)
                                                                                                                                              tGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTcc      >  1:428153/1‑150 (MQ=255)
                                                                                                                                               gcgcAGCGCCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCa     <  2:3257632/150‑1 (MQ=255)
                                                                                                                                               gcgcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGTCTGCTGCTGGGCGGCACTGTCCa     >  2:3795440/1‑150 (MQ=255)
                                                                                                                                               gcgcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCa     <  2:531476/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCT‑TTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACGA  >  NC_002947/4160290‑4160585

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