Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,051,108 T→A 5.1% intergenic (‑63/+23) PP_1825 ← / ← PP_1826 hypothetical protein/isochorismatase family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,051,1080TA 5.1% 409.5 / 11.5 136intergenic (‑63/+23)PP_1825/PP_1826hypothetical protein/isochorismatase family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (57/72);  new base A (3/4);  total (60/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

Only 100 of 137 total aligned reads displayed.

CGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACG  >  NC_002947/2050959‑2051255
                                                                                                                                                     |                                                                                                                                                   
cGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCt                                                                                                                                                     <  2:1581179/150‑1 (MQ=255)
   cAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGcc                                                                                                                                                  <  1:1604283/150‑1 (MQ=255)
      gccggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCggg                                                                                                                                               >  2:2532777/1‑150 (MQ=255)
       ccggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATGCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGc                                                                                                                                              <  2:1431987/150‑1 (MQ=255)
        cggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCg                                                                                                                                             >  2:903839/1‑150 (MQ=255)
        cggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCg                                                                                                                                             >  2:1219218/1‑150 (MQ=255)
        cggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCg                                                                                                                                             >  2:1539626/1‑150 (MQ=255)
        cggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCg                                                                                                                                             <  2:2505167/150‑1 (MQ=255)
                 ggcgggcgCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCaa                                                                                                                                    <  1:1608781/150‑1 (MQ=255)
                   cgggcgCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGc                                                                                                                                  <  1:1003641/150‑1 (MQ=255)
                     ggcgCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCg                                                                                                                                <  1:515010/150‑1 (MQ=255)
                       cgcCCCAGTCACGGGGTTCGCGCGGGCGAAGGCGCTTGGCGTCGAACCTGCCGTTAAGCTTCATTCCCCTCTCCTCGGTGTATATCCGGAATTCGACCTGCTGCTGCCCAGCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCt                                                                                                                              <  2:1343504/150‑1 (MQ=255)
                          cccAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGc                                                                                                                                                   <  2:2895527/126‑1 (MQ=255)
                          cccAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAg                                                                                                                           >  2:727503/1‑150 (MQ=255)
                           ccAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCAGCTCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGa                                                                                                                          <  2:1964484/150‑1 (MQ=255)
                            cAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGAt                                                                                                                         <  1:1812626/150‑1 (MQ=255)
                                           cgcgAGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTggg                                                                                                          <  2:2318737/150‑1 (MQ=255)
                                            gcgGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTgggg                                                                                                         <  1:2103072/150‑1 (MQ=255)
                                            gcgGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTgggg                                                                                                         <  2:2421995/150‑1 (MQ=255)
                                             cgGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTggggg                                                                                                        <  2:2856275/150‑1 (MQ=255)
                                               ggCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGcc                                                                                                      <  1:1219218/150‑1 (MQ=255)
                                                 cgcAGGCGCTTGGTGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTTCCGGGCGGTCCTGCAAGCTGCTCAGATCAGGCTGGCGGTGGGGGCCAc                                                                                                    <  1:717748/150‑1 (MQ=255)
                                                     ggCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAgg                                                                                                                     <  1:1551318/129‑1 (MQ=255)
                                                     ggCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGCCTGAGCTGCAAGCCGCTCAGATCAGCCTGGCGGTGGGGGCCACGCaa                                                                                                <  1:2239545/150‑1 (MQ=255)
                                                        gcTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCg                                                                                             >  1:2002225/1‑150 (MQ=255)
                                                        gcTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCg                                                                                             >  2:526310/1‑150 (MQ=255)
                                                         cTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGa                                                                                            <  2:1075451/150‑1 (MQ=255)
                                                          ttGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGaa                                                                                           <  2:293879/150‑1 (MQ=255)
                                                           tGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAg                                                                                          <  2:227281/150‑1 (MQ=255)
                                                           tGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAg                                                                                          <  2:482928/150‑1 (MQ=255)
                                                                 cGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCa                                                                                    <  2:1018942/150‑1 (MQ=255)
                                                                 cGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCATGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCa                                                                                    <  2:431344/150‑1 (MQ=255)
                                                                 cGAAACGGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCa                                                                                    <  1:1973773/150‑1 (MQ=255)
                                                                  gAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTg                                                                                                                                       <  2:2030739/98‑1 (MQ=255)
                                                                   aaaCTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTAAGc                                                                                  <  1:2941605/150‑1 (MQ=255)
                                                                            ttAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGCTCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGc                                                                         <  2:1964536/150‑1 (MQ=255)
                                                                               aGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCAt                                                                      >  1:148399/1‑150 (MQ=255)
                                                                                 cctccttCCGCTCTCCTCGCTGGATATCCGGAATTCGACGTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCATGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCAttt                                                                    <  2:1208936/145‑1 (MQ=255)
                                                                                 cTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACgg                                                                          >  1:427560/1‑144 (MQ=255)
                                                                                 cTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACgg                                                                          <  2:427560/144‑1 (MQ=255)
                                                                                 cTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCAttt                                                                    <  2:2680739/150‑1 (MQ=255)
                                                                                  ttCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTc                                                                   >  1:2141788/1‑150 (MQ=255)
                                                                                      ttCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCAc                                                               >  1:552538/1‑150 (MQ=255)
                                                                                      ttCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCAc                                                               >  1:1109170/1‑150 (MQ=255)
                                                                                      ttCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCAc                                                               >  2:1029980/1‑150 (MQ=255)
                                                                                        ccGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTg                                                             <  1:727503/150‑1 (MQ=255)
                                                                                         cGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGc                                                                                              <  1:1041364/116‑1 (MQ=255)
                                                                                         cGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTgg                                                            >  2:323042/1‑150 (MQ=255)
                                                                                         cGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTgg                                                            >  2:500757/1‑150 (MQ=255)
                                                                                          gCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGt                                                           <  1:2889927/150‑1 (MQ=255)
                                                                                          gCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGt                                                           <  2:2086536/150‑1 (MQ=255)
                                                                                           ctctCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCAtt                                                                     <  1:2562667/139‑1 (MQ=255)
                                                                                             ctcctcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGCCGGGCggaa                                                                                                                                          >  2:3041567/1‑66 (MQ=255)
                                                                                             ctcctcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGCCGGGCggaa                                                                                                                                          <  1:3041567/68‑3 (MQ=255)
                                                                                              tcctcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAAt                                                       >  1:2770546/1‑150 (MQ=255)
                                                                                              tcctcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAAt                                                       <  1:1029980/150‑1 (MQ=255)
                                                                                                ctcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCaa                                                                                                <  1:1110067/107‑1 (MQ=255)
                                                                                                ctcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGAAACGCaa                                                                                                >  2:1110067/1‑107 (MQ=255)
                                                                                                 tcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCTATCACGGCCAATTCGCACCTGTGAATCTg                                                    >  2:2191553/1‑150 (MQ=255)
                                                                                                  cGGTGGATATCCGGAAGTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTgc                                                   <  2:1326321/150‑1 (MQ=255)
                                                                                                     tGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGcgcg                                                >  2:174266/1‑150 (MQ=255)
                                                                                                           tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTc                                                                                     >  1:2569189/1‑107 (MQ=255)
                                                                                                           tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTc                                                                                     <  2:2569189/107‑1 (MQ=255)
                                                                                                           tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCAc                                                               <  1:1902152/129‑1 (MQ=255)
                                                                                                           tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCAc                                                               >  2:1902152/1‑129 (MQ=255)
                                                                                                           tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTgg                                                            >  2:761602/1‑132 (MQ=255)
                                                                                                           tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTgc                                                   >  2:2087997/1‑141 (MQ=255)
                                                                                                           tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTgc                                                   <  1:2087997/141‑1 (MQ=255)
                                                                                                           tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGa                                          >  2:881821/1‑150 (MQ=255)
                                                                                                           tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGa                                          >  2:2734932/1‑150 (MQ=255)
                                                                                                            ccGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGaa                                         <  2:3014532/150‑1 (MQ=255)
                                                                                                              ggAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGaa                                                                                           <  1:866622/98‑1 (MQ=255)
                                                                                                              ggAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGATCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCa                                                                                    >  2:1507865/1‑105 (MQ=255)
                                                                                                                  ttCGACCTGCTGCTCCCCATCGTGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTTAGATCAGGCTGGCGGTGGGTGCCACGCAAGCCAAGTTGTCAGCAATCACGGTCATTTCGCACTAGTGAATCTGCGCGGCCGGAATGACGc                                   <  2:217844/150‑1 (MQ=255)
                                                                                                                      aCCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCCGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGaa                                         <  1:2329838/140‑1 (MQ=255)
                                                                                                                        ctgctgctgCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCt                             <  2:2002225/150‑1 (MQ=255)
                                                                                                                         tgctgctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGACGTTGTCAGCCCTCACGGCCATTTCCCAGCGGGGGATCTGCGGGGTGTTAATGACGCCGTccgt                            <  2:1572738/150‑3 (MQ=255)
                                                                                                                         tgctgctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCtt                            >  1:114521/1‑150 (MQ=255)
                                                                                                                          gctgctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTc                                                                                     <  1:751714/92‑1 (MQ=255)
                                                                                                                          gctgctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTc                                                                                     >  2:751714/1‑92 (MQ=255)
                                                                                                                               tgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCTATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGc                      <  1:2191553/150‑1 (MQ=255)
                                                                                                                                gcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGTGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCCCTGCTGAAGCGGCGCGGCCGGAATGACGCCGTCCTTGAATGcc                     >  2:2254858/1‑150 (MQ=255)
                                                                                                                                gcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGt                                                                                      <  1:1219198/85‑1 (MQ=255)
                                                                                                                                gcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGt                                                                                      >  2:1219198/1‑85 (MQ=255)
                                                                                                                                 cccATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCAc                                                               >  2:1628401/1‑107 (MQ=255)
                                                                                                                                 cccAGCGGGTTCCGCCCGGCTGCCGGGCGGTCGTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCa                    <  2:1075016/150‑1 (MQ=255)
                                                                                                                                   cATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAgg                  >  1:2605196/1‑150 (MQ=255)
                                                                                                                                     tCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTc                >  1:354909/1‑150 (MQ=255)
                                                                                                                                       gggTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCAc              >  1:2890491/1‑150 (MQ=255)
                                                                                                                                        aggaccGCCCGGCAGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTcc                              >  1:2901176/5‑133 (MQ=255)
                                                                                                                                        aggaccGCCCGGCAGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTcc                              <  2:2901176/129‑1 (MQ=255)
                                                                                                                                         gTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACgg            <  1:2774685/150‑1 (MQ=255)
                                                                                                                                         gTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACgg            <  1:903839/150‑1 (MQ=255)
                                                                                                                                          ttCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCAtt                                                                     >  1:1317096/1‑92 (MQ=255)
                                                                                                                                          gaccGCCCGGCAGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAAt                                        >  1:997737/3‑121 (MQ=255)
                                                                                                                                            ccGCCCGGCAGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTg         <  1:1385874/150‑1 (MQ=255)
                                                                                                                                               cccGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGt                                <  2:286364/124‑1 (MQ=255)
                                                                                                                                               cccGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGt                                >  1:286364/1‑124 (MQ=255)
                                                                                                                                                 cGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACgg            <  1:1704877/142‑1 (MQ=255)
                                                                                                                                                  ggCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGAc   <  1:548129/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACG  >  NC_002947/2050959‑2051255

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