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

Only 100 of 105 total aligned reads displayed.

GTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACG  >  NC_002947/4160290‑4160584
                                                                                                                                                   |                                                                                                                                                   
gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                                   <  1:4403063/150‑1 (MQ=255)
gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                                   <  1:4439333/150‑1 (MQ=255)
    tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCAGCTGCGGCGCCGCCGCCGTGCTGGCAGTCGAATCGGCGCGGCGCCGCGccc                                                                                                                                               >  2:3566399/1‑150 (MQ=255)
     ggTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccc                                                                                                                                              >  2:4274519/1‑150 (MQ=255)
     ggTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAGCGCAGCGcccc                                                                                                                                              >  2:4324800/1‑150 (MQ=255)
       tAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCCCGcccccc                                                                                                                                            >  1:4228347/1‑150 (MQ=255)
        aGTCACCAGTACCCTGCTGATCGGCGTGTGGTGTGGCATTAATGTGTTCACGCTCGACCGCGATACCGCCTTGCTGACCGCCGCGGGCTGGGCCATCTGCGGCGCCGCCGCGGTGCTGGCATTCTATTCGGGGATGCGCAGCGccgattc                                                                                                                                           <  2:4303078/150‑6 (MQ=255)
           cccctgTACCCTGCTGATCGGCCTGTGGTTTGGAATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGGTGGCTTTCGAATCGGGGCTGCGCAGCGCCCCCCATaa                                                                                                                                        <  2:614064/145‑1 (MQ=255)
            aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCAGCGcccccccgaaa                                                                                                                                       >  1:2828800/1‑145 (MQ=255)
            aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCACcccccccccaaaa                                                                                                                                       >  1:223617/1‑145 (MQ=255)
             ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATaaaa                                                                                                                                      >  2:823068/1‑150 (MQ=255)
             ccAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCGCCCCCCATaaaa                                                                                                                                      <  2:3268460/150‑1 (MQ=255)
              cAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAg                                                                                                                                     >  1:1301443/1‑150 (MQ=255)
                   cccTGGTGATGGGCCTGTGTTGTGGATTTAATGTTTTAAAGCTCGACGGCTATACCGCTTTGCTGACCGCCGCGGGCTGCGCATTCTGGGGCGCCGCGGCGGTGCTGTCTTTGGAATGGGGGCTGGGCTGCGCCCCCCATAAAAGCGCCa                                                                                                                                <  2:2473349/150‑1 (MQ=255)
                   cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCa                                                                                                                                <  2:1242986/150‑1 (MQ=255)
                    ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTGGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAt                                                                                                                               <  2:3901555/150‑1 (MQ=255)
                    ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCCAAAAAGCGCCaa                                                                                                                               >  1:3157659/1‑149 (MQ=255)
                    ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGccca                                                                                                                               >  1:2175138/1‑148 (MQ=255)
                    ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAACGccca                                                                                                                               >  1:3638095/1‑148 (MQ=255)
                    ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCAAAAAAGCGCCaa                                                                                                                               >  1:3838400/1‑149 (MQ=255)
                    ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCCCCGCGCCCCCCCAAAAAGCGCCag                                                                                                                               >  1:3106858/1‑149 (MQ=255)
                      tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTGCGCAGCGCCCCCCCTAAAAGCGCCCTgg                                                                                                                             >  2:4831259/1‑150 (MQ=255)
                      tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCGCCATgg                                                                                                                             >  2:5077066/1‑150 (MQ=255)
                             cGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCgg                                                                                                                      <  2:1900577/150‑1 (MQ=255)
                                 cTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCAAAAATCGCCCAGGCAGTCGGAAGc                                                                                                                  >  1:4159258/1‑150 (MQ=255)
                                  tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCAAAAAACGCCAGGGCCGTCGGAAGCg                                                                                                                 >  1:777312/1‑150 (MQ=255)
                                    gggtgtggCATGAACGGGTTCACGGGCGCCCGCGATCCCGCCTTGCTGCCCGCCGCCGGCAGCGCCAGCTGCGGCGCCGCCGCCGTGGTGGCATTCGAAGCGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa                                                                                                               <  2:5320814/149‑1 (MQ=255)
                                         tggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCCCCATGGCAGTCGGAAGCGGAGGGCt                                                                                                          >  2:2891813/1‑150 (MQ=255)
                                         tggCATGAAAGGGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCt                                                                                                          <  2:2764047/150‑1 (MQ=255)
                                            cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGACTCGGCGCGGCGCCGCGCCCCCCCTAAAAGCGCCAAGGCAGTCCGTAGCGGAGTGCTGtt                                                                                                       >  2:3896234/1‑150 (MQ=255)
                                             aTGAAAGTGTACAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTc                                                                                                      <  2:5359608/150‑1 (MQ=255)
                                               gAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGGGCAGCGCCCCCCAAAAAACCGCCCAGGCCCTCGGTAGCGTAGGTCTGTTCgg                                                                                                    >  1:2832814/1‑150 (MQ=255)
                                                 aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTAGGTa                                                                                                  >  2:1558206/1‑150 (MQ=255)
                                                 aaGTGTTAAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGAATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTa                                                                                                  <  2:1439654/150‑1 (MQ=255)
                                                    tgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCGAAAAGCGCCATGGCAGGCGGACGCGTAGTGCTGTTCGGTACCt                                                                                               >  2:1733022/1‑150 (MQ=255)
                                                      ttCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTg                                                                                             >  2:2321143/1‑150 (MQ=255)
                                                      ttCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGGACCTTg                                                                                             >  2:750238/1‑150 (MQ=255)
                                                        cacgCTCGACCGCGCGACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGGTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATAACAGCGCCCTGGCAGTCGGAAGCGTAGTGCGGTTCGGTACCTTGTc                                                                                           <  2:2956341/147‑1 (MQ=255)
                                                        cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCCGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTc                                                                                           >  2:4296637/1‑150 (MQ=255)
                                                        cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                           >  1:4870170/1‑150 (MQ=255)
                                                        cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCGAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCGTGTc                                                                                           >  2:4296764/1‑150 (MQ=255)
                                                         aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCCTAAAACCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCCGTACCGTGTCg                                                                                          >  1:5336889/1‑150 (MQ=255)
                                                          cgcgcGACCGCGATACCGCCTGGCGGACCGCCGCCGCCAGCGCCATCTGCGGCGCCGCCGCCGTGCGGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGa                                                                                         <  2:4547104/146‑1 (MQ=255)
                                                            cTCGACCGCGCTCCCGCCTCGCTCCCCGCCGCCGCCCGCGCCCGCGGGGGCGCCCCCGCCGTGGTTGCATTCTACTCGGCGCTGCGCACCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATg                                                                                       <  1:2427230/150‑1 (MQ=255)
                                                            cTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATg                                                                                       <  2:1037548/150‑1 (MQ=255)
                                                               gcccGCGATACCGCCTTGCGGACCGCCGCGGGCAGCGCCATCTGCGGCGCCGCGGCGGTGCTGGATTTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTc                                                                                    <  1:4831259/148‑1 (MQ=255)
                                                               gACCGCGATACCGCCTGGCGGACCGCCGCCGGCAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGCATGCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTc                                                                                    <  2:2010899/150‑1 (MQ=255)
                                                                cccGCTATACCGCGTTGCGGACCGCCGCCGGCCGCGCCCTCGGCGGCGCCGCGGCGGTGGGTGAATTCTACTGGGCGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTcc                                                                                   <  1:3247788/149‑1 (MQ=255)
                                                                    gatctccGCCTTGCTGACCGCCGCCGGCAGCGCCTTCTGGGGCGCCGCCGCGGTGGTGGAATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTa                                                                               <  2:5290101/145‑1 (MQ=255)
                                                                     gaacccGCCTGGCTGCCCGCCGCCGCCAGCGCCATCGGCGGCGCCGCCGCCGTGCGGGAATGCGACGCGGCGCGGCGCAGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAc                                                                              <  1:2789625/146‑1 (MQ=255)
                                                                      ctcccGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAcc                                                                             <  2:4708212/147‑1 (MQ=255)
                                                                         ccGCCTTGCGCCCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGGGGTGTCATTGCAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCAc                                                                          <  1:3274576/150‑1 (MQ=255)
                                                                           gCCTTGCTGACCGCCGCCGGCAGCGCCATGTGCGGCGCCGCCGCCGTGCTGGCATTCTAATCGGCGCTGCTCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTg                                                                        <  2:374031/150‑1 (MQ=255)
                                                                           gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCGCCCCCCCGAAAAGCGCCCTGGCAGTCGGAAGCGTAGTGCTGGTCGGTACCGTGTCGATGTTCCTGTACCCACTg                                                                        >  2:625637/1‑150 (MQ=255)
                                                                                  tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCACAAAAGCGCCCAGGCAGTCGGTAGCGGAGTGGTGTTGGGTACCGTGTCGATGTTCCTTTAACAACAGGCGATCa                                                                 >  1:4793636/1‑150 (MQ=255)
                                                                                  tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGGGCCGCGCCCCCCCCAAAAGCGCCCAGGCAGGCGGAAGCGGAGTGCTGTTAGGTACCGTGTCGCTGTTCCGTTAACCACAGGCGATCa                                                                 >  1:4793869/1‑150 (MQ=255)
                                                                                   gACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCCATGGCGCTCaa                                                                >  2:4406349/1‑150 (MQ=255)
                                                                                   ccccggcgccgGCCGCGCCTTGGGCGGCGCCGCGGCGTTGGTTGCATTCTATTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCaa                                                                <  1:1489628/148‑1 (MQ=255)
                                                                                      cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCCTTCGAATCGGCGCGGCGCAGCGCCCCCCCTCAAAGCGCCATGGCCGTCGGAAGCGTAGTGCTGGTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAGGc                                                             >  2:4264189/1‑150 (MQ=255)
                                                                                      cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGGGCCGCGCCCCCCCGAAAAGCGCCCAGGCCGGCGGAAGCGGGGTGGTGTTCGGTACCGTGTCGGTGTTGCGTTACCCACACGCGCTCAATGc                                                             >  1:3339828/1‑150 (MQ=255)
                                                                                      ccccgccgGCCGCGCCATCGGCGGCGCCGCCGCGGTGGTGGATGTCGACGCGGCGCGGCGCAGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGc                                                             <  1:3896234/148‑1 (MQ=255)
                                                                                       gccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGGAGTGGGGTTCGGTACCTTGTCGTTGTTCCTTTACCCACTGGCGATCAATGcc                                                            >  2:1053879/1‑150 (MQ=255)
                                                                                           ccgGCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGAATTCTAATCGGGGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct                                                        <  2:183917/150‑1 (MQ=255)
                                                                                                     cATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                                              >  2:5388340/1‑150 (MQ=255)
                                                                                                     cATCTGCGGCGCCGCCGCCGTGCGGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAACAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCt                                              <  2:2683681/150‑1 (MQ=255)
                                                                                                      aggggcggcgCCCCCGCGGTGGGGGCATTCGCCGCGGGGCGGCGCCGCGCCCCCCATAAACGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTg                                             <  1:743757/146‑1 (MQ=255)
                                                                                                        cTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCTAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCGGTTCGGTACCGTGTCGCTGTTCCGTTACCCACAGGCGATCAATGCCGGCTGGCTGCATCTGGa                                           >  2:2943409/1‑150 (MQ=255)
                                                                                                        cTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGa                                           <  2:3807988/150‑1 (MQ=255)
                                                                                                         tGCGGCGCCGCCGCCGTGCTGGCTTTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGAt                                          <  2:5058462/150‑1 (MQ=255)
                                                                                                         ggcggcgCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAACCGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGAt                                          <  2:5058438/149‑1 (MQ=255)
                                                                                                          gcggcgCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATa                                         <  1:5344317/150‑1 (MQ=255)
                                                                                                          gcggcgCCGCCGCCGGGGGGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATa                                         <  2:1656697/150‑1 (MQ=255)
                                                                                                           cggcgCCGCCGCGGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGTAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATAc                                        <  2:5171880/150‑1 (MQ=255)
                                                                                                             gcgcCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCa                                      <  2:5439398/150‑1 (MQ=255)
                                                                                                              cgccgccgcCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAt                                     <  2:4787450/150‑1 (MQ=255)
                                                                                                                ccgccgccgTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                                   <  2:3170387/150‑1 (MQ=255)
                                                                                                                ccgccgccgTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATgg                                   <  2:4843611/150‑1 (MQ=255)
                                                                                                                 cgccgccgTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATggg                                  <  2:2474155/150‑1 (MQ=255)
                                                                                                                       cgTGCTGGCATTCGACGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCgg                            <  2:3701594/150‑1 (MQ=255)
                                                                                                                        gTGCTGGCTTGCGACGCGGGGCGGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCggg                           <  1:4964321/150‑1 (MQ=255)
                                                                                                                          gCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgct                         <  2:218965/150‑1 (MQ=255)
                                                                                                                           cTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCGTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctg                        >  2:731154/1‑150 (MQ=255)
                                                                                                                            tGGCATTCGAATCGGCGCTGGGGAGCGCCCCCCATAAAATCGGCATGGCAGTCGGAAGCGTAGTGGTGTTGGGTTCCTTGTGGATGTTGCTTTTGCCAGTTGGGATCAATGCCGGGTGGCTTGATATGGATTCCATGGGCGGCGGggtgg                       >  2:5097551/1‑149 (MQ=255)
                                                                                                                                 ttCGAATCGGCGCTGCGCAGCGCCCCCCATACAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg                  <  2:4637749/150‑1 (MQ=255)
                                                                                                                                  tcgacgcGGCGCTGCGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg                 <  1:3566399/144‑1 (MQ=255)
                                                                                                                                  tCGAATCGGCGCTGCGCAGCGCCCCCCCTAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAATCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTgg                 >  2:101328/1‑150 (MQ=255)
                                                                                                                                     aaTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg              <  2:2175138/150‑1 (MQ=255)
                                                                                                                                      aTCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATTCCGGCTGGCTGCATCTGGATACCAAGGGCGCCGGGCTGCTGCTGggcgg             >  1:508974/1‑150 (MQ=255)
                                                                                                                                        cggcgcggcgCCGCGCCCCCCCTAACAGCGCCATGGCAGGCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa           <  2:4228347/150‑1 (MQ=255)
                                                                                                                                        cggcgcggcgCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa           <  2:34271/150‑1 (MQ=255)
                                                                                                                                           cgcTGCGCCGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTg        <  1:2288515/150‑1 (MQ=255)
                                                                                                                                            gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCg                             >  2:5381727/1‑128 (MQ=255)
                                                                                                                                            gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCg                             <  1:5381727/128‑1 (MQ=255)
                                                                                                                                            gcTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTGGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGt       >  2:3288285/1‑150 (MQ=255)
                                                                                                                                             cTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTc      >  1:1386973/1‑150 (MQ=255)
                                                                                                                                              tGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTcc     >  1:2681627/1‑150 (MQ=255)
                                                                                                                                                cgcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCAc   >  2:1270351/1‑150 (MQ=255)
                                                                                                                                                cgcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCAc   >  2:1270246/1‑150 (MQ=255)
                                                                                                                                                cgcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCAc   >  1:5236481/1‑150 (MQ=255)
                                                                                                                                                 gcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACg  >  1:4955076/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCCACG  >  NC_002947/4160290‑4160584

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