Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,160,437 A→C 16.4% 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.4% 282.4 / 9.2 147S161R (AGC→CGC) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base A (64/53);  new base C (11/17);  total (75/72)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

CGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCAC  >  NC_002947/4160288‑4160583
                                                                                                                                                     |                                                                                                                                                   
  gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAg‑c                                                                                                                                                  <  2:3869061/150‑1 (MQ=255)
     cTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTgcgccg‑cgcc                                                                                                                                               >  1:1072276/1‑150 (MQ=255)
      tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CAccc                                                                                                                                              <  1:2353127/150‑1 (MQ=255)
        gtagtaACCAGTACCCTGCTGATCGGCGTGTGGTGTGGCATTAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCTTCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGccccc                                                                                                                                            <  1:2608408/150‑1 (MQ=255)
        gTAGTCACCAGTACCCTGTTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGccgat                                                                                                                                            <  2:2644177/150‑4 (MQ=255)
        gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGTATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCG‑CGccccc                                                                                                                                            >  2:607358/1‑150 (MQ=255)
            tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAG‑CGCCCCCCaga                                                                                                                                        >  2:1685375/1‑148 (MQ=255)
            tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAG‑CGCCCCCCaaa                                                                                                                                        >  2:575598/1‑148 (MQ=255)
              aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑cccccccccaaaa                                                                                                                                      >  1:360038/1‑145 (MQ=255)
              aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCAaaaa                                                                                                                                      >  1:4014501/1‑150 (MQ=255)
              aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAG‑CGcccccccaaaa                                                                                                                                      >  1:3213771/1‑145 (MQ=255)
              aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATCCGAATCGGCGCGGCGCCC‑CGcccccccccaa                                                                                                                                      >  1:2683293/1‑145 (MQ=255)
               ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATaaaa                                                                                                                                     <  2:2536380/150‑1 (MQ=255)
               ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCGGCCGCCGTGCTGGCATTGGAATCGGCGCTGCGCAG‑GGCCCCCCATaaaa                                                                                                                                     <  2:2681179/150‑1 (MQ=255)
               ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCaaaaaa                                                                                                                                     >  2:3153425/1‑150 (MQ=255)
                cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCG‑CGCCCCCCATCAAAg                                                                                                                                    >  2:1607799/1‑150 (MQ=255)
                 aGTACCCTGCTTATCGGCCTGTGGTGTGGCATGAATGTGTTCTAGCTCGACCGCGTTACCGCCTTGCTTACCGCCGCCGGCTGCGCCTTCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAgc                                                                                                                                   <  2:2197059/150‑1 (MQ=255)
                        tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGGGCCC‑CGCCCCCCCAAAAAGcccccagg                                                                                                                            >  1:1615413/1‑146 (MQ=255)
                             atgggCGTGTGGTGTGGATTGAATGTGTTCACGCTCGCCCGCGATCCCGCGTTGCTGACCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCGGTGGTGGCTTTCGAATCGGCGCTGCGCCG‑CGCCCCCCATAAAAGCGCCATGGCAGTc                                                                                                                       <  1:4030331/147‑1 (MQ=255)
                              tgggCCTGTGGTGTGGCTTGAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACGGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCg                                                                                                                      <  2:398076/148‑1 (MQ=255)
                                 gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCTATACCGCGTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGAAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                                   <  2:3902853/150‑1 (MQ=255)
                                 gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGGGCAG‑CGCCCCCCAAAAAAACCCCAAGGGAGTCGGaa                                                                                                                   >  1:851719/1‑150 (MQ=255)
                                 gCCTGTGGTGGGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGaa                                                                                                                   <  2:2707822/150‑1 (MQ=255)
                                   cTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGAATCGGCGCTGCGCAG‑CGCCCCCCCGAAAAGCGCCAAGGCAGTCGGGAGc                                                                                                                 >  1:3004498/1‑150 (MQ=255)
                                    tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGCGCAG‑CGCCCCCCCAAAAAGCGCCAAGGCAGACGGAAGCg                                                                                                                >  2:391467/1‑150 (MQ=255)
                                      tggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                              <  2:3238684/150‑1 (MQ=255)
                                        gtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCG‑CGCCCCCCAAAAAAGCGCCCAGGCAGTCCGAAGCGTAGt                                                                                                            >  1:3696301/1‑150 (MQ=255)
                                         tgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTg                                                                                                           >  2:2065043/1‑150 (MQ=255)
                                         tgtggCATGAAAGTGTTAAAGCTCGACCGCGATACCGCCTTGGTGACCGCCGCCGGCAGCGCCATGTGGGGCGCCGCCGCCGTGGTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTg                                                                                                           <  2:4522191/150‑1 (MQ=255)
                                           tggCATTAATGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCTTCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCt                                                                                                         <  2:407877/150‑1 (MQ=255)
                                             gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCCGAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGt                                                                                                       >  2:309814/1‑150 (MQ=255)
                                              cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCAATCTGCGGCGCCGCCGCCGTGCTGGCTTTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGtt                                                                                                      <  2:1228576/150‑1 (MQ=255)
                                               aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAG‑CGCCCCCCCGAAAAGCGCCCTGGCAGTCGGAAGCGGAGTGCGGTTc                                                                                                     >  1:3899502/1‑150 (MQ=255)
                                                tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCCTAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCg                                                                                                    >  2:3739100/1‑150 (MQ=255)
                                                 gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGGGCAGCCCCCCCCCCAAAAAGCGCCAAGGGAGGCGGAAGCCAAgtggtggtgg                                                                                                    >  1:3824060/1‑146 (MQ=255)
                                                 gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCCAAAAAGCGCCAGGGCAGTCCGAACCCTAGTGCTGTTCgg                                                                                                   >  2:2446151/1‑150 (MQ=255)
                                                 gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCgg                                                                                                   >  2:391937/1‑150 (MQ=255)
                                                  aatgtgtACACGCTCGCCCGCGAGCCCGCCTTGCTGACCGCCGCGGGCAGCGCCTGCGGCGGCGCCGCCGCCGTGCTGTCATTCGCAGCGGCGCGGCCCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGt                                                                                                  <  1:2501280/147‑1 (MQ=255)
                                                  aaaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGt                                                                                                  <  2:2389293/150‑1 (MQ=255)
                                                   aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTa                                                                                                 <  2:3219170/150‑1 (MQ=255)
                                                     gtgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTAcc                                                                                               >  2:2968009/1‑150 (MQ=255)
                                                      tgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCt                                                                                              >  2:4502556/1‑150 (MQ=255)
                                                       gtcctaGCTCGACCGCGATACCGCTTTGCTGACCGCCGCCGGCAGCGCATTCTGGGGGGCCGCCGCCGTGGTGGCATTCGATCCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCtt                                                                                             <  1:1392481/145‑1 (MQ=255)
                                                       gtcaaaGCTCGCCGGCGATACCGCCTGGCTGCCCGCCGCCGGCCGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAAGCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCtt                                                                                             <  1:1392657/146‑1 (MQ=255)
                                                       gtTCACGCTCGCCCGCGAGCCCGCCTTGCTGCCCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCCGGGCTGGCATTCGACTCGGCGCGGCGCCG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCtt                                                                                             <  1:3234496/150‑1 (MQ=255)
                                                         taaaGCCCGACCGCGATACGGCTTTGCTGCCCGCCGCCGGCAGCGCATTGTGGGGCGCCGCCGCCGTGCTGTCTTTCTATTGGGGGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGt                                                                                           <  1:3578227/148‑1 (MQ=255)
                                                          cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                          >  2:1585072/1‑150 (MQ=255)
                                                          cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCAAAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                          >  2:483339/1‑150 (MQ=255)
                                                            aGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCG‑CGCCCCCCAAAAAAGCGCCCTGGCAGGCGGAAGCGGAGTGCTGTTCGGTACCGTGTCGa                                                                                        >  2:2916085/1‑150 (MQ=255)
                                                            aGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATTGCAGTCGGAATCGTATTGCTGTTCGGTACCTTTTCGa                                                                                        >  2:2251503/1‑150 (MQ=255)
                                                            aGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTAGa                                                                                        >  2:2251531/1‑150 (MQ=255)
                                                               tCGCCCGCTATACCGCCTTGCTGCCCGCCGCCGGCAGCGCCATCGGCGGCGCCGCCGGCGTGCTGGCATTCGTAGCGGCGCGGCGCCG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                     <  1:2835647/150‑1 (MQ=255)
                                                               tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                                     >  2:3418696/1‑150 (MQ=255)
                                                                   cggcgATACCGCCTTGGTGACCGCCGCGGGCAGCGCCATCGGCGGCGCGGCCGCGGGGGTGGCTGTCGACTGGGCGCTGCCCCG‑CCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCt                                                                                 <  1:1123296/148‑1 (MQ=255)
                                                                    cgcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCCTTCGAATCGGCGCGGCGCCG‑CGCCCCCCAGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCtt                                                                                >  2:3409745/1‑150 (MQ=255)
                                                                    cgcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCtt                                                                                <  2:3741192/150‑1 (MQ=255)
                                                                    cgcgATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTGGTCGATGTTCCtt                                                                                >  2:1821264/1‑150 (MQ=255)
                                                                          aCCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCa                                                                          >  2:2294600/1‑150 (MQ=255)
                                                                          aCCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCAAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTCGATGTTCCTTTACCCa                                                                          >  2:2294759/1‑150 (MQ=255)
                                                                             gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCG‑CGCCCCCCAAAAAAGCGCCCTGGCAGTCGGAAGCGTAGTGCTGGTCGGTACCTTGTAGATGTTCCTTTAACCAcag                                                                       >  2:1248331/1‑148 (MQ=255)
                                                                              ccTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTgg                                                                      <  2:3632773/150‑1 (MQ=255)
                                                                                 tGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGa                                                                   <  2:827588/150‑1 (MQ=255)
                                                                                  gCTGACCGCCGCGGGCAGCGCCATCTGCGGCGCCGCCGCGGTGGTGGAATTCGACTCGGGGCTGCGCCG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGAt                                                                  <  1:1773460/150‑1 (MQ=255)
                                                                                  gCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGAt                                                                  <  2:4092555/150‑1 (MQ=255)
                                                                                   cTGACCGCCGGCGGCTGCGCCATGTGGGGCGCGGCCGCGGTGGTGGATTTCTCTTGGGGGGTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATc                                                                 <  1:4276073/150‑1 (MQ=255)
                                                                                      aCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAt                                                              <  2:2867487/150‑1 (MQ=255)
                                                                                      aCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAG‑CGCCCCCCAAAAAAGCCCCAAGGCAGTCGGACGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCGTTACACACAGGCGATaaaa                                                              >  1:4443776/1‑149 (MQ=255)
                                                                                       ccgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATg                                                             <  2:681679/150‑1 (MQ=255)
                                                                                        cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCCAAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCGGTTGGGTTCCTTGTGGATGGTCCTTTAACCACTGGCGGTCAAAGc                                                            >  1:4505031/1‑150 (MQ=255)
                                                                                        cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGGTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCCATGc                                                            >  2:3482745/1‑150 (MQ=255)
                                                                                         gccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAG‑CGCCCCCCATAAACGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGcc                                                           <  2:3357694/150‑1 (MQ=255)
                                                                                             ccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct                                                       <  2:1444858/150‑1 (MQ=255)
                                                                                                          cggcggcgCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                          <  2:2821528/148‑1 (MQ=255)
                                                                                                          cTGCGGCGCCGCCGCCGTGCTGTCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                          <  2:851719/150‑1 (MQ=255)
                                                                                                                 gccgccgcgGTGGTGGATTTCGACGCGGGGCTGCGCCG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATg                                   <  1:2053837/150‑1 (MQ=255)
                                                                                                                 gccgccgccGTGCTGGCATTCGAATCGGCGCTGCGCCG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATg                                   <  2:392871/150‑1 (MQ=255)
                                                                                                                   cgccgccgTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATAAATTCCGGCTGGCTGCATCTGGATACCATggg                                 >  2:732056/1‑150 (MQ=255)
                                                                                                                    gccgccgTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCAAAAAATCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTAGATTTTCCATTACCCACAGGCGATCAAATCCGGCTGGCTGCATCTGGATAACAAGGgc                                >  2:1399273/1‑150 (MQ=255)
                                                                                                                    gccgccgGGGGGGAATTCTAAGCGGCGCTGCGCCG‑CCCCCCCCATAAAAGCGCCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGgc                                <  1:603196/150‑1 (MQ=255)
                                                                                                                         ggTGCTGGCATTCGAGCCGGCGGGGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCgg                           <  1:4177330/149‑1 (MQ=255)
                                                                                                                          gTGCTGGCATTCGATTCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCggg                          <  2:4187125/150‑1 (MQ=255)
                                                                                                                          gTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCggg                          <  2:3675863/150‑1 (MQ=255)
                                                                                                                          gTGCTGGATTTCGAAGCGGCGCGGCGCCG‑CGCCCCCCATAAAAGCGCCCTGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCggg                          <  1:3739100/150‑1 (MQ=255)
                                                                                                                              tGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCAACTGGATACCATGGGCGCCGGgctgc                      >  2:3964362/1‑150 (MQ=255)
                                                                                                                              tGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCGTTACCCACAGGCGATCAATGCCGGCTGGCAGCATCTGGATACCAAGGGCGCCGGgctgc                      >  2:2492826/1‑150 (MQ=255)
                                                                                                                              tGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCAATGCAGTCGGAAGCGTAGTGCTGTTCGGTACCATGTCGAAGTTCCTTTACCCACAGGCGATAAAATCCCGCTGGCTGCATATGGATAACATGGGCGCCGggcggc                      >  2:2903684/1‑150 (MQ=255)
                                                                                                                                gCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCAGAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCATGTCGATGTTCCTTTAACCACAGGCGATCAAAGCCGGCTGGCGGCAAATGGATACCATGGGCGCCCGgctgctg                    >  1:3721080/1‑150 (MQ=255)
                                                                                                                                 cATTCGAAGCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgc                   <  2:2606621/150‑1 (MQ=255)
                                                                                                                                  aTTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCATTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgct                  >  2:2134634/1‑150 (MQ=255)
                                                                                                                                   ttCGACTCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg                 <  2:4013392/150‑1 (MQ=255)
                                                                                                                                      gattCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc              <  1:1106471/147‑1 (MQ=255)
                                                                                                                                       caccGGCGCTGCCCCC‑CCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg             <  1:392063/147‑1 (MQ=255)
                                                                                                                                        aggggcgcggcgcCG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcgg            <  1:607358/147‑1 (MQ=255)
                                                                                                                                            gcgcTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt        >  2:546155/1‑150 (MQ=255)
                                                                                                                                            gcgcTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt        >  2:1003728/1‑150 (MQ=255)
                                                                                                                                            gcgcTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACt        >  1:1974655/1‑150 (MQ=255)
                                                                                                                                              gcTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt      >  2:2985681/1‑150 (MQ=255)
                                                                                                                                                tGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTcc    >  2:640558/1‑150 (MQ=255)
                                                                                                                                                ggcgcAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTcc    <  2:1143123/149‑1 (MQ=255)
                                                                                                                                                 gcgcAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCa   <  2:2703510/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAG‑CGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCAC  >  NC_002947/4160288‑4160583

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