Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,282,627 C→T Q922* (CAG→TAG)  narG → nitrate reductase A subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,282,6270CT83.3% 100.0 / 0.1 49Q922* (CAG→TAG) narGnitrate reductase A subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (8/0);  new base T (13/27);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 5.39e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCAC  >  NC_000913/1282382‑1282870
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          
gCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGAAAAAGGATTGCCAGGAACACCCagaac                                                                                                                                                                                                                                                       >  1:498916/1‑249 (MQ=255)
        ggTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGGAAAACCGCGCCACACATTATGGTCGTAGAGCGCGAGTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGa                                                                                                                                                                                                                                               <  2:393255/251‑1 (MQ=255)
            aaGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGCAAAGGGATTTCCTGGAACACCTAAAACGAGATGGaaacg                                                                                                                                                                                                                                           >  1:291107/1‑247 (MQ=255)
                   cGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAATGCGAGTGCGACCTGATCCCAGTTAAAACCGCGCCACACATTATTGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGACGGATAAAATCGGTAAGGGCGGTAAAGGGATTGCCTGGAACACCTAGCGCGAGATGGATCTGTTGCGTa                                                                                                                                                                                                                                    <  2:585823/251‑1 (MQ=255)
                   cGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGGTTGGCCGGAACACCCAGAGCGAGGTGGATCTGGTGCgg                                                                                                                                                                                                                                     >  1:291401/1‑249 (MQ=255)
                    gACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCCGGAACCCCTAGAGCGGGATGGATatgcttccga                                                                                                                                                                                                                                    >  1:436703/1‑245 (MQ=255)
                                cTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGGAAAAGGATTGCCTTGAACACCCCGAGCGAGATGGATCTGCTTCGTCAGCTCAAATacac                                                                                                                                                                                                                       >  1:424038/1‑251 (MQ=255)
                                     tATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAAAGCGAGGTGGAACTGCTGCGTAAGCCCCACTACAc‑ccaa                                                                                                                                                                                                                  >  1:299014/1‑247 (MQ=255)
                                      aTCCAGCATGACGCTGCCGCTGAACGGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGc                                                                                                                                                                                                               <  2:118072/251‑1 (MQ=255)
                                            cATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGAAAAGGGGATGCCTGGAACACCTAAAGCGGGGGGGATCTGCTGAGACAGCTCAAATAAAC‑CA‑AAGCGGCA‑gg                                                                                                                                                                                                        >  1:497864/1‑251 (MQ=255)
                                                                   gcAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAACGGGGTTTCCTGGGACACCTCGAGCGAGAAGGGGCGGTTGCGGAAGCTCAACTACAC‑CA‑AAGCGGGA‑GGAACGGCGCAAGG‑CCAGCCGATgg                                                                                                                                                                                >  2:19175/1‑250 (MQ=255)
                                                                       ccGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGaa                                                                                                                                                                           <  1:469633/251‑1 (MQ=255)
                                                                               tgtgAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGa                                                                                                                                                                                                            <  2:32922/213‑1 (MQ=255)
                                                                               tgtgAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGGTGGATCTGGTGCGGAAGCTCAACTACAC‑CA‑AAGCggg                                                                                                                                                                                                            >  1:32922/1‑212 (MQ=255)
                                                                                          ggAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCGTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGc                                                                                                                                                         <  1:578423/250‑1 (MQ=255)
                                                                                           gAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCGCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGTA‑GTTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCag                                                                                                                                                       <  1:226942/251‑1 (MQ=255)
                                                                                               aaaGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGATTGCCTGGGAAACCTAAAACGGGATGGATCTGCTGGGGAAGCTCAACCA‑ACACA‑AAACGGGA‑GGGCCGGCGAAAGGGCCACCCG‑TG‑CTAAACACCGCCATTGATGCGGccaga                                                                                                                                                     >  1:219549/1‑245 (MQ=255)
                                                                                                aaGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGGTTGCCTGGAACACCTAGAGCGAGATGGGTCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCCGCCGATG‑CTGGACACCGCAATTGATGCGGCAGAgatg                                                                                                                                                  >  1:249470/1‑251 (MQ=255)
                                                                                                      aGTGCGACCTGATCCAAGGTAAAACCGCGCCCCACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGACTGGAACACCCAGAGCGAGGGGGAGCTGCTGCGGAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGCAG‑CGGAAAAACACAATTTATGCGGCAGAgatgagcccg                                                                                                                                            >  2:70104/1‑246 (MQ=255)
                                                                                                       gTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGa                                                                                                                                           <  1:378581/251‑1 (MQ=255)
                                                                                                       gTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGa                                                                                                                                           <  1:49680/251‑1 (MQ=255)
                                                                                                         gCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGac                                                                                                                                          <  1:241629/250‑1 (MQ=255)
                                                                                                          cGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGacac                                                                                                                                        <  1:14453/251‑1 (MQ=255)
                                                                                                            cccTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTg                                                                                                                                      <  2:161908/250‑1 (MQ=255)
                                                                                                                    ccAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGAATGCCTTGGACACCAAGGACGGGGAGGAACGGATGCGCTAGATCAACTACAC‑CC‑AAACGGAAAGGACCG‑CGCAAGG‑CCAGACGAAG‑GTGAAAACCGCAATTGATGCGGgagaagagagccttgcacgggccccggg                                                                                                                              >  2:322786/1‑227 (MQ=255)
                                                                                                                          aaaaCCGCGCCACACATTATGCTCGTAGAGCGCGATTATCAGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGGAATGGCGGAAAAGGGGTTTCCTGGAACACCCCGAGAGAGATGGGTTTGGTGCGTAAGGTCAACAAAAC‑CCCAACAGGAA‑G‑ACCGGCGAAAGG‑CCAACCGATT‑CTGAACACCGCAAATAATGCGACAGAGATGGTACTGGCACAGG‑CACCGGGAACCa                                                                                                                        >  2:487306/1‑250 (MQ=255)
                                                                                                                            aaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGGATGCCTGGGACACCTTGAGAGAGGTGGGTCTGCTGCGGAAGCCCCACTAAAC‑CA‑AAGCGGGA‑AGTCCGGCGACAGG‑GCAGCCGGTGGCT‑AACACCCCAATTGGTGCGGCCGGGAAGAACCTGAcaccgg‑caccggggaacac                                                                                                                       >  2:250944/1‑242 (MQ=255)
                                                                                                                                gcgacacacaGTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGATAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGATCACCGCAATTGATGCGGCAGAGATGAGACTGACACTGG‑CACCGGAAACCAACGG‑TCa                                                                                                                <  2:270753/247‑1 (MQ=255)
                                                                                                                                 cgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAAACAACGG‑TCAg                                                                                                               <  2:251518/251‑1 (MQ=255)
                                                                                                                                   ccACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGt                                                                                                             <  1:25564/251‑1 (MQ=255)
                                                                                                                                           tATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTg                                                                                                      <  1:18741/250‑1 (MQ=255)
                                                                                                                                                tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGGACACCCAGAGCGCGAAGGGTCTGCTGCGTAAACTCCACTACCC‑CA‑AAACGGGA‑GGGCCGGCGGAAGG‑CCAGCCCAAG‑CTGAACAAACCAATTGATGCTGGAGGGATGATCCTGCCCCTGGGCACGGGAAAcaaaggtccgggaagcggggaaagc                                                                                                  >  2:262862/1‑227 (MQ=255)
                                                                                                                                                        agcgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGAGCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGac                                                                                                                                          <  2:33837/202‑1 (MQ=255)
                                                                                                                                                             ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAACTAGAGCGGGATGGATCTGCTGCGTAAGCTCAACAACAC‑CA‑AAACGGAA‑GGTCCGGCGAAAGG‑CCAGCCGAAG‑CTGGAAACCGCCATTGGTGCGGCAGAGAAGGTCCTGACAACGG‑GACCCGAAACCAACGG‑TCAGGTAGCCGTGGAAGCCTGGGCTGccccgg                                                                                   >  1:134991/1‑248 (MQ=255)
                                                                                                                                                                   ggatacgtACGAACGCCTTCCCTCTATCTGGCCGCTGATGGAGAGAACCGGTAATGGCGGGAAAGTGAATGCCTGGAACACCTGGAGCGTGATTGATCTGCTCCGTGAGCGCAACTGCCC‑CT‑AATCTGAA‑GGTCCGGCGTAAGT‑CTTTCCGATG‑CTGATCTCCGCAATGGATGGGGCGCAGATGTTCCTAACACTGG‑CGCCGTAGACCAACGG‑TCAGTTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAAt                                                                             <  1:547588/244‑1 (MQ=255)
                                                                                                                                                                      gACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGATTGCCCGGAAAACCTAGAGAGAGAAGGGTCTTCTGCGTAAGCTCAACTACAC‑CC‑AAGCGGGA‑GGTACGGCGGAAGG‑GCAGCCGATG‑CCGAAAACCGCAAATGATGGGGCAGAGATGATCCTGGCACTGG‑CACCGGGAAACAACGGGCCGGGAAGCCG‑GAAAGCCTGGGCTGGCCTgggagcattt                                                                           >  1:171501/1‑241 (MQ=255)
                                                                                                                                                                       ccTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTAc                                                                         <  2:527351/250‑1 (MQ=255)
                                                                                                                                                                               aCGCTTGACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTg                                                                 <  2:295719/251‑1 (MQ=255)
                                                                                                                                                                                            aacggtccGCTGACGGAGAAAATCGGTTATGGCGGTAAAGGGCTGGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCTTAAGCTCAACTACTC‑CA‑ATGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAGCACCGTAATTGATGCGGCAGAGAGGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGTCCATACGCATCt                                                    <  1:194091/245‑1 (MQ=255)
                                                                                                                                                                                            aTCGGCCCGTTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGAGCTGCTGCGTAAGCTCAACTACAC‑CA‑AATCGGAA‑GGTCCGGAGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACGGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATca                                                    <  2:199302/251‑2 (MQ=255)
                                                                                                                                                                                            aTCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATc                                                     <  2:249470/250‑1 (MQ=255)
                                                                                                                                                                                                                          ggCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAGACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATcc                       <  1:250944/250‑1 (MQ=255)
                                                                                                                                                                                                                           gtgggaaaGGGATTGCCGGGAACACCTAGATCGAGATGGATCTGCGGCGTAAGCTTAACAACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCTATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCTTGAGCGAATTTACCGGACGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCga                     <  2:456402/246‑2 (MQ=255)
                                                                                                                                                                                                                                     gATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGat            >  1:386301/1‑250 (MQ=255)
                                                                                                                                                                                                                                         gCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTc        >  2:187981/1‑250 (MQ=255)
                                                                                                                                                                                                                                              gAACACCTAGATCGAGATGGATCTGCTGCTTAAGCTCAACTCCAC‑CA‑AATCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGcac  <  2:286407/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGATTCCGCGATATTCAGGcac  >  2:354255/1‑251 (MQ=255)
                                                                                                                                                                                                                                                acacCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑ACGCGGAC‑GGTCCGGCGAAAGG‑CAAGCCGATG‑CTGAACACCGCCATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGcac  >  1:279625/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 cacCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑TCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGAc                                 <  1:388619/217‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          
GCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACAC‑CA‑AAGCGGAA‑GGTCCGGCGAAAGG‑CCAGCCGATG‑CTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGG‑CACCGGAAACCAACGG‑TCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCAC  >  NC_000913/1282382‑1282870

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑