Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,015,825 C→T 19.3% intergenic (+52/+58) 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,8250CT19.3% 202.6 / 39.0 99intergenic (+52/+58)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base C (27/52);  new base T (7/11);  total (36/63)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 105 total aligned reads displayed.

AGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGG  >  NC_002947/6015676‑6015973
                                                                                                                                                     |                                                                                                                                                    
aGGGCGGCCGTGCCGAACTGGACAGTCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtc                                                                                                                                                      >  1:1789138/1‑150 (MQ=255)
  ggcggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtctc                                                                                                                                                    <  2:907162/150‑1 (MQ=255)
  ggcggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtctc                                                                                                                                                    <  2:1866910/150‑1 (MQ=255)
  ggcggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGtctc                                                                                                                                                    <  2:1866610/150‑1 (MQ=255)
     ggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGaa                                                                                                                                                 <  2:1102301/150‑1 (MQ=255)
     ggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGaa                                                                                                                                                 <  2:1102238/150‑1 (MQ=255)
     ggcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGaa                                                                                                                                                 <  1:4093001/150‑1 (MQ=255)
      gcCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAAc                                                                                                                                                <  2:1326144/150‑1 (MQ=255)
         gTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCt                                                                                                                                             <  1:1028266/150‑1 (MQ=255)
            ccGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCtgtg                                                                                                                                          <  1:4960258/150‑1 (MQ=255)
             cGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                                         <  1:1323612/150‑1 (MQ=255)
             cGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGc                                                                                                                                         <  1:1323813/150‑1 (MQ=255)
              gAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCCCAGGGGGAGCCCAGGGGTCGAAACCGGTGcc                                                                                                                                        >  1:5316000/1‑150 (MQ=255)
              gAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGGGCAGCACAGGTCTCGAACGCTGTGcc                                                                                                                                        >  2:2695986/1‑150 (MQ=255)
               aaCTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCACCCCCACGGGCAGCACAGGGGTCCAAACCTGTGccc                                                                                                                                       >  1:3976630/1‑149 (MQ=255)
                 cTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTa                                                                                                                                     <  1:3098841/150‑1 (MQ=255)
                   ggACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGATCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTAcc                                                                                                                                   <  2:3463496/150‑1 (MQ=255)
                     aCAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCCCGCCCCCACGTGCACCACAGGTCTCGAACGCAGTGCCGTAccgg                                                                                                                                 >  2:2835160/1‑148 (MQ=255)
                      cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGTGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                                <  2:3263051/150‑1 (MQ=25)
                      cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                                <  1:3642335/150‑1 (MQ=255)
                      cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                                <  1:3642153/150‑1 (MQ=255)
                      cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                                <  1:838364/150‑1 (MQ=255)
                      cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgt                                                                                                                                <  1:4486983/150‑1 (MQ=25)
                      cAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCACTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgt                                                                                                                                <  1:564402/150‑1 (MQ=255)
                       aGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCtgtg                                                                                                                               <  2:321281/150‑1 (MQ=25)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTTGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                              >  2:2409126/1‑150 (MQ=255)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTTGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                              <  1:2409126/150‑1 (MQ=255)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGGAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTgg                                                                                                                              <  2:2907960/150‑1 (MQ=21)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtcg                                                                                                                              <  1:2938683/150‑3 (MQ=255)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                              <  1:3115730/150‑1 (MQ=255)
                        gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTgg                                                                                                                              <  2:3592536/150‑1 (MQ=25)
                         ccAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTggg                                                                                                                             <  1:4797730/150‑1 (MQ=21)
                          cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGTGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGa                                                                                                                            <  2:1140045/150‑1 (MQ=21)
                          cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCCGCACAGGTCTCGAACCCGGTGCCCTACCAGTgggg                                                                                                                            >  2:2778571/1‑149 (MQ=255)
                          cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGa                                                                                                                            <  2:1140059/150‑1 (MQ=21)
                          cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGa                                                                                                                            <  2:3880504/150‑1 (MQ=21)
                              gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                              <  1:1074152/144‑1 (MQ=255)
                              gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTgg                                                                                                                              >  2:1074152/1‑144 (MQ=255)
                              gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCgg                                                                                                                        <  2:3279997/150‑1 (MQ=21)
                               ttACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCggg                                                                                                                       <  2:4851399/150‑1 (MQ=255)
                                tACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGc                                                                                                                      <  1:1019854/150‑1 (MQ=255)
                                  ccGCGTTGACTCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCTAATAACCGCGCCACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGt                                                                                                                    <  2:1720511/150‑1 (MQ=255)
                                     cgTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGc                                                                                                                      <  2:3333809/145‑1 (MQ=18)
                                     cgTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGc                                                                                                                      >  1:3333809/1‑145 (MQ=18)
                                     cgTTGAAGCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGATTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGcc                                                                                                                 <  1:4465265/150‑1 (MQ=255)
                                                    ccGGGTGAAGTCGCTGGAAGGTGGTGTGTTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCTAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCaa                                                                                                  <  2:2637498/150‑1 (MQ=255)
                                                    ccGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAACCGCGCCCCCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCaa                                                                                                  <  2:2637289/150‑1 (MQ=255)
                                                    ccGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCaa                                                                                                  >  2:5374921/1‑150 (MQ=255)
                                                     cGGGTGAAGTCGCTGGAAGGTGGGGTGGTCAGCTTCATCGAATAAACGCGACGCCGGCTTCGCGGGCACGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAAc                                                                                                 <  2:3238723/150‑1 (MQ=255)
                                                          gAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTTTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCttt                                                                                            <  2:2203099/150‑1 (MQ=255)
                                                           ttgTCGGTGGATGGTGGTGTGGTCAGCTTCCGCGAATCACCGCGCCCCCGGCTTCGCGGGCCCGCGCGCTCCCCCCGGTGCCGCACCGGTCTCGAACGCTGTGCGGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCtttt                                                                                           <  2:5095876/148‑1 (MQ=25)
                                                           atgtggCTGGATGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACCCGTGCTTCGCGGGCCCGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGGGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCtttt                                                                                           <  2:2366520/145‑1 (MQ=38)
                                                                gCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGTGAATGGCCGCAACCCTTTTATCTg                                                                                      >  2:109121/1‑150 (MQ=255)
                                                                           gtgtggtCAGCTTCATCGAATACCCGCGCCCCCGGCTTCGCGGGCACGCCCGCTCCCACTGGTGCCGCACAGGTCTCGACCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAg                                                                           <  2:296744/150‑1 (MQ=37)
                                                                                 tCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACTGGTGCAGCACAGGTCTCGAACGCTGTGCCGAACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGc                                                                     <  1:3130861/150‑1 (MQ=255)
                                                                                    gCTTCACCGAAAAACCGCGACACCTGCTTCGCGGGCCCGCCCGCTCCCACAGGTGCCGCACCGGTCCCGACCGCTGTGCCGTACGTGGGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaa                                                                  <  1:397417/150‑1 (MQ=38)
                                                                                        aaTCTAATAAACGCGACCCCTGCTTGGGGGGCCCGCCCGCTCCCACTGTTGCAGCACTGGTCTCGAACGGTGGGCCGTACCTGGGGGAGGGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCcaacaac                                                              <  1:4647230/149‑1 (MQ=255)
                                                                                          tCGAATAAACGCGACACCTGCTTCGCGGGCCCGCCCCCGCCCCCCCGGGCGCCCCCCCGCGCGGACGCGGTGCCGGACCAGAGGGTGCGGGCGGGCCCGCCCATTGCCCCACACCGTTTTTCTGGGTGGCGCCCGCCCCCCcaacaacaa                                                            >  2:3670283/1‑150 (MQ=255)
                                                                                               tAACCGCTACACCTGCTTCGGGGGCACGCCCGCTCCCACTGGTGAAGCACAGGGCTCGAACGGTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCt                                                       <  1:4809475/150‑1 (MQ=38)
                                                                                                aaaCGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTGCAGCACTGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                      <  2:182879/150‑1 (MQ=255)
                                                                                                     cgACACCTGGTGCGCGGCCACCCGCGCACCAACTGTTGCAGCACTGTTCTCCACGGCTGTGCGTTCCCTGTGGGCGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTaaa                                                 <  1:2748856/150‑1 (MQ=255)
                                                                                                         acCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtg                                                                                                                               >  2:2752865/1‑68 (MQ=255)
                                                                                                         acCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCtgtg                                                                                                                               <  1:2752865/68‑1 (MQ=255)
                                                                                                          ccTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCTTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttct                                            <  2:5399755/150‑1 (MQ=255)
                                                                                                                 cgcgGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTaa                                     <  1:1795773/150‑1 (MQ=255)
                                                                                                                     ggCACGCCCGCTCCCACCGGTGCCGCACCGTTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGccc                                 <  1:4050280/150‑1 (MQ=255)
                                                                                                                      gCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGcccc                                <  1:5374921/150‑1 (MQ=255)
                                                                                                                        aCGCCCGCTCCCACAGGTACGGCACAGCGTTCGAGACCTGTGCAGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCAGGTTGCCGCCACCCCCGCCAACAACAAGAGCTGGTAAAGTTCATCTTTAAGccccag                              >  1:435954/1‑148 (MQ=18)
                                                                                                                           cccGCTCCCACCGGTGCAGCACAGGTCTCGACCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCgg                           <  2:983782/150‑1 (MQ=255)
                                                                                                                           cccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTgcc                                                                                >  1:298187/1‑97 (MQ=255)
                                                                                                                           cccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTgcc                                                                                <  2:298187/97‑1 (MQ=255)
                                                                                                                            ccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGc                          >  2:5108937/1‑150 (MQ=255)
                                                                                                                                 cccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                     >  1:1262981/1‑150 (MQ=255)
                                                                                                                                  ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCaa                                                                                                  <  1:1709640/72‑1 (MQ=255)
                                                                                                                                  ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCaa                                                                                                  >  2:1709640/1‑72 (MQ=255)
                                                                                                                                  ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                    >  1:1933468/1‑150 (MQ=255)
                                                                                                                                  ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                    >  2:2237075/1‑150 (MQ=255)
                                                                                                                                  ccACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                    >  1:955789/1‑150 (MQ=255)
                                                                                                                                  ccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                    >  1:1105688/1‑150 (MQ=25)
                                                                                                                                  ccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                    >  1:1560357/1‑150 (MQ=25)
                                                                                                                                  ccACAGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                    >  1:1105666/1‑150 (MQ=25)
                                                                                                                                   cacaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGTTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCa                   >  2:1555264/1‑150 (MQ=21)
                                                                                                                                    acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                  >  2:273569/1‑150 (MQ=255)
                                                                                                                                    acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                  >  2:979327/1‑150 (MQ=255)
                                                                                                                                    acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCCGCGGCTGCCCCaa                  >  1:4070885/1‑150 (MQ=255)
                                                                                                                                    acaGGTACGGCACAGCGTTCGAGACCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTACAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCaa                  >  2:608593/1‑150 (MQ=25)
                                                                                                                                     caGTTGCAGCACAGGTCTCGAACGGTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAAc                 <  2:2909957/150‑1 (MQ=255)
                                                                                                                                      aGGTGCAGCACAGGTCTCTAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACt                <  2:1454939/150‑1 (MQ=255)
                                                                                                                                      aGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCGGTGGGCGCGGGCGTGCCCGCGAATGGCCGCAACCCTTGTATCTGGTTGCCGCCAGCCCCGCCCACAACCAGCGCTGGTACAGTTCTTCGTTAAGCCCCTGCGGCTGCCCCAAcg                >  2:2923478/1‑149 (MQ=255)
                                                                                                                                        gTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGc              >  1:140946/1‑150 (MQ=255)
                                                                                                                                         tGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCa             >  2:434837/1‑150 (MQ=255)
                                                                                                                                              cacaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACCAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc                      >  1:1272520/1‑136 (MQ=255)
                                                                                                                                              cacaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc                      <  2:1272520/136‑1 (MQ=255)
                                                                                                                                              cacaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc                      >  1:1821544/1‑136 (MQ=255)
                                                                                                                                              cacaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcc                      <  2:1821544/136‑1 (MQ=255)
                                                                                                                                               acaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTAGTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGcc       >  2:3440928/1‑150 (MQ=255)
                                                                                                                                                   gTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAg   >  1:981954/1‑150 (MQ=255)
                                                                                                                                                   gTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAg   >  1:3773432/1‑150 (MQ=255)
                                                                                                                                                    tctcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAgg  >  2:1368771/1‑150 (MQ=255)
                                                                                                                                                    tctcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAgg  >  1:1705443/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGG  >  NC_002947/6015676‑6015973

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: