Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,160,432 T→G 22.3% L159R (CTG→CGG)  PP_3661 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,160,4320TG22.3% 147.0 / 9.4 101L159R (CTG→CGG) PP_3661membrane protein
Reads supporting (aligned to +/- strand):  ref base T (33/42);  new base G (15/10);  total (49/52)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

GTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCC  >  NC_002947/4160286‑4160581
                                                                                                                                                  |                                                                                                                                                     
gtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTgcg                                                                                                                                                    <  1:4150014/150‑1 (MQ=255)
 tCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTgcgc                                                                                                                                                   >  1:4547104/1‑150 (MQ=255)
   ggTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTgcgccg                                                                                                                                                 >  1:4843611/1‑148 (MQ=255)
   ggTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCCGCAGCGCCATCTGCGGCGCCGCCGCCGAGCTGGCATACGAAACCCCGCTgcgccc                                                                                                                                                 >  2:3310025/1‑148 (MQ=255)
    gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                                <  1:4439333/150‑1 (MQ=255)
    gTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgc                                                                                                                                                <  1:4403063/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)
                                                                                       gACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGAATCGGCGCGGCGCCGCGCCCCCCCGACAAGCGCCCAGGCAGGCGGACGCGGCGTGCGGGGCGGTACCGGGGCGAGGGTCCGGTACCCACGGGCGCTCac                                                             >  2:3047128/1‑149 (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)
                                                                                               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)
                                                                                                           tCTGCGGCGCCGCCGCCGTGCTGGCATTCGACGCGGCGCGGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTgg                                         <  2:1199538/150‑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)
                                                                                                                    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)
                                                                                                                                                  |                                                                                                                                                     
GTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCACTGTCC  >  NC_002947/4160286‑4160581

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