Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,051,108 T→A 6.2% 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 6.2% 292.1 / 8.8 96intergenic (‑63/+23)PP_1825/PP_1826hypothetical protein/isochorismatase family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (41/49);  new base A (3/3);  total (44/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

GGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGAC  >  NC_002947/2050964‑2051254
                                                                                                                                                |                                                                                                                                                  
ggccggccGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCgg                                                                                                                                               <  2:495385/150‑1 (MQ=255)
 gccggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCggg                                                                                                                                              >  1:1674389/1‑150 (MQ=255)
     gccgccAGGCGGGCGCCCCAGTTACGCGGGTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTGTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCggg                                                                                                                                          >  2:1875401/1‑149 (MQ=255)
      ccgccAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTc                                                                                                                                         <  2:1569628/150‑1 (MQ=255)
      ccgccAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTc                                                                                                                                         >  1:607577/1‑150 (MQ=255)
      ccgccAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGCCGGGCggaa                                                                                                                                         <  2:866805/150‑3 (MQ=255)
        gccAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGCCGGGCggaa                                                                                                                                         >  1:1384238/1‑146 (MQ=255)
        gccAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGCCGGGCggaa                                                                                                                                         <  2:1384238/148‑3 (MQ=255)
        gccAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCCTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCt                                                                                                                                       >  2:501811/1‑150 (MQ=255)
         ccAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCtt                                                                                                                                      <  1:67709/150‑2 (MQ=255)
         ccAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTg                                                                                                                                      >  1:864148/1‑150 (MQ=255)
            ggcgggcgCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCaa                                                                                                                                   <  1:500143/150‑1 (MQ=255)
                ggcgCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCg                                                                                                                               >  1:151271/1‑150 (MQ=255)
                     cccAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTTCGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAg                                                                                                                          >  2:1927614/1‑150 (MQ=255)
                     cccAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAg                                                                                                                          >  2:1886204/1‑150 (MQ=255)
                      ccAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGa                                                                                                                         >  1:1700482/1‑150 (MQ=255)
                        aGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATc                                                                                                                       <  2:880207/150‑1 (MQ=255)
                        aGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCCCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATc                                                                                                                       <  2:278451/150‑1 (MQ=255)
                           tACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCAGCAAGCCGCTCAGATCAgg                                                                                                                    <  1:1035215/150‑1 (MQ=255)
                                     gcgcgGGCGCAGGCGCTTGGCGTCGCAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTgg                                                                                                          <  2:1700482/150‑1 (MQ=255)
                                      cgcgGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTggg                                                                                                         <  2:607577/150‑1 (MQ=255)
                                       gcgGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGCgggg                                                                                                        <  2:1176585/150‑1 (MQ=255)
                                          ggCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGcc                                                                                                     <  1:1813018/150‑1 (MQ=255)
                                           gcgcAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCa                                                                                                    <  1:1420162/150‑1 (MQ=255)
                                              cAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGc                                                                                                 >  2:1223657/1‑150 (MQ=255)
                                                     ttGGCGTCGAAACTGCCTTTAAGCTTCATGCCGCTCTCCTCGGTGGATACCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCAGCCCGGCTGCCGTGCGGTCCTGCAAGCCGCTTAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGaa                                                                                          <  2:840679/150‑1 (MQ=255)
                                                             gAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAg                                                                                  >  2:367880/1‑150 (MQ=255)
                                                                aCTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCCCGCAAGCGAAGTTGTCAGCCa                                                                               >  1:166104/1‑150 (MQ=255)
                                                                  tGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCg                                                                                                                                            >  2:274384/1‑87 (MQ=255)
                                                                  tGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCg                                                                                                                                            <  1:274384/87‑1 (MQ=255)
                                                                  tGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGt                                                                                     <  1:30871/142‑1 (MQ=255)
                                                                  tGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGt                                                                                     >  2:30871/1‑142 (MQ=255)
                                                                  tGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATc                                                                             >  1:935776/1‑150 (MQ=255)
                                                                   gCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGtcatca                                                                            >  1:178616/1‑150 (MQ=255)
                                                                      tttAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCt                                                                                                                                       <  2:1743451/88‑1 (MQ=255)
                                                                      tttAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCt                                                                                                                                       >  1:1743451/1‑88 (MQ=255)
                                                                      tttAAGCTTCATTCCACTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACgg                                                                         <  2:1347997/150‑1 (MQ=255)
                                                                       ttAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGc                                                                        <  2:935776/150‑1 (MQ=255)
                                                                        tAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGcc                                                                       <  1:1223657/150‑1 (MQ=255)
                                                                           gCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCAtt                                                                    <  1:1746265/150‑1 (MQ=255)
                                                                            cTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTgggg                                                                                                        <  1:988287/113‑1 (MQ=255)
                                                                            cTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTgggg                                                                                                        >  2:988287/1‑113 (MQ=255)
                                                                                aTTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCa                                                               <  2:39759/150‑1 (MQ=255)
                                                                                    cGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCg                                                                                                                                            <  2:254739/69‑1 (MQ=255)
                                                                                    cGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCg                                                                                                                                            >  1:254739/1‑69 (MQ=255)
                                                                                    cGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTgg                                                           >  1:241560/1‑150 (MQ=255)
                                                                                     gCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGt                                                          <  1:113419/150‑1 (MQ=255)
                                                                                     gCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGt                                                          <  2:1252600/150‑1 (MQ=255)
                                                                                           ctcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCt                                                    >  1:462275/1‑150 (MQ=255)
                                                                                            tcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTg                                                   >  1:246162/1‑150 (MQ=255)
                                                                                             cGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTgc                                                  <  1:501811/150‑1 (MQ=255)
                                                                                              ggTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTgcg                                                 >  1:1286490/1‑150 (MQ=255)
                                                                                               gTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACgg                                                                         <  1:699726/125‑1 (MQ=255)
                                                                                               gTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACgg                                                                         >  2:699726/1‑125 (MQ=255)
                                                                                               gTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTgcgc                                                <  2:1816218/150‑1 (MQ=255)
                                                                                                  gATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGc                                             <  1:367880/150‑1 (MQ=255)
                                                                                                     atCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGt                                                                                        >  1:412873/1‑104 (MQ=255)
                                                                                                     atCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGt                                                                                        <  2:412873/104‑1 (MQ=255)
                                                                                                     atCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCgg                                          >  1:1778826/1‑150 (MQ=255)
                                                                                                      tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGaa                                                                                          >  2:343580/1‑101 (MQ=255)
                                                                                                      tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGaa                                                                                          <  1:343580/101‑1 (MQ=255)
                                                                                                      tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATc                                                                             >  1:1250616/1‑114 (MQ=255)
                                                                                                      tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATc                                                                             <  2:1250616/114‑1 (MQ=255)
                                                                                                      tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGa                                         >  2:1597155/1‑150 (MQ=255)
                                                                                                      tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGCCGGGCggaa                                                                                                                                         >  2:1858175/1‑52 (MQ=38)
                                                                                                      tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGCCGGGCggaa                                                                                                                                         <  1:1858175/54‑3 (MQ=255)
                                                                                                        cGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCAt                                                                              >  2:570876/1‑111 (MQ=255)
                                                                                                        cGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCAt                                                                              <  1:570876/111‑1 (MQ=255)
                                                                                                         ggAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATg                                      >  2:1334946/1‑150 (MQ=255)
                                                                                                         ggAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCCGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATg                                      >  2:1242716/1‑150 (MQ=255)
                                                                                                                  ccTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGc                                                                                                                                             >  2:1703144/1‑38 (MQ=255)
                                                                                                                  ccTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGc                                                                                                                                             <  1:1703144/38‑1 (MQ=255)
                                                                                                                  ccTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTggggg                                                                                                       >  2:815846/1‑76 (MQ=255)
                                                                                                                  ccTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTggggg                                                                                                       <  1:815846/76‑1 (MQ=255)
                                                                                                                   ctgctgctgCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCt                            <  2:1122215/150‑1 (MQ=255)
                                                                                                                       tgctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAg                                                                                         >  2:1862942/1‑85 (MQ=255)
                                                                                                                       tgctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAg                                                                                         <  1:1862942/85‑1 (MQ=255)
                                                                                                                       tgctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGaa                        <  2:1840290/150‑1 (MQ=255)
                                                                                                                        gctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAAt                       <  2:864148/150‑1 (MQ=255)
                                                                                                                         ctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGaa                                                                                          <  2:58233/82‑1 (MQ=255)
                                                                                                                         ctgcCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGaa                                                                                          >  1:58233/1‑82 (MQ=255)
                                                                                                                          tgcCCATCGGGTTGCGCCCGGCTGCCGGGCGCTAGCCCGTGGCTCATAGGGGCTGGTCGCGGTGGAGGCCACGCAAGCTAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGc                     <  1:1886204/150‑1 (MQ=255)
                                                                                                                            cccATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCa                   <  2:241560/150‑1 (MQ=255)
                                                                                                                            cccATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCa                   <  1:1811900/150‑1 (MQ=255)
                                                                                                                            cccATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCa                   >  1:1786868/1‑150 (MQ=255)
                                                                                                                            cccATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCa                   <  2:1786868/150‑1 (MQ=255)
                                                                                                                              cATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAgg                 >  1:1311757/1‑150 (MQ=255)
                                                                                                                                tCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTAGGGGCCACGCAAGCGAAGTTGTCAGCGATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCAAGGTc               >  2:465063/1‑150 (MQ=255)
                                                                                                                                   ggTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGa                                                        >  2:145884/1‑106 (MQ=255)
                                                                                                                                   ggTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGa                                                        <  1:145884/106‑1 (MQ=255)
                                                                                                                                   ggTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACg            >  2:698246/1‑150 (MQ=255)
                                                                                                                                   aggaccGCCCGGCAGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACg            >  1:1176908/5‑150 (MQ=255)
                                                                                                                                    gTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTAAACGTTGTTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACgg           <  1:1107326/150‑1 (MQ=255)
                                                                                                                                       ccGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTg        <  1:1927614/150‑1 (MQ=255)
                                                                                                                                         gCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGc      <  1:1242716/150‑1 (MQ=255)
                                                                                                                                             ggCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGAc  <  2:151271/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGAC  >  NC_002947/2050964‑2051254

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: