Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,348,607 T→A 5.3% L20H (CTC→CAC) ‡ PP_4703 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,348,6070TA 5.3% 509.3 / 11.8 151L20H (CTC→CAC) ‡PP_4703hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (70/73);  new base A (6/2);  total (76/75)
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.88e-01

Only 100 of 157 total aligned reads displayed.

CCGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCCCCTGGCCCTG  >  NC_002947/5348458‑5348755
                                                                                                                                                     |                                                                                                                                                    
ccGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCt                                                                                                                                                      <  2:3604355/150‑1 (MQ=255)
ccGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCt                                                                                                                                                      <  2:3604564/150‑1 (MQ=255)
 cGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTc                                                                                                                                                     <  1:1781191/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)
       gAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCGAGGGCCAg                                                                                                                                               <  1:2465116/150‑1 (MQ=255)
       gAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAg                                                                                                                                               <  1:1756296/150‑1 (MQ=255)
               tACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCt                                                                                                                                       <  2:4810188/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:5434415/1‑150 (MQ=255)
                    ggCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAAt                                                                                                                                  >  1:1844863/1‑150 (MQ=255)
                       cTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGc                                                                                                                               <  2:3004864/150‑1 (MQ=255)
                       cTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGc                                                                                                                               <  2:5086033/150‑1 (MQ=255)
                            cGACAAACGCCTCAACCAGTAACCGGCGGCCACAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAg                                                                                                                          <  1:858015/150‑1 (MQ=255)
                             gACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAgg                                                                                                                         >  1:3451160/1‑150 (MQ=255)
                              aCAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCCGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAggg                                                                                                                        >  2:3146802/1‑150 (MQ=255)
                               cAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCTTGGGc                                                                                                                       <  1:264340/150‑1 (MQ=255)
                               cAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGc                                                                                                                       <  2:2852643/150‑1 (MQ=255)
                               cAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGc                                                                                                                       <  1:953479/150‑1 (MQ=255)
                                aaaCGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCt                                                                                                                                                      <  2:2203461/118‑1 (MQ=255)
                                aaaCGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCg                                                                                                                      >  2:2424758/1‑150 (MQ=255)
                                aaaCGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCg                                                                                                                      >  2:1528104/1‑150 (MQ=255)
                                aaaCGCCCCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTg                                                                                                                                      >  1:2865632/1‑134 (MQ=255)
                                 aaCGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGa                                                                                                                     <  1:1744942/150‑1 (MQ=255)
                                  aCGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGAc                                                                                                                    >  1:1807573/1‑150 (MQ=255)
                                      cTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCTTGGGCGattcca                                                                                                                <  1:2283077/150‑6 (MQ=255)
                                      cTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCg                                                                                                                >  2:2361807/1‑150 (MQ=255)
                                       tCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCTTGGGCGattccag                                                                                                               <  1:1336532/150‑7 (MQ=255)
                                       tCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCgg                                                                                                               >  2:802847/1‑150 (MQ=255)
                                        cAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGc                                                                                                              <  1:2580988/150‑1 (MQ=255)
                                               tAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCg                                                                                                       <  1:802847/150‑1 (MQ=255)
                                                   cggcggCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGATACCCATGTGTTGAGCCTCTCCGGGCTGGCGCTCGGCCAGCTCGACCTGGACTCGCCCAAGGGCGACCCCGGCCTGTTCGGccc                                                                                                   <  2:2796735/150‑1 (MQ=255)
                                                   cggcggCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGccc                                                                                                   <  1:311858/150‑1 (MQ=255)
                                                          cGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCTTGGGCGat                                                                                                                    >  2:928596/1‑125 (MQ=255)
                                                          cGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCt                                                                                                            <  1:5307155/134‑1 (MQ=255)
                                                          cGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGt                                                                                            >  2:4366182/1‑150 (MQ=255)
                                                           gCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGcc                                                                                                                                                 <  1:3810596/96‑1 (MQ=255)
                                                           gCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTa                                                                                           >  1:2576585/1‑150 (MQ=255)
                                                                tCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGc                                                                                      >  1:1545681/1‑150 (MQ=255)
                                                                tCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGc                                                                                      <  2:1427752/150‑1 (MQ=255)
                                                                  tGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTg                                                                                    >  1:5183538/1‑150 (MQ=255)
                                                                   ggCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTgg                                                                                   <  2:1385466/150‑1 (MQ=255)
                                                                        cccTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGAGCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGt                                                                              <  2:588466/150‑1 (MQ=255)
                                                                          cTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTCTCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGACTCGCCCAAGGGCGACCCCGGCCTGGTCGGCCCGCACAGGCTTAGCTGGCGCGTcc                                                                            >  2:3004246/1‑150 (MQ=255)
                                                                          cTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTcc                                                                            >  2:1493474/1‑150 (MQ=255)
                                                                              aCAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGAcccg                                                                                                                                      <  1:2938091/88‑3 (MQ=255)
                                                                              aCAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATgg                                                                        <  1:1481092/150‑1 (MQ=255)
                                                                              aCAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATgg                                                                        >  2:1481092/1‑150 (MQ=255)
                                                                               cagcagGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCGGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGc                                                                       <  2:3648073/150‑1 (MQ=255)
                                                                                 gcagGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGa                                                                     <  2:4070579/150‑1 (MQ=255)
                                                                                  cagGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGAc                                                                    <  2:4981776/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)
                                                                                              gAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACGTCCCGAGCCTg                                                        >  2:2869211/1‑150 (MQ=255)
                                                                                                 gCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTg                                                     <  2:2878387/150‑1 (MQ=255)
                                                                                                 gCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCATCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTg                                                     >  2:1158506/1‑150 (MQ=255)
                                                                                                   ccTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCt                                                                                                                                       <  2:2146417/66‑1 (MQ=255)
                                                                                                   ccTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCt                                                                                                                                       >  1:2146417/1‑66 (MQ=255)
                                                                                                   ccTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTACCGACCagcatgct                                                   <  2:1863131/150‑8 (MQ=255)
                                                                                                    cTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTgggc                                                  >  2:4962330/1‑148 (MQ=255)
                                                                                                       cgtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGt                                               <  1:4968091/150‑1 (MQ=255)
                                                                                                        gtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTg                                              <  1:1026951/150‑1 (MQ=255)
                                                                                                        gtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTg                                              <  1:2525462/150‑1 (MQ=255)
                                                                                                        gtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTg                                              >  2:5152537/1‑150 (MQ=255)
                                                                                                        gtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTg                                              >  2:5003469/1‑150 (MQ=255)
                                                                                                        gtcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTg                                              <  1:3767531/150‑1 (MQ=255)
                                                                                                         tcgCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTgg                                             >  1:587050/1‑150 (MQ=255)
                                                                                                                  cGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAgcgc                                    <  1:3730443/150‑1 (MQ=255)
                                                                                                                   gAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGc                                                                                                              >  2:1658747/1‑75 (MQ=255)
                                                                                                                       cccAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGa                               <  1:795286/150‑1 (MQ=255)
                                                                                                                        ccAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGAt                              >  2:4766974/1‑150 (MQ=255)
                                                                                                                          aGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGAtgc                            <  2:3951455/150‑1 (MQ=255)
                                                                                                                           ggTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGAtgct                           <  2:2015271/150‑1 (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)
                                                                                                                               aTGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAg                       <  1:4932256/150‑1 (MQ=255)
                                                                                                                                tGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAgt                      >  2:3616689/1‑150 (MQ=255)
                                                                                                                                 gAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAgtt                     >  2:1100896/1‑150 (MQ=255)
                                                                                                                                    ccTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAgttgtt                  >  2:3852203/1‑150 (MQ=255)
                                                                                                                                    ccTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAgttgtt                  <  2:4117246/150‑1 (MQ=255)
                                                                                                                                    ccTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAgttgtt                  <  2:627991/150‑1 (MQ=255)
                                                                                                                                      tCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGc                <  1:5387235/150‑1 (MQ=255)
                                                                                                                                        gtcgagCTGGCCGAGGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACATCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCaa              >  1:3176892/6‑149 (MQ=255)
                                                                                                                                        aCCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCAt              >  1:5198144/1‑150 (MQ=255)
                                                                                                                                        aCCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCAt              >  1:5286532/1‑150 (MQ=255)
                                                                                                                                            agCTGGCCGAGGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATcccc          >  1:4446561/2‑150 (MQ=255)
                                                                                                                                            agCTGGCCGAGGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATcccc          >  1:4446657/2‑150 (MQ=255)
                                                                                                                                               tGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctgg       <  1:473902/150‑1 (MQ=255)
                                                                                                                                                   ccTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccct   >  1:824603/1‑150 (MQ=255)
                                                                                                                                                   ccTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccct   >  2:1960830/1‑150 (MQ=255)
                                                                                                                                                    cTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccctg  >  2:2730888/1‑150 (MQ=255)
                                                                                                                                                    cTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCccctggccctg  >  2:5403232/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGACCCGAGCGTGGTACAGGGCCTGATCGACAAACGCCTCAACCAGTAACCGGCGGCCGCAAATCTGGCCGCCCTGTACAGCAGGTCGCCATGGAAGCCCTTCGTCGCCGTATCGAAACCCAGGTCATGAGCCTCACCGGGCTGGCCCTCGGCCAGCTCGACCTGGAATCGCCCAAGGGCGACCCCGGCCTGTTCGGCCCGCACAGTATTAGCTGGCGCGTCCATGGCGACTTCCCGAGCATGCTGGTCGGTGGCATCAGCGCGCTGATGCTGCAGTTGTTGCATCCCCTGGCCCTG  >  NC_002947/5348458‑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: