Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 513,502 C→T 5.1% S309S (AGC→AGT ybaT → putative amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913513,5020CT 5.1% 250.1 / 2.3 78S309S (AGC→AGTybaTputative amino acid transporter
Reads supporting (aligned to +/- strand):  ref base C (37/36);  new base T (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

GCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGTTGCCAGCGCCACCTTTTTGATT  >  NC_000913/513356‑513649
                                                                                                                                                  |                                                                                                                                                   
gCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAgc                                                                                                                                                     <  8:332471/147‑1 (MQ=255)
gCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAgc                                                                                                                                                     >  7:332471/1‑147 (MQ=255)
 ccGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAgcgaa                                                                                                                                                  >  8:142831/1‑149 (MQ=255)
 ccGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAgcgaa                                                                                                                                                  <  2:135178/149‑1 (MQ=255)
    aTACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCTAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAgcgaacgc                                                                                                                                               <  6:133367/149‑1 (MQ=255)
        cgccgTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAgcgaacgcgaac                                                                                                                                           >  3:377497/1‑149 (MQ=255)
          ccgTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTg                                                                                                                                         <  2:154913/149‑1 (MQ=255)
               gcgcAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGaa                                                                                                                                    <  3:197410/149‑1 (MQ=255)
               gcgcAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGa                                                                                                                                     <  6:330801/148‑1 (MQ=255)
                        gctTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGaa                                                                                                                           >  3:8065/1‑149 (MQ=255)
                         ctTCTCCGCTGCTCGGGCATGGGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGCCGGCTTCAGCCATTAACGCGAACCTGTGCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGa                                                                                                                            <  5:75303/147‑1 (MQ=255)
                          ttCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAAt                                                                                                                          >  4:257721/1‑148 (MQ=255)
                               cGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAAt                                                                                                                              >  7:5973/1‑139 (MQ=255)
                               cGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAAt                                                                                                                              <  8:5973/139‑1 (MQ=255)
                               cGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATcc                                                                                                                    <  7:261167/149‑1 (MQ=255)
                               cGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATcc                                                                                                                    <  7:142831/149‑1 (MQ=255)
                                gctgctCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATccc                                                                                                                   >  7:2837/1‑149 (MQ=255)
                                    ctCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCtgtg                                                                                                               <  3:257721/149‑1 (MQ=255)
                                     tCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGTGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTgg                                                                                                              >  1:328810/1‑149 (MQ=255)
                                       gggCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTggcg                                                                                                            <  7:37877/149‑1 (MQ=255)
                                         gCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGTGAACGCGAACTGCCGAAGCTAAt                                                                                                                              >  1:316099/1‑129 (MQ=255)
                                         gCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGTGAACGCGAACTGCCGAAGCTAAt                                                                                                                              <  2:316099/129‑1 (MQ=255)
                                            tgtgGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACCTGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCaga                                                                                                       <  1:36941/149‑1 (MQ=255)
                                            tgtgGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCa                                                                                                         >  2:8457/1‑147 (MQ=255)
                                                 gTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACTCGAAACTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCatcttacc                                                                                                  >  2:330761/1‑142 (MQ=255)
                                                  ttATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCt                                                                                                 >  1:191290/1‑149 (MQ=255)
                                                  gtATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTGCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCt                                                                                                 <  2:241164/148‑1 (MQ=255)
                                                    aTGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTGAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACTCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTAccgg                                                                                                <  4:15409/148‑3 (MQ=255)
                                                     tgtgATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTa                                                                                                    >  6:250416/1‑143 (MQ=255)
                                                     tgtgATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTa                                                                                                    >  6:244701/1‑143 (MQ=255)
                                                     tgtgATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTa                                                                                                    <  5:250416/143‑1 (MQ=255)
                                                     tgtgATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTa                                                                                                    <  5:244701/143‑1 (MQ=255)
                                                     tgagATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTggg                                                                                              >  6:245633/4‑149 (MQ=255)
                                                      gtgATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTggg                                                                                              <  2:191290/148‑1 (MQ=255)
                                                       tgATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGc                                                                                            <  4:377497/149‑1 (MQ=255)
                                                         aTCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCaa                                                                                          >  8:205605/1‑149 (MQ=255)
                                                            gtcgtcATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACat                                                                                       >  7:287039/1‑149 (MQ=255)
                                                                tcATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTgg                                                                                               <  2:338874/137‑1 (MQ=255)
                                                                tcATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTgg                                                                                               >  1:338874/1‑137 (MQ=255)
                                                                  aTCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGt                                                                                 >  2:194651/1‑149 (MQ=255)
                                                                  aTCGGCGCTTTACTCGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACCTGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACattatt                                                                                   <  6:107859/147‑1 (MQ=255)
                                                                      gcgcTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCgtg                                                                             <  1:19694/149‑1 (MQ=255)
                                                                        gcTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCgtggt                                                                           <  8:2837/149‑1 (MQ=255)
                                                                            tACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGTACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGtt                                                                        <  5:192372/148‑1 (MQ=255)
                                                                               tGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGAt                                                                     >  4:326807/1‑148 (MQ=255)
                                                                                 gCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGAtt                                                                    <  5:188970/147‑1 (MQ=255)
                                                                                  cGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAAt                                                                                                                              <  3:22585/88‑1 (MQ=255)
                                                                                  cGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAAt                                                                                                                              >  4:22585/1‑88 (MQ=255)
                                                                                   gACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGc                                                                >  7:212702/1‑149 (MQ=255)
                                                                                       gCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCtgat                                                            <  5:245633/149‑1 (MQ=255)
                                                                                           cAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGTGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGAc                                                         <  2:328810/148‑1 (MQ=255)
                                                                                            aGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGAcgg                                                       >  6:67387/1‑149 (MQ=255)
                                                                                             gCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCgtg                                                                             >  7:184593/1‑126 (MQ=255)
                                                                                             gCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCgtg                                                                             <  8:184593/126‑1 (MQ=255)
                                                                                             gCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCACTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCTTGGTGTTGATTATGCTGATGAcggc                                                      >  4:70678/1‑149 (MQ=255)
                                                                                                aTTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGAcggcgg                                                    >  6:147265/1‑148 (MQ=255)
                                                                                                aTTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGAcggcgg                                                    >  1:225965/1‑148 (MQ=255)
                                                                                                 ttAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCa                                                  >  6:8672/1‑149 (MQ=255)
                                                                                                    aCGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACt                                                >  6:234870/1‑148 (MQ=255)
                                                                                                    aCGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCACTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACt                                                >  6:73923/1‑148 (MQ=255)
                                                                                                       agAACCTGTTCGCCGTGTTTAACATCATGGACAACATTGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTCTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTGGATTATGCTGATGACGGCGGCACTGAAt                                            <  2:280247/148‑1 (MQ=255)
                                                                                                          aCCTGTCCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTa                                         <  2:100376/149‑1 (MQ=255)
                                                                                                           ccTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAg                                        <  5:192570/149‑1 (MQ=255)
                                                                                                            cTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAgg                                       <  3:326807/149‑1 (MQ=255)
                                                                                                              gTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCt                                     >  1:306361/1‑149 (MQ=255)
                                                                                                              gTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCt                                     >  7:55693/1‑149 (MQ=255)
                                                                                                                   ccGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTc                                <  4:375847/149‑1 (MQ=255)
                                                                                                                        tttAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGc                                                                >  6:240004/1‑112 (MQ=255)
                                                                                                                        tttAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGc                                                                <  5:240004/112‑1 (MQ=255)
                                                                                                                           aaCATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGt                        >  7:152292/1‑149 (MQ=255)
                                                                                                                            acttcatGGACAACATGTGCAGCGAACGCGAACTGCCGAAGCTAATGAATAACTCCCTGTGGCGGCAGAGTACCTTGGTCAACATTATTGTCGTGGGGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGtt                       <  1:330761/146‑1 (MQ=255)
                                                                                                                               tcatGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGCATTTAGGCTCACTCGCCAGCGTTGc                     >  1:50636/1‑148 (MQ=255)
                                                                                                                               tcatGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGTTGcc                    <  3:70678/149‑1 (MQ=255)
                                                                                                                                 atGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGTTGCCAg                  >  7:409979/1‑149 (MQ=255)
                                                                                                                                  tGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGTTGCCAgc                 <  2:306361/149‑1 (MQ=255)
                                                                                                                                   ggACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAACGTTGCCAgcg                >  8:427350/1‑149 (MQ=255)
                                                                                                                                     acaacaTGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGTTGCCAGCGcc              >  8:31503/1‑149 (MQ=255)
                                                                                                                                        acaTGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTa                                                                   >  1:207955/1‑93 (MQ=255)
                                                                                                                                        acaTGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTa                                                                   <  2:207955/93‑1 (MQ=255)
                                                                                                                                         caTGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGTTGCCAGCGCCACCt          >  3:187390/1‑149 (MQ=255)
                                                                                                                                            gggCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGTTGCCAGCGCCACCtttt       >  1:213305/1‑149 (MQ=255)
                                                                                                                                             gtcAGAGAAAGCGAACTGCCGAAGCTAATTAATAAATCCCTGTGGCGTCAGAGTACCAGTGGCAACATTATTGTCGTTGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGTTGCCAGTGCCACCtttt       <  4:303675/146‑1 (MQ=255)
                                                                                                                                             ggCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGTTGCCAGCGCCACCttttt      >  6:14077/1‑149 (MQ=255)
                                                                                                                                             ggCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGTTGCCAGCGCCACCttttt      <  8:409979/149‑1 (MQ=255)
                                                                                                                                                 gcgaacgcgaacTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTg                                                                       <  6:151401/80‑1 (MQ=255)
                                                                                                                                                 gcgaacgcgaacTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTg                                                                       >  5:151401/1‑80 (MQ=255)
                                                                                                                                                 gcgaacgcgaacTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGTTGCCAGCGCCACCTTTTTGAtt  <  8:112178/149‑1 (MQ=255)
                                                                                                                                                 gcgaacgcgaacTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGTTGCCAGCGCCACCTTTTTGAtt  <  4:8065/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGGCTCACTCGCCAGCGTTGCCAGCGCCACCTTTTTGATT  >  NC_000913/513356‑513649

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: