Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,031,616 G→A 5.6% G84S (GGT→AGT)  uacT → uric acid permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,031,6160GA 5.6% 229.1 / 2.6 71G84S (GGT→AGT) uacTuric acid permease
Reads supporting (aligned to +/- strand):  ref base G (35/32);  new base A (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

TCTGCAACACGTCCTTGTCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATCGCCGCGGGTTTTATCAC  >  NC_000913/3031468‑3031761
                                                                                                                                                    |                                                                                                                                                 
tCTGCAACACGTCCTTGTCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATgtatcg                                                                                                                                                   <  7:173899/149‑1 (MQ=255)
tCTGCAACACGTCCTTGTCATGTAAGCAGGTGCAGTCTCTGTTCATCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATgtatcg                                                                                                                                                   <  3:9837/149‑1 (MQ=255)
  tGCAACACGTCCTTGTCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATgtatcgg                                                                                                                                                  >  2:46056/1‑148 (MQ=255)
    cAACACGTCCTTGTCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATgtatcggtat                                                                                                                                               >  7:282582/1‑149 (MQ=255)
    cAACACGTCCTTGTCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATgtatcg                                                                                                                                                   <  1:210531/145‑1 (MQ=255)
    cAACACGTCCTTGTCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATgtatcg                                                                                                                                                   >  2:210531/1‑145 (MQ=255)
        acGTCCTTGTCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGc                                                                                                                                           >  1:119288/1‑149 (MQ=255)
          gTCCTTGTCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATgtatcggtatcg                                                                                                                                             >  3:281323/1‑145 (MQ=255)
          gTCCTTGTCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATgtatcggtatcg                                                                                                                                             <  4:281323/145‑1 (MQ=255)
          gTCCTTGTCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCg                                                                                                                                         >  6:163608/1‑149 (MQ=255)
          gTCCTTGTCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCg                                                                                                                                         >  4:55967/1‑149 (MQ=255)
              ttGTCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCtt                                                                                                                                      >  7:264346/1‑148 (MQ=255)
                 tCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATg                                                                                                                                  <  2:119288/149‑1 (MQ=255)
                   aTGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATggg                                                                                                                                <  4:2999/149‑1 (MQ=255)
                         gCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCg                                                                                                                          >  7:246310/1‑149 (MQ=255)
                           aGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGcc                                                                                                                        >  5:216194/1‑149 (MQ=255)
                            ggTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCt                                                                                                                       <  8:264346/149‑1 (MQ=255)
                              tGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGc                                                                                                                     >  8:314883/1‑149 (MQ=255)
                                 aGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGACGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCgg                                                                                                                  <  3:988/149‑1 (MQ=255)
                                  gTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGt                                                                                                                 <  2:212629/149‑1 (MQ=255)
                                   tCGCTGTTCCTCGTATGAGTGGTGACCGCCTGGGCCTCTCAAAAGAAGCGATTGCGATGATCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTGACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCAGGt                                                                                                                 >  6:143725/1‑148 (MQ=255)
                                        gTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTAt                                                                                                           >  3:330909/1‑149 (MQ=255)
                                            ctctTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTc                                                                                                        <  2:178864/148‑1 (MQ=255)
                                             tctTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCgg                                                                                                      >  7:363544/1‑149 (MQ=255)
                                                tATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCAGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGa                                                                                                   >  5:256026/1‑149 (MQ=255)
                                                tATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCAGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGa                                                                                                   >  7:154624/1‑149 (MQ=255)
                                                 aTGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGt                                                                                                     <  8:282582/146‑1 (MQ=255)
                                                 aTGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGt                                                                                                     <  4:330909/146‑1 (MQ=255)
                                                  tGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGAcc                                                                                                 <  8:363544/149‑1 (MQ=255)
                                                    aTTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCGCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCt                                                                                                <  4:20610/148‑1 (MQ=255)
                                                                   ggCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGTATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAaca                                                                                  <  2:340722/147‑1 (MQ=255)
                                                                   ggCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAaca                                                                                  <  1:46056/147‑1 (MQ=255)
                                                                     ccTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCaa                                                                              >  7:246252/1‑149 (MQ=255)
                                                                      ctctCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAAt                                                                             <  7:161522/149‑1 (MQ=255)
                                                                         tcAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGa                                                                           >  6:197408/1‑148 (MQ=255)
                                                                           aaaaGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATgtatcggtatc                                                                                                                                              >  4:339992/1‑79 (MQ=255)
                                                                           aaaaGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATTGGGATCCGCCTGCCGGTGATTATGTCGGt                                                                                                     <  4:51602/120‑1 (MQ=255)
                                                                           aaaaGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATTGGGATCCGCCTGCCGGTGATTATGTCGGt                                                                                                     >  3:51602/1‑120 (MQ=255)
                                                                           aaaaGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGATCTTTGCTGCTGTAACACCAATGATAg                                                                        >  3:93998/1‑149 (MQ=255)
                                                                           aaaaGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAg                                                                        >  4:157708/1‑149 (MQ=255)
                                                                           aaaaGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCAGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATa                                                                         <  8:154624/148‑1 (MQ=255)
                                                                           aaaaGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCAGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATa                                                                         <  6:256026/148‑1 (MQ=255)
                                                                            aaaGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTc                                                                                                        <  5:351399/116‑1 (MQ=255)
                                                                            aaaGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTc                                                                                                        >  6:351399/1‑116 (MQ=255)
                                                                             aagaagCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATACGCCTGCCGGTGATTATGTCGGAGACCTTTGCTGCTGTAACACCAATGATAGc                                                                       >  6:78002/1‑148 (MQ=255)
                                                                                 agCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTg                                                                  <  7:314883/149‑1 (MQ=255)
                                                                                  gCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTc                                                                                                        >  6:96754/1‑110 (MQ=255)
                                                                                  gCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTc                                                                                                        <  5:96754/110‑1 (MQ=255)
                                                                                    tATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGa                                                               >  4:250485/1‑149 (MQ=255)
                                                                                     aTTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGAt                                                              >  3:7329/1‑149 (MQ=255)
                                                                                      ttccGATGCTCATTAGCTAGGATCTCTTTTGCTGCGGGATCGGCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGa                                                               <  6:216194/144‑1 (MQ=255)
                                                                                        gCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGAt                                                                                                              <  8:256285/98‑1 (MQ=255)
                                                                                        gCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGAt                                                                                                              >  7:256285/1‑98 (MQ=255)
                                                                                         cGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAAc                                                          >  1:328828/1‑149 (MQ=255)
                                                                                         cGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAAc                                                          <  3:265936/149‑1 (MQ=255)
                                                                                          gATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAAcc                                                         >  5:273088/1‑149 (MQ=255)
                                                                                              cTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCgg                                                      <  1:382550/148‑1 (MQ=255)
                                                                                                cATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGat                                                    >  4:302898/1‑148 (MQ=255)
                                                                                                    aGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGATCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCg                                                <  4:93998/148‑1 (MQ=255)
                                                                                                      cTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGcc                                             >  6:141329/1‑149 (MQ=255)
                                                                                                               ctTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTgggg                                     <  7:15010/148‑1 (MQ=255)
                                                                                                               ctTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTgggg                                     >  8:15010/1‑148 (MQ=255)
                                                                                                                ttttGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTgg                                       <  3:55967/145‑1 (MQ=255)
                                                                                                                   tgctgcGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGatat                                 <  6:273088/148‑1 (MQ=255)
                                                                                                                        cGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTgg                             <  5:141329/147‑1 (MQ=255)
                                                                                                                        cGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGt                            <  4:7329/148‑1 (MQ=255)
                                                                                                                        cGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTg                           <  5:197408/149‑1 (MQ=255)
                                                                                                                         gggATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGc                          <  8:246310/149‑1 (MQ=255)
                                                                                                                          ggATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGc                          >  7:275756/1‑148 (MQ=255)
                                                                                                                                tCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATc                   >  7:57806/1‑149 (MQ=255)
                                                                                                                                 cacaTTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATcg                  >  3:295199/1‑149 (MQ=255)
                                                                                                                                    attattGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGAGGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATcgccg               >  2:251195/1‑149 (MQ=255)
                                                                                                                                    attattGCAATGTATCGGTATCGGCCGCTTTATGAGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGACCCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATcgccg               >  5:25005/1‑149 (MQ=255)
                                                                                                                                        ttGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATCGCCGCgg            <  3:157708/148‑1 (MQ=255)
                                                                                                                                          gCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATCGCCGCGGGtt         >  5:282775/1‑149 (MQ=255)
                                                                                                                                          gCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATCGCCGCGGGt          <  8:275756/148‑1 (MQ=255)
                                                                                                                                           cAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATCGCCGCGGGttt        >  7:352003/1‑149 (MQ=255)
                                                                                                                                            aaTGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGAGAGCGGCCAGCTGGGGAGATTTGGTGCCACTATCGCCGCGGGgtt        >  8:403032/1‑148 (MQ=255)
                                                                                                                                            aaTGTATCGGTATCGGCCGCGTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATCGCCGCGGGtttt       <  7:403032/149‑1 (MQ=255)
                                                                                                                                              tGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATCGCCGCGGGTTTTAt     >  8:410480/1‑149 (MQ=255)
                                                                                                                                                tatcggtatcggCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATCGCCGCGGGTTTTATc    <  3:250485/148‑1 (MQ=255)
                                                                                                                                                  tcggtatcggCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGGCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGACGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATCGCCGCGGGTTTTATcac  <  1:251195/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCTGCAACACGTCCTTGTCATGTACGCAGGTGCAGTCGCTGTTCCTCTTATGATTGGTGACCGACTGGGCCTCTCAAAAGAAGCTATTGCGATGCTCATTAGCTCGGATCTCTTTTGCTGCGGGATCGTCACATTATTGCAATGTATCGGTATCGGCCGCTTTATGGGGATCCGCCTGCCGGTGATTATGTCGGTGACCTTTGCTGCTGTAACACCAATGATAGCCATTGGGATGAACCCGGATATCGGCCTGCTGGGGATATTTGGTGCCACTATCGCCGCGGGTTTTATCAC  >  NC_000913/3031468‑3031761

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: