Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,015,814 G→A 10.3% intergenic (+41/+69) 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,015,8140GA10.3% 327.3 / 23.3 118intergenic (+41/+69)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (40/63);  new base A (6/6);  total (49/69)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.34e-01

Only 100 of 128 total aligned reads displayed.

AGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCAT  >  NC_002947/6015667‑6015963
                                                                                                                                                   |                                                                                                                                                     
 gTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTgcag                                                                                                                                                    <  1:3326728/150‑1 (MQ=255)
  tCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTgcagc                                                                                                                                                   >  1:3352801/1‑150 (MQ=255)
   cGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTgcagca                                                                                                                                                  <  2:132891/150‑1 (MQ=255)
    gTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGcac                                                                                                                                                 <  1:1667599/150‑1 (MQ=255)
     ttCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGcaca                                                                                                                                                <  2:3928733/150‑1 (MQ=255)
      tCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAg                                                                                                                                               <  1:1181318/150‑1 (MQ=255)
         aGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtc                                                                                                                                            >  1:3915560/1‑150 (MQ=255)
          gggCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGATACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtct                                                                                                                                           <  1:3889755/150‑1 (MQ=255)
          gggCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtct                                                                                                                                           <  2:3729185/150‑1 (MQ=255)
           ggcggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtctc                                                                                                                                          <  2:2370310/150‑1 (MQ=255)
           ggcggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtctc                                                                                                                                          <  2:4217853/150‑1 (MQ=255)
           ggcggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtctc                                                                                                                                          <  1:4275427/150‑1 (MQ=255)
            gcggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCg                                                                                                                                         <  1:52290/150‑1 (MQ=255)
               gcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAc                                                                                                                                      <  1:1279114/150‑1 (MQ=255)
                ccGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCGCCCCCAGGTGCAGCACAGGGGTCGAACg                                                                                                                                     >  2:1165158/1‑150 (MQ=255)
                    gCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCtgt                                                                                                                                 <  2:802464/150‑1 (MQ=255)
                    gCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCtgt                                                                                                                                 <  2:3915560/150‑1 (MQ=255)
                    gCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGACACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCtgt                                                                                                                                 <  1:1817999/150‑1 (MQ=255)
                     ccGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCtgtg                                                                                                                                <  1:1438957/150‑1 (MQ=255)
                      cGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                               <  1:3252970/150‑1 (MQ=255)
                      cGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                               <  1:3030874/150‑1 (MQ=255)
                      cGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCGGCTTCGCGGGCACGCCCGCTCCCACAGGGGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                               >  2:2606409/1‑150 (MQ=255)
                         aCTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGt                                                                                                                            <  2:1245326/150‑1 (MQ=255)
                          cTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTa                                                                                                                           <  1:4401235/150‑1 (MQ=255)
                          cTGGACAGCCAGGGTTACCGCGTTGAAGCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTa                                                                                                                           <  1:2860912/150‑1 (MQ=255)
                           tGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                          <  2:2992162/150‑1 (MQ=255)
                           tGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                          <  1:1143433/150‑1 (MQ=255)
                               cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGGACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                      <  2:1486978/150‑1 (MQ=25)
                               cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                      <  2:3284492/150‑1 (MQ=25)
                               cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                      <  2:2639645/150‑1 (MQ=255)
                                aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtg                                                                                                                     <  1:4355715/150‑1 (MQ=255)
                                aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgtg                                                                                                                     <  2:811260/150‑1 (MQ=25)
                                aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgtg                                                                                                                     <  1:230153/150‑1 (MQ=25)
                                 gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTTGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                    <  2:1299461/150‑1 (MQ=255)
                                 gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                    <  2:3934881/150‑1 (MQ=255)
                                 gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                    <  1:2544936/150‑1 (MQ=255)
                                 gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                    <  1:2073399/150‑1 (MQ=255)
                                 gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTgg                                                                                                                    <  1:2063179/150‑1 (MQ=25)
                                   cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGa                                                                                                                  <  1:4069183/150‑1 (MQ=255)
                                   cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCAGTGCCGTACCAGTGGGa                                                                                                                  >  2:1469558/1‑150 (MQ=255)
                                    aGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAg                                                                                                                 <  1:2265306/150‑1 (MQ=255)
                                       gTTACCGCGTTGAATTGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCgg                                                                                                              <  2:3814256/150‑1 (MQ=255)
                                       gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCgg                                                                                                              <  1:4120558/150‑1 (MQ=255)
                                       gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCgg                                                                                                              <  1:1751221/150‑1 (MQ=255)
                                       gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACACTGTGCCGTACCTGTGGGAGCgg                                                                                                              <  1:2130086/150‑1 (MQ=255)
                                       gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACACTGTGCCGTACCTGTGGGAGCgg                                                                                                              <  1:2130039/150‑1 (MQ=255)
                                         tcccGCGTTGAATCGCTGGCCCGGGTTAAGTCGCTGGAAGGTGGTGTGGTCTGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCAACAGGTGCAGCACAGGTATCGAACGCTGTGCCGTACCTGTGGGAGCGGGc                                                                                                            <  1:2010870/148‑1 (MQ=255)
                                             gcgTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                          <  1:637134/132‑1 (MQ=255)
                                             gcgTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                          >  2:637005/1‑132 (MQ=255)
                                             gcgTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                          <  1:637005/132‑1 (MQ=255)
                                                       ccTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCACCTTCATCGAATAAACCCGACACCTGCTTCGCGGGCACCCCCGCACCCACAGGTGCACCACAGCACTCCAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATgg                                                                                              <  2:2835930/149‑1 (MQ=38)
                                                         tGGCCGGGGTGACGTCGCTGGATGGTGGTGTGGTCAGCTTAATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCCCCGGTGCAGCACAGGTCTCTAACGCTGTGCCGAACGTGTGGGAGCGGGCGTGCCCGCGAATGGcc                                                                                            <  1:1837130/150‑1 (MQ=37)
                                                             ccGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCCCTCCCAAACGGTCAGCCCACGTCTAGAACGCGGTGCCGTCCCGGTGGGAGCGGGCGTGCCCGCGCATGGCCGcca                                                                                        >  2:1572178/1‑148 (MQ=255)
                                                                ggTGAAGTCGGGTGAAGGTGGTGTGGTAAGATTCATGGAATAAACGCGACACCTGCTTCGCGGGAACGCCCGCTCCCACCGTTGCCGCACTGGTCTCGACCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAAccc                                                                                     <  2:1215922/150‑1 (MQ=37)
                                                                    aaGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCtttt                                                                                 <  2:1850068/150‑1 (MQ=255)
                                                                         gCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACCCCTGCTTCGGGGGCACGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTg                                                                            <  2:2692190/150‑1 (MQ=255)
                                                                         gCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTTGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTg                                                                            <  2:3495659/150‑1 (MQ=255)
                                                                          cTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACCCCCGCGCCCACAGGGGCAGCCCAGGTCTCGAAAGCCGTGCCGGACCGGGGGGAGCGGGCGGGCCCCCGCAAGGCCGCCACCCCTGTAtaggg                                                                           >  1:3735745/1‑146 (MQ=255)
                                                                           tGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCCCGCCCGCGACCCCACGTTCAGCACAGGTCTCGAACCCAGTGCCCTAACAGTTGGAGCCGGCGTGCCCCCCAAAGGCCGCAACCCGTTTaacaggg                                                                          >  2:1967067/1‑144 (MQ=255)
                                                                              aaGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTgc                                                                       <  2:3409142/150‑1 (MQ=255)
                                                                                        ggtCAGCTTCATCGAATAACCGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCattcgc                                                             <  1:4511102/150‑6 (MQ=255)
                                                                                                  aTCGAATAAACGCGACACCTGCTTCGCGGGCACCCCCGCCCCCACAGGGGCAGCCCAGGTCGCCCACCCAGTGCCGGACCCGTGGGGGCGGGCGGGCCCCCCAAAGGCCCCAACCCCTTTAACAGGTTTCCCCCCCCCCCCCccacaaca                                                   >  2:613757/1‑150 (MQ=255)
                                                                                                   tCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCGCCCACACGGGCAGCAAAGGTGTCGAACGCGGTGGCGTAACAGTGGGAGCGGGCGTGCCCGCGGAAGGCCGCAAACCATTTTTATGGTTGGCGCCAGCCCCCCcaacaacaa                                                  >  2:619466/1‑150 (MQ=25)
                                                                                                   tCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCCCCCCCCCAGGTGCAGCACAGGTGAGGAACGCAGGGGCGGACCAGTGGGGGCCGGCGTGCCCCCCAAATGCCGCCACCCATTTATCTGGTTGGCCCCACCCCCCCccacaacaa                                                  >  2:619475/1‑150 (MQ=255)
                                                                                                    cGAATAACCGCGACACCGGCTTCGCGGGCACGCCCGCTCCCACCGGTGCAGCACTGGTCTCGAACGCGGTGCCGGACCTGTGGGCGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAg                                                 <  1:1469558/150‑1 (MQ=38)
                                                                                                     gAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGGGCCGCACCGGTCTCGAACGCTGTGCCGGACCGGTGGGCGCGGGCGGGCCCCCGAAAGGCCGCAACCCGTTTATCTGGTTGCCGCCACCCCCGCCAACAACCAgc                                                >  2:307989/1‑150 (MQ=37)
                                                                                                       aTAAACGCCAAACCTGCTTCGCGGGAACCCCCGCTCCCACAGGTGCACCACAGGTCTCGAACGATTTGCCGTACCTGTGGCAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAgcgc                                              <  1:2981865/150‑1 (MQ=255)
                                                                                                                   ccTGCTTCGCGGGCACGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttct                                  <  2:4171482/150‑1 (MQ=255)
                                                                                                                    atGCTTCGCGGGCACGCCCGCTCCAACAGGTTCAGCACAGGTCTAGAACCCTGTGCCGAACCAGTGGGAGCGTGCGTGCCCGCGAATGGCAGCAAACCTTTTATCTGGTTGCCCCCATCCCCGCCAAAAACAAGAGCTTGTAAAgatatt                                 >  2:1462705/2‑145 (MQ=21)
                                                                                                                      gcgtCGGGGGCACGCCCGCTCCCACTGGTGCCGCACAGGTCCCCCACGCTGTGCCGCACCTGGGGGCGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTtcttct                               <  1:2781048/147‑1 (MQ=255)
                                                                                                                                   gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCg                  >  2:1111400/1‑150 (MQ=255)
                                                                                                                                     ccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGc                >  1:1532717/1‑150 (MQ=255)
                                                                                                                                     ccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGc                >  1:2406864/1‑150 (MQ=255)
                                                                                                                                     ccGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGc                >  2:3057628/1‑150 (MQ=21)
                                                                                                                                      cGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCt               >  1:2760103/1‑150 (MQ=255)
                                                                                                                                      cGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCt               >  2:336991/1‑150 (MQ=255)
                                                                                                                                      cGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCt               >  1:2182024/1‑150 (MQ=255)
                                                                                                                                      cGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCt               >  2:336975/1‑150 (MQ=255)
                                                                                                                                         ccccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGCTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc            <  1:968236/149‑1 (MQ=255)
                                                                                                                                           ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc          >  2:848625/1‑150 (MQ=255)
                                                                                                                                           ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc          >  1:1184201/1‑150 (MQ=255)
                                                                                                                                           ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc          >  2:2789250/1‑150 (MQ=255)
                                                                                                                                           ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCCCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc          >  2:3903441/1‑150 (MQ=255)
                                                                                                                                           ccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc          >  1:1034907/1‑150 (MQ=25)
                                                                                                                                           ccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc          >  2:1408592/1‑150 (MQ=25)
                                                                                                                                           ccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc          >  2:1381194/1‑150 (MQ=25)
                                                                                                                                            cacaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCa         >  1:1610073/1‑150 (MQ=255)
                                                                                                                                            cacaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCa         >  1:1610125/1‑150 (MQ=255)
                                                                                                                                             acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa        >  2:179880/1‑150 (MQ=255)
                                                                                                                                             acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa        >  2:1092660/1‑150 (MQ=255)
                                                                                                                                             acaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa        >  1:2753235/1‑150 (MQ=25)
                                                                                                                                              caGGTGCTGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAAc       <  2:3988903/150‑1 (MQ=255)
                                                                                                                                               aGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACt      >  1:4373083/1‑150 (MQ=255)
                                                                                                                                               aGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACt      >  1:616772/1‑150 (MQ=255)
                                                                                                                                                 gTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAATTGc    >  1:2744245/1‑150 (MQ=255)
                                                                                                                                                 gTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGc    >  1:76377/1‑150 (MQ=255)
                                                                                                                                                 gTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGCGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGc    >  2:2754952/1‑150 (MQ=255)
                                                                                                                                                  tGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCa   >  1:1754896/1‑150 (MQ=255)
                                                                                                                                                   gcagcaCAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCAt  >  1:617906/1‑150 (MQ=255)
                                                                                                                                                   gcagcaCAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCAt  >  2:3398079/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCAT  >  NC_002947/6015667‑6015963

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: