Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,018,851 T→G 18.3% intergenic (+51/+59) xpt → / ← PP_5266 xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,018,8510TG18.3% 358.1 / 73.0 162intergenic (+51/+59)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (56/73);  new base G (11/19);  total (67/95)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

Only 100 of 164 total aligned reads displayed.

AGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAG  >  NC_002947/6018703‑6018999
                                                                                                                                                    |                                                                                                                                                     
aGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTACCAAAGGTGCAGCACAGGtg                                                                                                                                                      >  2:3443776/1‑149 (MQ=255)
  ggcggcCGTGTCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCGCCCACACGTGCACCAAAGGgctc                                                                                                                                                    >  1:7857186/1‑150 (MQ=255)
     ggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGaa                                                                                                                                                 <  2:7429986/150‑1 (MQ=255)
      gcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAc                                                                                                                                                <  1:4769034/150‑1 (MQ=255)
      gcCGTGCCGAACTGGACAGCCAGGGTTACCGCGCTGAATCGCTGGCCCGGGTGAAGTCGCCGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGCCACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAc                                                                                                                                                <  2:1812462/150‑1 (MQ=255)
        cGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGc                                                                                                                                              >  2:1983520/1‑150 (MQ=255)
         gTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCt                                                                                                                                             <  1:2885629/150‑1 (MQ=255)
            ccGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCtgtg                                                                                                                                          <  1:3132819/150‑1 (MQ=255)
            ccGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCtgta                                                                                                                                          <  2:8305980/150‑2 (MQ=255)
             cGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                                         <  2:6794248/150‑1 (MQ=255)
             cGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGGCGCGAACGCAGTGc                                                                                                                                         >  2:7457811/1‑150 (MQ=255)
                  tGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                                    <  2:4935019/150‑1 (MQ=255)
                  tGGACAGCCAGGGTTACCGCATTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                                    <  1:7254929/150‑1 (MQ=255)
                     aCAGCCAGGGTTACCGCGTTGAATCGCTGGCCGGGGTGAAGTCGCTGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                                <  1:7271586/150‑1 (MQ=255)
                      cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                                <  1:3097443/150‑1 (MQ=255)
                      cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                                <  1:4311068/150‑1 (MQ=255)
                      cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                                <  1:6509583/150‑1 (MQ=25)
                      cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                                <  1:2562345/150‑1 (MQ=25)
                      cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                                <  2:3822404/150‑1 (MQ=25)
                      cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                                <  1:5819839/150‑1 (MQ=25)
                       aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtg                                                                                                                               <  2:7894284/150‑1 (MQ=255)
                       aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtg                                                                                                                               <  1:6363129/150‑1 (MQ=255)
                       aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgtg                                                                                                                               <  2:5203860/150‑1 (MQ=25)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                              <  1:1149973/150‑1 (MQ=255)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                              <  1:8516993/150‑1 (MQ=255)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                              <  2:6023575/150‑1 (MQ=255)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                              <  2:6537121/150‑1 (MQ=255)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                              <  2:6941503/150‑1 (MQ=255)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                              >  2:8406886/1‑150 (MQ=255)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                              <  1:2887672/150‑1 (MQ=255)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTgg                                                                                                                              <  2:7250803/150‑1 (MQ=25)
                          cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGa                                                                                                                            <  1:4990744/150‑1 (MQ=255)
                          cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGa                                                                                                                            <  1:5009108/150‑1 (MQ=21)
                              gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCgg                                                                                                                        <  1:6162561/150‑1 (MQ=255)
                              gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCgg                                                                                                                        <  2:2069516/150‑1 (MQ=255)
                              gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCgg                                                                                                                        <  2:7087837/150‑1 (MQ=255)
                              gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCgg                                                                                                                        <  1:1352229/150‑1 (MQ=21)
                              gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCgg                                                                                                                        <  1:6951597/150‑1 (MQ=21)
                                tACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGc                                                                                                                     <  2:4228390/150‑1 (MQ=21)
                                  ccGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGt                                                                                                                   <  2:7905776/150‑1 (MQ=255)
                                      attGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAACGGTGGGAGCGGGCGTGacc                                                                                                               >  2:4277746/2‑150 (MQ=255)
                                                  gCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGGACGCCCGCTCCCAAAGGTGGGGGACAGGTGTGGTAAGCTGTTCCGTTAATGTTGGTGGGGGGGTGGCCGCGAATGGgcgc                                                                                                   >  1:7973763/1‑150 (MQ=255)
                                                  gCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCCCGCCCGCGCCCACACGGGGCGCACACGGGTCGGACGCTGTGCCGGACCAGTGGGGGCGGGCCGGGCCCCCGAAGGCCGc                                                                                                   >  1:7108366/1‑150 (MQ=255)
                                                         tgacggcGCTGGATGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGTTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCtt                                                                                            <  1:8730256/144‑1 (MQ=255)
                                                          gAAGTCGCTGGACGGTGGTGTGGTCAGCTTCATCGAATAAACGCGCCCCCTGGTTCGCGGGCACGCCCGCTCCCACTGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGCGCGGGCGTGCCCGCGAATGGCCGCAACCCttt                                                                                           <  2:1181326/150‑1 (MQ=38)
                                                          gAAGTCGCTGGAAGGTGGTGTGGTCAGCTACAGCTAACAAACGCGCCACCGGCTTGGCGGGCACGCCCGCTCCCCCAGGTGCAGCACTGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCttt                                                                                           <  2:6613923/150‑1 (MQ=37)
                                                           acggcGCTGGATGGTGGTGTGGTCAGCTCCATCGAATAAACGCGCCACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTGCAGCACTGTTCTCGAAGGCTGTGCCGTACCTGTGGGCGCGGGCGTGCCCGCGAATGGCCGCAACCCtttt                                                                                          <  2:2546843/146‑1 (MQ=37)
                                                           aaGTCGCTGGACGGTGTTGTTGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGGGGGCACGCCCGCTCCCCCCGGTGCAGCACTGGAGCGGAACGCTGTGGCTTACGTGTGTGCGCGGGAGTTCCGGCTAATGGCCGCAACCCtttt                                                                                          <  2:4730641/150‑1 (MQ=255)
                                                                 cTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCACCCACACGTGCAGCACAGGTCTCCAAAGCTGTGCCGTACCTCAGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCAATTATCTgg                                                                                    >  2:3591232/1‑150 (MQ=38)
                                                                                gtCAGCTTCATCGAATAAACGCTACACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTGCTGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGATCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCg                                                                     <  2:7334424/150‑1 (MQ=255)
                                                                                          tCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaacaacaa                                                           <  2:5789192/150‑1 (MQ=255)
                                                                                          tCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCGCCCCCCCGGGCCGCCCCGGTCTCGCACCCGGTGCCGGACCCGTGGGGGCGGGCGGGCCCGCGCATGGCCGCCACCCGATTATAAGGTTTGCGCCCGCCCCCCccaccaaca                                                           >  2:4808683/1‑145 (MQ=255)
                                                                                             aaTAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCCGCCCCGCCAACAACAAgcg                                                        >  2:8743801/1‑150 (MQ=255)
                                                                                              aTAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAgcgc                                                       <  2:4549826/150‑1 (MQ=255)
                                                                                                 aaCGCGACACCTGCTTCGCGGGCACGCCCGCGCCCACACGTGGAGCAAAAGTCTCGCACCCGGTGCCGTACCGGGGGGTCCGGGCGTGCCCCCGAATTGCCCCAAACCAATTATAAGTTTTGCGCCCGCCCCCCcaacaaaaacaggtgg                                                    >  2:1616464/1‑143 (MQ=255)
                                                                                                   cgcgACACCTGCTTCGCGGGCACCCCCCCGCCCCCACGGGCCGCCCACGTCTCGAAACCGGGGGCCTACCAGTGGGAGCCGGCGGGCCCCCGGATGGCCGCAACCCGTTGAGCTGGTTGCCGCCACCCCCCCCACCCACCAgccccggga                                                  >  2:5435811/1‑143 (MQ=255)
                                                                                                     ctacacCTGCTTCGCGGGCACGCCCGCTCCCACTGGTGCAGAACAGGTATCTACCGGTGTGCCGTACATGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTaaa                                                <  1:7889315/148‑1 (MQ=38)
                                                                                                     cgACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTaaa                                                <  2:3255801/150‑1 (MQ=255)
                                                                                                     cgACACCGGCTTCGGGGGCACGCCCGCTCCCACCGGTGCTGCACTGGTCTCGATCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTaaa                                                <  1:4250019/150‑1 (MQ=38)
                                                                                                         acCTGTGTCGGGGGCCCGCCCGCCCCCACTGGTGCCGCACTGGTCGCGACCGCTTCGCGGGACGTGTGGTCGGGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttc                                            <  1:123565/150‑1 (MQ=255)
                                                                                                           cTGCTTCGCGGGCACGCCCGCTCCCACAGGTGGAGCACAGGTCTCGAAAGCTGTGCCGTACCGGTGGGAGCGGGCGTGGCCGCGAAAGGCCGCAACCCATTTATCTGGTTTCCGCCACCCCCGCCCACAACAAGCCCTGGGAAAGttctt                                          >  2:5586301/1‑150 (MQ=37)
                                                                                                           cTGCTTCGCGGGCACCCCCGCGCCCACAGGTGCAGCACAGGTCACGAACCCGGTGCCGTACCAGGGGGAGCGGGCGTGCCCGCGAATTGCCGCAACCCGTTTATCTGGTTGCCGCCAGCCCCGCCAACAACCAGCGCGGGGAAAGttctt                                          >  2:5586309/1‑150 (MQ=37)
                                                                                                               ttCGCGGGCACGCCCGATCCCACAGTTGCACCACAGCTCTATAACGCTGTGCCCTACCTGTGGTAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCttt                                      <  1:3443776/150‑1 (MQ=255)
                                                                                                               ggcgcgGGCCCGCCCGCCCCCCCCGGGGCCGAACCGGGCCGGACCGCGGGGCGGTCCGGGGGGGCGCGGGCGGGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCttt                                      <  1:3371891/148‑1 (MQ=255)
                                                                                                               ctCGCGGGCACGCCCGCTCCCACAGGAGCAGCACAGGACACGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCttt                                      <  2:545173/149‑1 (MQ=255)
                                                                                                                 ggggggCACGCCCGCTCCCACAGTTGCTGAACTGTTCACGAACGCTGTGCCTTCCCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTaa                                    <  1:8406886/147‑1 (MQ=38)
                                                                                                                 cgcgGGCACGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTaa                                    <  2:7064118/150‑1 (MQ=255)
                                                                                                                  gcgGGCACGCCCGCTCCCACAGGTACGGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                                    >  1:903918/1‑54 (MQ=18)
                                                                                                                   ggggCACGCCCGCTCCCCCGGGTGCCGCACCGGTCGCGCCCGGTGTGCGGTACCTGGGGGCGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGc                                  <  2:2905738/149‑1 (MQ=255)
                                                                                                                    gggCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCAGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTAACTGGTTGCCGCCAGCCCCGCCAACAACAAGCCCTGGTAAAGTTCATCTTTAAGcc                                 >  1:2729198/1‑150 (MQ=255)
                                                                                                                        aCGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTg                             <  2:3334743/150‑1 (MQ=255)
                                                                                                                            ccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGc                         >  2:663270/1‑150 (MQ=255)
                                                                                                                            ccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGc                         >  1:4667333/1‑150 (MQ=255)
                                                                                                                            ccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGc                         >  1:4444618/1‑150 (MQ=255)
                                                                                                                            ccGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGc                         >  2:8239068/1‑150 (MQ=21)
                                                                                                                              gCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTg                       <  2:7173013/150‑1 (MQ=255)
                                                                                                                                tCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGTTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc                     >  1:405889/1‑150 (MQ=255)
                                                                                                                                tCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc                     >  1:7356578/1‑150 (MQ=255)
                                                                                                                                ccccACCGGTGCAGCACCGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc                     <  1:7337405/149‑1 (MQ=255)
                                                                                                                                  ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                     <  1:1582824/116‑1 (MQ=255)
                                                                                                                                  ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                     >  2:1582824/1‑116 (MQ=255)
                                                                                                                                  ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                   >  1:6695682/1‑150 (MQ=255)
                                                                                                                                  ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                   >  1:7511516/1‑150 (MQ=255)
                                                                                                                                  ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                   >  1:888008/1‑150 (MQ=255)
                                                                                                                                  ccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                   >  2:2192786/1‑150 (MQ=25)
                                                                                                                                   cacaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCa                  >  1:536618/1‑150 (MQ=255)
                                                                                                                                   cacaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCa                  >  2:6803507/1‑150 (MQ=25)
                                                                                                                                   cacaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCa                  >  1:73436/1‑150 (MQ=25)
                                                                                                                                   cacaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCa                  >  1:4365102/1‑150 (MQ=25)
                                                                                                                                    gcaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                 >  2:3352082/2‑150 (MQ=255)
                                                                                                                                    acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                 >  1:5515652/1‑150 (MQ=255)
                                                                                                                                    acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                 >  2:7667651/1‑150 (MQ=255)
                                                                                                                                    acaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                 >  2:2999666/1‑150 (MQ=25)
                                                                                                                                    acaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                 >  1:459806/1‑150 (MQ=25)
                                                                                                                                        gTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGc             >  1:2228230/1‑150 (MQ=255)
                                                                                                                                        gTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTGGTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGc             >  1:6033250/1‑150 (MQ=255)
                                                                                                                                            agcaCAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATcc         >  1:608452/1‑150 (MQ=255)
                                                                                                                                                 aGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGcccc    >  1:4395670/1‑150 (MQ=255)
                                                                                                                                                 aGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGcccc    >  2:223009/1‑150 (MQ=255)
                                                                                                                                                 aGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGcccc    >  2:6283381/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAG  >  NC_002947/6018703‑6018999

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: