Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,018,868 T→C 16.6% intergenic (+68/+42) 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,8680TC16.6% 361.7 / 58.6 158intergenic (+68/+42)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (44/81);  new base C (12/14);  total (63/95)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.09e-01

Only 100 of 161 total aligned reads displayed.

CTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTC  >  NC_002947/6018720‑6019017
                                                                                                                                                    |                                                                                                                                                      
cTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTa                                                                                                                                                       <  2:4554125/150‑1 (MQ=255)
cTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTa                                                                                                                                                       <  1:468647/150‑1 (MQ=255)
cTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTa                                                                                                                                                       <  1:200313/150‑1 (MQ=255)
cTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTa                                                                                                                                                       <  1:109518/150‑1 (MQ=255)
cTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTa                                                                                                                                                       <  2:6295085/150‑1 (MQ=255)
 tGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                                                      <  1:2715097/150‑1 (MQ=255)
 tGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                                                      <  1:942294/150‑1 (MQ=255)
 tGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAc                                                                                                                                                      >  1:5025616/1‑150 (MQ=255)
  ggACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCCGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAcc                                                                                                                                                     <  2:242213/150‑1 (MQ=255)
   gACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCt                                                                                                                                                    <  1:5347909/150‑1 (MQ=255)
    aCAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtg                                                                                                                                                   <  2:146912/150‑1 (MQ=255)
     cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                                                  <  2:495369/150‑1 (MQ=255)
     cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                                                  <  1:741871/150‑1 (MQ=25)
     cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                                                  <  1:1822967/150‑1 (MQ=25)
      aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtg                                                                                                                                                 <  2:5878932/150‑1 (MQ=255)
      aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtg                                                                                                                                                 <  1:4257707/150‑1 (MQ=255)
       gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                                                <  2:3928538/150‑1 (MQ=255)
       gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                                                <  1:3260858/150‑1 (MQ=255)
       gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                                                <  1:3234769/150‑1 (MQ=255)
       gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                                                <  2:5028483/150‑1 (MQ=255)
       gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                                                <  2:6564690/150‑1 (MQ=255)
       gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTgg                                                                                                                                                <  2:4559612/150‑1 (MQ=25)
       gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTgg                                                                                                                                                <  1:4696991/150‑1 (MQ=25)
       gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTgg                                                                                                                                                <  2:5022463/150‑1 (MQ=25)
       gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTgg                                                                                                                                                <  2:1899581/150‑1 (MQ=25)
        ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCCCGCCCACAGGTGCAGCACAGGTATCGAACGCGGTGCCCTACCAGTggg                                                                                                                                               >  2:5941695/1‑150 (MQ=255)
         cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACCCGGGCAGCACAAGGCGCGGACCCGGGGGCGGACCAGTGGGa                                                                                                                                              >  1:6193276/1‑150 (MQ=255)
            ggTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCg                                                                                                                                          <  1:899689/150‑1 (MQ=255)
             gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCgg                                                                                                                                         <  1:3556380/150‑1 (MQ=255)
                  cgcgTTGAATCGCTGGCCCGGGTGACGTCGCTGGAATGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTg                                                                                                                                    <  2:4821718/150‑1 (MQ=255)
                         aaTCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGGGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTACGGCACAGCGTTCGAGACCTGTGCGGCACCGGTGGGAGCGGGCGTGCCCGCGa                                                                                                                             <  2:6381179/150‑1 (MQ=18)
                          aTCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGCACGCTGTGCCGTACCAGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                            >  2:5533375/1‑150 (MQ=255)
                           tCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCCCGCCCACAGGTGCACCACAGGTCACGAACGCAGTGCCGTACCAGTGGGAGCGGGCGGGCCCCCGaaa                                                                                                                           >  2:3196541/1‑149 (MQ=255)
                               tGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTACAGCGAATAAACGCGACCCCTGCTTCGCGGGCACGCCCGCTCCCCCCGGTGCAGCACTGGTCTCGAAGGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGcc                                                                                                                       <  2:4847596/150‑1 (MQ=37)
                                         gAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAACCGCGACACCGGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCttt                                                                                                             <  1:5122153/150‑1 (MQ=255)
                                             tCGCTGGACGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATc                                                                                                         <  2:2053485/150‑1 (MQ=255)
                                             tCGCTGGAAGGGGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATc                                                                                                         <  2:5205363/150‑1 (MQ=255)
                                               gCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGGACCAGTGGGAGCGGGCGGGCCCCCGAATGGCCGCAACCCGTTTATCTg                                                                                                       >  2:1667310/1‑150 (MQ=38)
                                               gCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGGCTCGAACCCAGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGCCGCAACCCGTTTATATg                                                                                                       >  2:457186/1‑150 (MQ=37)
                                                gtGGATGGTGGTGGGGTCAGCTACATCGAATAATCGCGACACCGGCTTCGCGGGCACGCCCGCTCCCCCTGGTGCAGCACCGGTCCCTACCGGTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTgg                                                                                                      <  2:3069350/149‑1 (MQ=37)
                                                        tggtgtggtCAGCTTCATTGAATAAACGCGACACCTGCTTCGGGGGCACGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTgccgcc                                                                                              <  2:2374967/150‑1 (MQ=255)
                                                        gggggtggtCAGCTTGATTGACTAAACGCGACACCGGCTTCGCGGGCAGGCCCGCTCCAACAGGTGAAGCACTGGTCACGACCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCTAATGGCCGCAAACCTTTTATCTGGTTgccgcc                                                                                              <  2:5942547/146‑1 (MQ=37)
                                                                     ttCATCGAATAACCGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACTGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaaca                                                                                 <  1:6685538/150‑1 (MQ=255)
                                                                     ttCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTGCAGCACAGGTCACGAACGCTGTGCCGTCCCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaaca                                                                                 <  1:4546182/150‑1 (MQ=255)
                                                                     ttCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCCCTGCCACCGGTTCACCACACCGCTAGAACGCGGTGCCCGACCTGTGGGAGCGGGCCAGGCACCCAAAGGCCCCCACCCATTTATCTTGTTTCCGCCACCCCCCCccacc                                                                                 >  2:2552637/1‑149 (MQ=255)
                                                                               aatcgcgACACGTGGTTCGCGGGCCGGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCGGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                       <  1:3801593/147‑1 (MQ=38)
                                                                               aaccgcgACACCTGCTTGGCGGGCACGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                       <  1:3320212/147‑1 (MQ=255)
                                                                               aaaCGCGACACCTGCTTCGGGGGCACGCCCGCTCCCACTGGTGCAGCACATGTCTCTAAGGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                       <  1:2068920/150‑1 (MQ=38)
                                                                               aaaCGCGACACCTGCTTCGGGGGCACGCCCGCTCCCACTGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                       <  2:2567746/150‑1 (MQ=255)
                                                                               aaaCGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                       <  2:4155588/150‑1 (MQ=255)
                                                                               aaaCGCGACACCTGCTGCGCGGGCACGCCCGCTCCCACAGGTGCAGCACTGTTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                       <  2:746209/150‑1 (MQ=255)
                                                                                    cgACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTaaa                                                                  <  2:5609211/150‑1 (MQ=255)
                                                                                           tGCTTCGCGGGCACGCCCGCTCCCCCACGTGCAGCACAGGTCTCCCACGCGGTTCCGGACCTGTGGGAGCGGGCGGGCCCCCCAATGGCCGCAACCCGTTTATCTGGTTTCCCCCCCCCCCCCCACCAACCAGCGCGGGGAAAGttaatg                                                           >  2:786506/1‑146 (MQ=255)
                                                                                           tGCTTCGCGGGCACCCCCGCTCCCCCAGGGGCAGCACAGGTCTAGAACGCAGTGCCCTACCAGTGGGACCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCAGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAActtcttg                                                           >  2:4761551/1‑149 (MQ=255)
                                                                                                cgcgGGCACGCCCGCTCCCACAGGTGCAGCACCGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACCAGCGCTGGTAAAGTTCTTCTTTaa                                                      <  2:1032100/150‑1 (MQ=255)
                                                                                                   gggCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGt                                                                                                                                     >  2:1591046/1‑68 (MQ=255)
                                                                                                   gggCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCATTTATCTGGTTGCCGCCAGCCCCGCCCACAACAAGCGCTGGTAAAGTTCTTCTTTAAGcc                                                   >  1:4477205/1‑150 (MQ=255)
                                                                                                     gCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGGGCTGCCCCGGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCGTGTATCTGGTTGCCGCCAGCCCCGCCCACCACAAGCGCTGGTAAAGTTCGTCTGTAAGcccc                                                 >  2:1056387/1‑150 (MQ=17)
                                                                                                        cGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                              <  1:2827450/150‑1 (MQ=255)
                                                                                                         gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCg                                             >  2:5934131/1‑150 (MQ=255)
                                                                                                           ccGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGt                                                                     >  2:4953517/1‑124 (MQ=17)
                                                                                                           ccGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGt                                                                     <  1:4953241/124‑1 (MQ=17)
                                                                                                             gCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTg                                         >  2:2198079/1‑150 (MQ=21)
                                                                                                              cTCCCACCGGTGCAGCACAGGTCTCGAACGCTGGGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                                        <  1:1638564/150‑1 (MQ=255)
                                                                                                               tCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                                        >  2:452483/1‑150 (MQ=255)
                                                                                                               tCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGGCGTGCCCGCGAATGGCCGCAACCCTGTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                                        >  2:6215788/1‑150 (MQ=255)
                                                                                                               tCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCtt                                                                                                              >  2:2504126/1‑79 (MQ=255)
                                                                                                               tCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCtt                                                                                                              <  1:2504126/79‑1 (MQ=255)
                                                                                                               tCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc                                       >  2:504869/1‑150 (MQ=21)
                                                                                                                cccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                                      >  2:5775647/1‑150 (MQ=255)
                                                                                                                cccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                                      <  2:5550302/150‑1 (MQ=255)
                                                                                                                 ccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                                     >  1:5075976/1‑150 (MQ=25)
                                                                                                                   acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                                   >  2:3871459/1‑150 (MQ=255)
                                                                                                                   acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                                   >  1:6728261/1‑150 (MQ=255)
                                                                                                                   acaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                                   >  1:5503753/1‑150 (MQ=25)
                                                                                                                   acaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                                   >  1:1612793/1‑150 (MQ=25)
                                                                                                                     aGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACt                                 >  2:4004003/1‑150 (MQ=255)
                                                                                                                     aGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACt                                 <  2:3950462/150‑1 (MQ=255)
                                                                                                                       gTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGc                               >  1:4955334/1‑150 (MQ=255)
                                                                                                                       gTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGc                               >  2:6706682/1‑150 (MQ=255)
                                                                                                                             ctgaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGc                         <  1:4399097/147‑1 (MQ=255)
                                                                                                                             cacaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTAACTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGc                         >  2:5920148/1‑150 (MQ=255)
                                                                                                                                aGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGcccc                      >  1:1815272/1‑150 (MQ=255)
                                                                                                                                  gTCTCGAACGCTTGGCCGTACCTGTGGGCGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAg                    <  1:5990792/150‑1 (MQ=255)
                                                                                                                                  gTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAg                    >  2:6547703/1‑150 (MQ=255)
                                                                                                                                   tctcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAgg                   >  1:1489505/1‑150 (MQ=255)
                                                                                                                                     tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttc                                                              <  1:4339763/105‑1 (MQ=255)
                                                                                                                                     tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttc                                                              >  2:4339763/1‑105 (MQ=255)
                                                                                                                                      cgagaccTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTgg                >  2:979537/7‑150 (MQ=255)
                                                                                                                                      cGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTgg                <  2:408348/150‑1 (MQ=255)
                                                                                                                                      cGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTgg                >  1:6030699/1‑150 (MQ=255)
                                                                                                                                       gAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAgc                                                                           <  2:2785912/90‑1 (MQ=255)
                                                                                                                                       gAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAgc                                                                           >  1:2785912/1‑90 (MQ=255)
                                                                                                                                       gAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGc               >  2:2178614/1‑150 (MQ=255)
                                                                                                                                        aaCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCt              <  2:2421808/150‑1 (MQ=255)
                                                                                                                                                 ccGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATAccc     >  1:5510963/1‑150 (MQ=255)
                                                                                                                                                   gTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGTGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTc  >  1:3849697/1‑150 (MQ=255)
                                                                                                                                                    tACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                                        >  2:3899554/1‑112 (MQ=255)
                                                                                                                                                    tACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                                        <  1:3899554/112‑1 (MQ=255)
                                                                                                                                                    tACCTGGGGGCGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTc  <  1:355156/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTC  >  NC_002947/6018720‑6019017

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: