Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,657,957 A→C 11.4% intergenic (‑58/+103) PP_4961 ← / ← PP_4962 lipoprotein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,657,9570AC11.4% 528.9 / 4.8 212intergenic (‑58/+103)PP_4961/PP_4962lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (93/83);  new base C (3/29);  total (96/116)
Fisher's exact test for biased strand distribution p-value = 3.64e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 225 total aligned reads displayed.

GTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTT  >  NC_002947/5657808‑5658106
                                                                                                                                                     |                                                                                                                                                     
gTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCcaca                                                                                                                                                       <  2:3003217/150‑1 (MQ=255)
gTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCcaca                                                                                                                                                       <  1:3838425/150‑1 (MQ=255)
  ccAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAgg                                                                                                                                                     <  2:4414410/150‑1 (MQ=255)
        gTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGt                                                                                                                                               <  2:1920560/150‑1 (MQ=255)
        gTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGt                                                                                                                                               <  1:5079928/150‑1 (MQ=255)
        gTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGt                                                                                                                                               <  2:4050545/150‑1 (MQ=255)
         tCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTg                                                                                                                                              <  2:2754063/150‑1 (MQ=255)
         tCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTg                                                                                                                                              <  1:3033273/150‑1 (MQ=255)
         tCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTg                                                                                                                                              <  2:3727375/150‑1 (MQ=255)
           tGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCa                                                                                                                                            <  1:1692179/150‑1 (MQ=255)
           tGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCa                                                                                                                                            <  1:3802631/150‑1 (MQ=255)
           tGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCa                                                                                                                                            <  2:716077/150‑1 (MQ=255)
            gCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCAAAGGTACAGTGCAg                                                                                                                                           <  2:656703/150‑1 (MQ=255)
             ccGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAGGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGt                                                                                                                                          >  2:5234755/1‑150 (MQ=255)
             ccGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGt                                                                                                                                          >  1:4971699/1‑150 (MQ=255)
                 ggCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGTAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTctct                                                                                                                                      <  1:4209431/150‑1 (MQ=255)
                 ggCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTctct                                                                                                                                      <  2:1224022/150‑1 (MQ=255)
                   ccGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCCGTCTCtgt                                                                                                                                    >  2:3516814/1‑150 (MQ=255)
                   ccGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCtgt                                                                                                                                    <  1:4269927/150‑1 (MQ=255)
                     ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                                  <  1:3817362/150‑1 (MQ=255)
                     ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                                  >  1:3112690/1‑150 (MQ=255)
                       ttGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGa                                                                                                                                <  2:2135563/150‑1 (MQ=255)
                       ttGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGa                                                                                                                                <  2:3903982/150‑1 (MQ=255)
                       ttGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGa                                                                                                                                <  2:1807282/150‑1 (MQ=255)
                       ttGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGa                                                                                                                                <  1:4388255/150‑1 (MQ=255)
                         gcgcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCGCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGc                                                                                                                              <  2:717398/150‑1 (MQ=25)
                         gcgcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCGCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGc                                                                                                                              <  2:2443325/150‑1 (MQ=25)
                          cgcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCg                                                                                                                             <  1:4466059/150‑1 (MQ=25)
                           gcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCGCCCACAGGTACAGTGCAGTCTCTGTGGGAGCgg                                                                                                                            <  2:1419089/150‑1 (MQ=255)
                             cGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGc                                                                                                                          >  2:4925513/1‑150 (MQ=255)
                             cGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGc                                                                                                                          <  1:921308/150‑1 (MQ=255)
                               gCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGa                                                                                                                        <  1:4076337/150‑1 (MQ=255)
                                 tGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGc                                                                                                                      >  1:1782361/1‑150 (MQ=255)
                                      cagTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCcgcg                                                                                                                 >  1:1254174/1‑150 (MQ=255)
                                             cATCCGGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTGTCATTCCGTTTCCTCGATTGGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGGGAAGTCTCTGTGGGTGCGGGCGAGCCCGCGAAGGGcc                                                                                                          <  2:4872213/150‑1 (MQ=38)
                                               tCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCACTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGCGGGAGCGGGCGAGCCCGCGAAGGGCCGc                                                                                                        >  1:1572873/1‑150 (MQ=255)
                                                   gCCTGTGGCGCCAGGGCCGCCACGGCCAAAAGCGCTTTCATCCAGTTCCCTCGATTGGGGGGTGCTGTGCCGCCCTTCTCGGGCTCGCCCGCGCCCACCGGTACCGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAg                                                                                                    <  2:4272950/150‑1 (MQ=37)
                                                   gCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCATTTGCCGCGATTGGGGGCGGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACCGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAg                                                                                                    <  2:751988/150‑1 (MQ=38)
                                                    cctgggtCGCCCGGGCCGTCCCGGCCAACCGCGCTTGCCTTCCGTGGCCGCGCTTGGGGGCTGCGGGGCAGCCCTTCTCGGGCTCGCCGGCGCCCCCCGGTCCCGGCCCGTCGCGGGGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGc                                                                                                   <  2:4281460/145‑1 (MQ=255)
                                                    ccAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGc                                                                                                   >  2:1177699/1‑150 (MQ=255)
                                                        ggtCGCCCGGGCCGGCACGGCAAAACGCGCTTGCTGTCAGTTCCCGCGATTTGGGGGTGGGTTGCCGCCCTACTGGGGCTCGCCCGCGCCCACCGGTACTGGGCAGTGCCTGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCACAGCGGcc                                                                                               <  2:1881108/149‑1 (MQ=255)
                                                         gtCGCCTGGGCCGTCACGGCAAAAAGCGGTTTCATTCCGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTGGGGCTCGCCGGCTCCCACCGGTACTGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGccc                                                                                              <  2:3807174/150‑1 (MQ=255)
                                                            gCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCACGGGCTCGCCCGCACCCCCAGGTACAGTGGAGTCTCGGGGGGAGCGGGCGAGCCCCCGAAGGGCCGCACCCCGGCCCctg                                                                                           >  1:1481535/1‑149 (MQ=255)
                                                              cAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCcttct                                                                                         >  2:5063603/1‑150 (MQ=255)
                                                              cAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTCTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGGGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCcttct                                                                                         >  1:2834895/1‑150 (MQ=255)
                                                               aGGGCCGTCACGGCAAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTACTCGGGCACGCCCGCTCCCCCAGGTACAGTGCAGTCACTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCcttctt                                                                                        <  1:4686167/150‑1 (MQ=38)
                                                                 ggCCGTCACGGCCAAAAGCGCTTTCTGTCAGTTGCCTCGATTTGGGGCTGCTGTGCAGCCCTTCGCGGGCGCGCCCGCTCCCACCGGTCCAGTGCAGTCACTGGGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCAGccccaaatcga                                                                                      <  2:4230272/150‑8 (MQ=21)
                                                                    cGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACCGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGa                                                                                   <  2:5461141/150‑1 (MQ=255)
                                                                    cGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGa                                                                                   <  2:3674257/150‑1 (MQ=255)
                                                                     gTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGAc                                                                                  >  2:2523496/1‑150 (MQ=255)
                                                                              aaaaGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTGTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCa                                                                         >  2:3923687/1‑150 (MQ=255)
                                                                                    gcTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACCGCGGCCCCTTCTTCAAGACCACAAATCAACGTTc                                                                   >  1:5183464/1‑150 (MQ=255)
                                                                                    gcTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTACCACAGGTACACTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAACGGCCGCACCGCGGCCCCTTCTTCAAGACCACAAAACAACGtta                                                                   >  1:636881/1‑149 (MQ=39)
                                                                                      tttCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCCCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGCAGGGCCGCACAGCGGCCCCGTCATAAAGACCACAAAACAACCTTCAg                                                                 >  1:1382739/1‑150 (MQ=37)
                                                                                      ttcctttctGTTCCCTCTTTTGGGGGCTGCTGTGCCCCCCTACGCGGCCTCGCCCGCCCCCACTGTAACTGTGCTTCCTGGTGGGGTGGGGGCTAGCCCGCGAAGGGCCGCACAGCGGTCCCTTCTTCAAGACCACAAATCAACGTTCAg                                                                 <  1:3316812/145‑1 (MQ=255)
                                                                                          attcaGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCCAGACCACAAAACAACCTTCAGGTTg                                                             >  1:4926631/1‑150 (MQ=255)
                                                                                          attcaGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCGCCCACAGGTACAGTGCAGTCGCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAAACAACGTTAAGGTTg                                                             >  1:2606242/1‑150 (MQ=255)
                                                                                           ttcaGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGc                                                            <  2:1842081/150‑1 (MQ=255)
                                                                                             caGATTCCTCGAATTGGGGCTGCTGTGCAGCCCGACTCGGGCTCGCCCGCTCCCACCCGTACAGAGCCGGCTCACTGGGAGCGGGCGCCCCCGCGAAGGCCCGCACAGCGGCCCCTTCTTCAAGACCACAAAACAACGGTCAGGTAGCCt                                                          >  2:3655232/1‑150 (MQ=255)
                                                                                                    cTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAccc                                                   >  2:4682171/1‑150 (MQ=255)
                                                                                                    cTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCCCTCCCCCAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCACAGCGGCCCCATCTTCAAGACCACAAATAAACGTTCAGGTTTCCGATTAccc                                                   >  1:5270980/1‑150 (MQ=38)
                                                                                                        tttgggggCTGCTGTGCAGCCCTTCGCGGGCGCGCCCCCACCCCCCGGAACTGTGCAGTCTCTGGGGGTGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATc                                               <  1:3089997/146‑1 (MQ=255)
                                                                                                           tGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACCGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAg                                            <  2:2911818/150‑1 (MQ=255)
                                                                                                            ggggCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCATCTTCAAGACCACAAAACAACGTTCAGGTTGCCTATTACCCGATCAAgg                                           >  1:5364586/1‑150 (MQ=255)
                                                                                                             gggCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAggg                                          >  2:1231207/1‑150 (MQ=255)
                                                                                                               gctgctGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACCAATCAACGTTCAGGTTGCCTCTTACCCGATCAAGGGCa                                        >  2:4550838/1‑150 (MQ=255)
                                                                                                                ctgctgTTCTGCCCTTCTCGGGCTGGCCCGCTCCCACCTGTACAGTGCTTTCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCaa                                       <  2:1421164/150‑1 (MQ=38)
                                                                                                                 tgctgTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAAc                                      >  2:3706490/1‑150 (MQ=255)
                                                                                                                 tgctgTGCAGCCCATCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGGTCAGGTTTCCTATTACCCGATCAAGGGCAAc                                      >  1:4577921/1‑150 (MQ=255)
                                                                                                                  gctgTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAAcc                                     >  2:2236718/1‑150 (MQ=255)
                                                                                                                      tgCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGgc                                 <  1:5054661/150‑1 (MQ=255)
                                                                                                                           cccTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAACCGATCAAGGGCAACCCGGCGCCtt                            >  1:1252341/1‑150 (MQ=255)
                                                                                                                           cccTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGCTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCtt                            >  2:2951449/1‑150 (MQ=255)
                                                                                                                             cTTCTCGGGCTCGCCCGCTCCCCCAGGTACCGTGCAGTAGCTGTGGGAGCGGGCGGGCCCGCGAAGGGCCGCACAGCGGCCCCTGCTTGAAGACCACAAATCCACGTTCACGTTGCCGATTAACCGATCAATGGCAACCCGGCGCCTTgc                          >  1:3424865/1‑150 (MQ=255)
                                                                                                                              ttCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTgcg                         >  2:624676/1‑150 (MQ=255)
                                                                                                                              ttCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTgcg                         >  2:4149510/1‑150 (MQ=255)
                                                                                                                                 tcGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtg                      <  2:2128500/150‑1 (MQ=255)
                                                                                                                                 gcGGGCCCCCCCCCCCCCACCGTTACCGTCCAGTCTCTGTGGGCGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtg                      <  1:3289177/149‑1 (MQ=255)
                                                                                                                                  cGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGCGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtgg                     <  2:1677470/150‑1 (MQ=255)
                                                                                                                                    ggCTCGCCCGCCCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtggtg                   <  2:1557123/150‑1 (MQ=255)
                                                                                                                                     gCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtggtgg                  >  1:813621/1‑150 (MQ=255)
                                                                                                                                       tCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGCTCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTc                >  2:3628899/1‑150 (MQ=255)
                                                                                                                                       tCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTCCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTc                >  1:1418017/1‑150 (MQ=255)
                                                                                                                                          cccGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtc             >  2:680254/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCCCCCGGTACAGTGCCGTCTCTGGGGGAGCGGGCGCGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcg            <  2:2043236/150‑1 (MQ=39)
                                                                                                                                           ccGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcg            >  1:2272750/1‑150 (MQ=25)
                                                                                                                                           ccGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCGGGCGAGCCCGAGAAGGGCTGCACAGCAGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcg            >  1:2364974/1‑150 (MQ=21)
                                                                                                                                           ccGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCGGGCGAGCCCGAGAAGGGCTGCACAGCAGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcg            >  1:2958855/1‑150 (MQ=21)
                                                                                                                                             gcccccACCGGTCCTGTGCAGTCTCGGTGGGTGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcgtc          <  1:3500555/147‑1 (MQ=255)
                                                                                                                                             gCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcgtc          >  1:4598664/1‑150 (MQ=255)
                                                                                                                                               tCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTcgtcgt        >  1:1089617/1‑150 (MQ=255)
                                                                                                                                               tCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTcgtcgt        >  1:2145246/1‑150 (MQ=255)
                                                                                                                                                 ccACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTAc      >  1:4563798/1‑150 (MQ=255)
                                                                                                                                                 ccACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTAc      >  1:1862783/1‑150 (MQ=255)
                                                                                                                                                 ccACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTAc      >  1:133640/1‑150 (MQ=255)
                                                                                                                                                  cacaGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACt     >  1:4502341/1‑150 (MQ=255)
                                                                                                                                                   acaGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACCGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACtt    >  2:3251230/1‑150 (MQ=255)
                                                                                                                                                   acaGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACtt    >  2:245422/1‑150 (MQ=255)
                                                                                                                                                   acaGGTACAGTGCAGTCTCTATGGGAGCGGGCGAGCCCGCGAATGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACtt    >  1:728976/1‑150 (MQ=255)
                                                                                                                                                     aGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACtttt  >  2:1037258/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCCAGTTGTCTGCCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTT  >  NC_002947/5657808‑5658106

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