Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,348,608 C→G 5.4% L20L (CTC→CTG) ‡ PP_4703 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,348,6080CG 5.4% 496.9 / 13.3 149L20L (CTC→CTG) ‡PP_4703hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (70/71);  new base G (6/2);  total (76/73)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 151 total aligned reads displayed.

CGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCCCCTGGCCCTG  >  NC_002947/5348459‑5348755
                                                                                                                                                     |                                                                                                                                                   
cGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTc                                                                                                                                                     <  1:1781191/150‑1 (MQ=255)
 gACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCg                                                                                                                                                    <  1:3701317/150‑1 (MQ=255)
 gACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCg                                                                                                                                                    <  2:301266/150‑1 (MQ=255)
 gACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCg                                                                                                                                                    <  1:2820476/150‑1 (MQ=255)
 gACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCg                                                                                                                                                    <  1:2167226/150‑1 (MQ=255)
  aCCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCgg                                                                                                                                                   >  2:3788398/1‑150 (MQ=255)
    ccGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGcc                                                                                                                                                 >  2:3457267/1‑150 (MQ=255)
      gAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAg                                                                                                                                               <  1:1756296/150‑1 (MQ=255)
      gAGCGTGGTACAGGGCCTGACCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAg                                                                                                                                               <  1:5124767/150‑1 (MQ=255)
         cGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTc                                                                                                                                            <  1:522879/150‑1 (MQ=255)
              tACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCt                                                                                                                                       <  2:5143912/150‑1 (MQ=255)
              tACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCt                                                                                                                                       <  2:4810188/150‑1 (MQ=255)
                cAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTgg                                                                                                                                     <  2:2306641/150‑1 (MQ=255)
                cAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTgg                                                                                                                                     <  1:2771273/150‑1 (MQ=255)
                cAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTgg                                                                                                                                     <  2:4653468/150‑1 (MQ=255)
                  gggCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGaa                                                                                                                                   <  2:4700417/150‑1 (MQ=255)
                   ggCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAAt                                                                                                                                  >  1:1844863/1‑150 (MQ=255)
                   ggCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAAt                                                                                                                                  >  1:5434415/1‑150 (MQ=255)
                         aTCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCa                                                                                                                            >  2:3730443/1‑150 (MQ=255)
                            gACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCGAGGGCCAGCCCg                                                                                                                                           <  1:3293064/132‑1 (MQ=255)
                            gACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAgg                                                                                                                         >  1:5382632/1‑150 (MQ=255)
                             aCAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCCGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAggg                                                                                                                        >  2:3146802/1‑150 (MQ=255)
                              cAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCTTGGGc                                                                                                                       <  1:264340/150‑1 (MQ=255)
                              cAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGc                                                                                                                       <  1:953479/150‑1 (MQ=255)
                               aaaCGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCg                                                                                                                      >  2:1528104/1‑150 (MQ=255)
                               aaaCGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCg                                                                                                                      >  2:2424758/1‑150 (MQ=255)
                               aaaCGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCg                                                                                                                      >  1:1787519/1‑150 (MQ=255)
                               aaaCGCCCCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTg                                                                                                                                      >  1:2865632/1‑134 (MQ=255)
                               aaaCGCCCCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTg                                                                                                                                      <  2:2865632/134‑1 (MQ=255)
                                aaCGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGa                                                                                                                     <  1:1744942/150‑1 (MQ=255)
                                   gCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGAccc                                                                                                                  >  1:1863131/1‑150 (MQ=255)
                                     cTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGcc                                                                                                                                                 >  2:2780019/1‑117 (MQ=255)
                                     cTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCTTGGGCGattcca                                                                                                                <  1:2283077/150‑6 (MQ=255)
                                      tCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCTTGGGCGattccag                                                                                                               <  1:1336532/150‑7 (MQ=255)
                                      tCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCgg                                                                                                               >  2:802847/1‑150 (MQ=255)
                                       cAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGc                                                                                                              <  1:2580988/150‑1 (MQ=255)
                                           cAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGt                                                                                                          >  2:1468861/1‑150 (MQ=255)
                                            aGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGtt                                                                                                         >  2:1261249/1‑150 (MQ=255)
                                                 agggcggcCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGcc                                                                                                    >  1:820206/3‑150 (MQ=255)
                                                  cggcggCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGATACCCATGTGTTGAGCCTCTCCGGGCTGGCGCTCGGCCAGCTCGACCTGGACTCGCCCAAGGGCGACCCCGGCCTGTTCGGccc                                                                                                   <  2:2796735/150‑1 (MQ=255)
                                                  cggcggCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGccc                                                                                                   <  1:4812571/150‑1 (MQ=255)
                                                         cGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCTTGGGCGat                                                                                                                    <  1:928596/126‑2 (MQ=255)
                                                         cGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCTTGGGCGat                                                                                                                    >  2:928596/1‑125 (MQ=255)
                                                         cGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCt                                                                                                            <  1:5307155/134‑1 (MQ=255)
                                                         cGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGt                                                                                            >  2:4366182/1‑150 (MQ=255)
                                                          gCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTa                                                                                           >  1:2576585/1‑150 (MQ=255)
                                                          gCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTa                                                                                           >  1:3801300/1‑150 (MQ=255)
                                                            aaaTCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTAtt                                                                                         >  2:1494727/1‑150 (MQ=255)
                                                               tCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGc                                                                                      >  1:1545681/1‑150 (MQ=255)
                                                                  ggCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTgg                                                                                   <  2:1385466/150‑1 (MQ=255)
                                                                       cccTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGAGCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGt                                                                              <  2:588466/150‑1 (MQ=255)
                                                                        ccTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTc                                                                             >  1:3951455/1‑150 (MQ=255)
                                                                         cTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTCTCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGACTCGCCCAAGGGCGACCCCGGCCTGGTCGGCCCGCACAGGCTTAGCTGGCGCGTcc                                                                            >  2:3004246/1‑150 (MQ=255)
                                                                             aCAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGAcccg                                                                                                                                      <  1:2938091/88‑3 (MQ=255)
                                                                             aCAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTg                                                                                                                                      >  2:2938091/1‑88 (MQ=255)
                                                                             aCAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATgg                                                                        <  1:1481092/150‑1 (MQ=255)
                                                                             aCAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATgg                                                                        >  2:1481092/1‑150 (MQ=255)
                                                                              cagcagGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCg                                                                                                                                <  1:2488928/93‑1 (MQ=255)
                                                                              cagcagGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCg                                                                                                                                >  2:2488928/1‑93 (MQ=255)
                                                                                gcagGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGa                                                                     <  2:4070579/150‑1 (MQ=255)
                                                                                        ccATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTCGCTGGCGCGTCCATGGCGACTTCCCGa                                                             >  2:3540579/1‑150 (MQ=255)
                                                                                         cATGGACGCCCTTCGTCGCCGGATCGAAACCCCGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGCCCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGGCCATGGCGACGTCCCgcg                                                            >  2:5363685/1‑148 (MQ=255)
                                                                                          aTGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGc                                                           >  2:795286/1‑150 (MQ=255)
                                                                                             gAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATg                                                        <  1:1067172/150‑1 (MQ=255)
                                                                                             gAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATg                                                        >  1:3007397/1‑150 (MQ=255)
                                                                                             gAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACGTCCCGAGCCTg                                                        >  2:2869211/1‑150 (MQ=255)
                                                                                                gCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCATCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTg                                                     >  2:1158315/1‑150 (MQ=255)
                                                                                                  ccTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCt                                                                                                                                       <  2:2146417/66‑1 (MQ=255)
                                                                                                  ccTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTACCGACCagcatgct                                                   <  2:1863131/150‑8 (MQ=255)
                                                                                                   cTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTgggc                                                  >  2:4962330/1‑148 (MQ=255)
                                                                                                      cgtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGt                                               <  1:5440573/150‑1 (MQ=255)
                                                                                                      cgtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGt                                               <  1:4968091/150‑1 (MQ=255)
                                                                                                       gtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTg                                              <  1:1026951/150‑1 (MQ=255)
                                                                                                       gtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTg                                              <  1:2525462/150‑1 (MQ=255)
                                                                                                       gtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTg                                              <  1:3767531/150‑1 (MQ=255)
                                                                                                       gtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTg                                              >  2:5152537/1‑150 (MQ=255)
                                                                                                       gtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTg                                              >  2:5003469/1‑150 (MQ=255)
                                                                                                       gtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGGCCATGGCGACTTCCCGAGCATGCTGGGCgggg                                              >  2:5003400/1‑148 (MQ=255)
                                                                                                                      cccAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGa                               <  1:795286/150‑1 (MQ=255)
                                                                                                                         aGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGAtgc                            <  2:3951455/150‑1 (MQ=255)
                                                                                                                          ggTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGAtgct                           <  2:2015271/150‑1 (MQ=255)
                                                                                                                           gTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGGTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGAtgctg                          >  2:1843513/1‑150 (MQ=255)
                                                                                                                            tCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGAtgctgc                         >  1:3491188/1‑150 (MQ=255)
                                                                                                                             cATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCACGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCa                        <  1:4715530/150‑1 (MQ=255)
                                                                                                                             cATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCa                        <  1:1494727/150‑1 (MQ=255)
                                                                                                                              aTGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAg                       <  1:4932242/150‑1 (MQ=255)
                                                                                                                               tGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAgt                      >  2:3616689/1‑150 (MQ=255)
                                                                                                                                   ccTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAgttgtt                  <  2:4117246/150‑1 (MQ=255)
                                                                                                                                   ccTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAgttgtt                  <  2:627991/150‑1 (MQ=255)
                                                                                                                                   ccTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAgttgtt                  >  2:3852203/1‑150 (MQ=255)
                                                                                                                                     tCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGc                <  2:3614812/150‑1 (MQ=255)
                                                                                                                                     tCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGc                <  2:4729668/150‑1 (MQ=255)
                                                                                                                                       gtcgagCTGGCCGAGGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACATCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCaa              >  1:3176892/6‑149 (MQ=255)
                                                                                                                                       gtcgagCTGGCCGAGGGCCAGCCCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCAt              >  1:784218/6‑150 (MQ=255)
                                                                                                                                       aCCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCAt              >  1:5198144/1‑150 (MQ=255)
                                                                                                                                           agCTGGCCGAGGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATcccc          >  1:4446561/2‑150 (MQ=255)
                                                                                                                                           agCTGGCCGAGGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGGATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATcccc          >  1:4446729/2‑150 (MQ=255)
                                                                                                                                                  ccTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccct   >  1:824603/1‑150 (MQ=255)
                                                                                                                                                  ccTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccct   >  2:1960830/1‑150 (MQ=255)
                                                                                                                                                   cTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccctg  >  2:5403232/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCCCCTGGCCCTG  >  NC_002947/5348459‑5348755

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: