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

GCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCAC  >  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)
                                     tATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAAAGCGAGGTGGAACTGCTGCGTAAGCCCCACTACAcccaa                                                                                                                                                                                                                  >  1:299014/1‑247 (MQ=255)
                                      aTCCAGCATGACGCTGCCGCTGAACGGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGc                                                                                                                                                                                                               <  2:118072/251‑1 (MQ=255)
                                            cATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGAAAAGGGGATGCCTGGAACACCTAAAGCGGGGGGGATCTGCTGAGACAGCTCAAATAAACCAAAGCGGCAgg                                                                                                                                                                                                        >  1:497864/1‑251 (MQ=255)
                                                                   gcAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAACGGGGTTTCCTGGGACACCTCGAGCGAGAAGGGGCGGTTGCGGAAGCTCAACTACACCAAAGCGGGAGGAACGGCGCAAGGCCAGCCGATgg                                                                                                                                                                                >  2:19175/1‑250 (MQ=255)
                                                                       ccGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGaa                                                                                                                                                                           <  1:469633/251‑1 (MQ=255)
                                                                               tgtgAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGa                                                                                                                                                                                                            <  2:32922/213‑1 (MQ=255)
                                                                               tgtgAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGGTGGATCTGGTGCGGAAGCTCAACTACACCAAAGCggg                                                                                                                                                                                                            >  1:32922/1‑212 (MQ=255)
                                                                                          ggAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCGTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGc                                                                                                                                                         <  1:578423/250‑1 (MQ=255)
                                                                                           gAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCGCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGTAGTTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCag                                                                                                                                                       <  1:226942/251‑1 (MQ=255)
                                                                                               aaaGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGATTGCCTGGGAAACCTAAAACGGGATGGATCTGCTGGGGAAGCTCAACCAACACAAAACGGGAGGGCCGGCGAAAGGGCCACCCGTGCTAAACACCGCCATTGATGCGGccaga                                                                                                                                                     >  1:219549/1‑245 (MQ=255)
                                                                                                aaGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGGTTGCCTGGAACACCTAGAGCGAGATGGGTCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCCGCCGATGCTGGACACCGCAATTGATGCGGCAGAgatg                                                                                                                                                  >  1:249470/1‑251 (MQ=255)
                                                                                                      aGTGCGACCTGATCCAAGGTAAAACCGCGCCCCACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGACTGGAACACCCAGAGCGAGGGGGAGCTGCTGCGGAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGCAGCGGAAAAACACAATTTATGCGGCAGAgatgagcccg                                                                                                                                            >  2:70104/1‑246 (MQ=255)
                                                                                                       gTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGa                                                                                                                                           <  1:378581/251‑1 (MQ=255)
                                                                                                       gTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGa                                                                                                                                           <  1:49680/251‑1 (MQ=255)
                                                                                                         gCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGac                                                                                                                                          <  1:241629/250‑1 (MQ=255)
                                                                                                          cGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGacac                                                                                                                                        <  1:14453/251‑1 (MQ=255)
                                                                                                            cccTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTg                                                                                                                                      <  2:161908/250‑1 (MQ=255)
                                                                                                                    ccAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGAATGCCTTGGACACCAAGGACGGGGAGGAACGGATGCGCTAGATCAACTACACCCAAACGGAAAGGACCGCGCAAGGCCAGACGAAGGTGAAAACCGCAATTGATGCGGgagaagagagccttgcacgggccccggg                                                                                                                              >  2:322786/1‑227 (MQ=255)
                                                                                                                          aaaaCCGCGCCACACATTATGCTCGTAGAGCGCGATTATCAGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGGAATGGCGGAAAAGGGGTTTCCTGGAACACCCCGAGAGAGATGGGTTTGGTGCGTAAGGTCAACAAAACCCCAACAGGAAGACCGGCGAAAGGCCAACCGATTCTGAACACCGCAAATAATGCGACAGAGATGGTACTGGCACAGGCACCGGGAACCa                                                                                                                        >  2:487306/1‑250 (MQ=255)
                                                                                                                            aaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGGATGCCTGGGACACCTTGAGAGAGGTGGGTCTGCTGCGGAAGCCCCACTAAACCAAAGCGGGAAGTCCGGCGACAGGGCAGCCGGTGGCTAACACCCCAATTGGTGCGGCCGGGAAGAACCTGAcaccggcaccggggaacac                                                                                                                       >  2:250944/1‑242 (MQ=255)
                                                                                                                                gcgacacacaGTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGATAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGATCACCGCAATTGATGCGGCAGAGATGAGACTGACACTGGCACCGGAAACCAACGGTCa                                                                                                                <  2:270753/247‑1 (MQ=255)
                                                                                                                                 cgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAAACAACGGTCAg                                                                                                               <  2:251518/251‑1 (MQ=255)
                                                                                                                                   ccACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGt                                                                                                             <  1:25564/251‑1 (MQ=255)
                                                                                                                                           tATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTg                                                                                                      <  1:18741/250‑1 (MQ=255)
                                                                                                                                                tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGGACACCCAGAGCGCGAAGGGTCTGCTGCGTAAACTCCACTACCCCAAAACGGGAGGGCCGGCGGAAGGCCAGCCCAAGCTGAACAAACCAATTGATGCTGGAGGGATGATCCTGCCCCTGGGCACGGGAAAcaaaggtccgggaagcggggaaagc                                                                                                  >  2:262862/1‑227 (MQ=255)
                                                                                                                                                        agcgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGAGCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGac                                                                                                                                          <  2:33837/202‑1 (MQ=255)
                                                                                                                                                             ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAACTAGAGCGGGATGGATCTGCTGCGTAAGCTCAACAACACCAAAACGGAAGGTCCGGCGAAAGGCCAGCCGAAGCTGGAAACCGCCATTGGTGCGGCAGAGAAGGTCCTGACAACGGGACCCGAAACCAACGGTCAGGTAGCCGTGGAAGCCTGGGCTGccccgg                                                                                   >  1:134991/1‑248 (MQ=255)
                                                                                                                                                                   ggatacgtACGAACGCCTTCCCTCTATCTGGCCGCTGATGGAGAGAACCGGTAATGGCGGGAAAGTGAATGCCTGGAACACCTGGAGCGTGATTGATCTGCTCCGTGAGCGCAACTGCCCCTAATCTGAAGGTCCGGCGTAAGTCTTTCCGATGCTGATCTCCGCAATGGATGGGGCGCAGATGTTCCTAACACTGGCGCCGTAGACCAACGGTCAGTTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAAt                                                                             <  1:547588/244‑1 (MQ=255)
                                                                                                                                                                      gACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGATTGCCCGGAAAACCTAGAGAGAGAAGGGTCTTCTGCGTAAGCTCAACTACACCCAAGCGGGAGGTACGGCGGAAGGGCAGCCGATGCCGAAAACCGCAAATGATGGGGCAGAGATGATCCTGGCACTGGCACCGGGAAACAACGGGCCGGGAAGCCGGAAAGCCTGGGCTGGCCTgggagcattt                                                                           >  1:171501/1‑241 (MQ=255)
                                                                                                                                                                       ccTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTAc                                                                         <  2:527351/250‑1 (MQ=255)
                                                                                                                                                                               aCGCTTGACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTg                                                                 <  2:295719/251‑1 (MQ=255)
                                                                                                                                                                                            aacggtccGCTGACGGAGAAAATCGGTTATGGCGGTAAAGGGCTGGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCTTAAGCTCAACTACTCCAATGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAGCACCGTAATTGATGCGGCAGAGAGGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGTCCATACGCATCt                                                    <  1:194091/245‑1 (MQ=255)
                                                                                                                                                                                            aTCGGCCCGTTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGAGCTGCTGCGTAAGCTCAACTACACCAAATCGGAAGGTCCGGAGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACGGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATca                                                    <  2:199302/251‑2 (MQ=255)
                                                                                                                                                                                            aTCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATc                                                     <  2:249470/250‑1 (MQ=255)
                                                                                                                                                                                                                          ggCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAGACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATcc                       <  1:250944/250‑1 (MQ=255)
                                                                                                                                                                                                                           gtgggaaaGGGATTGCCGGGAACACCTAGATCGAGATGGATCTGCGGCGTAAGCTTAACAACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCTATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCTTGAGCGAATTTACCGGACGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCga                     <  2:456402/246‑2 (MQ=255)
                                                                                                                                                                                                                                     gATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGat            >  1:386301/1‑250 (MQ=255)
                                                                                                                                                                                                                                         gCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTc        >  2:187981/1‑250 (MQ=255)
                                                                                                                                                                                                                                              gAACACCTAGATCGAGATGGATCTGCTGCTTAAGCTCAACTCCACCAAATCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGcac  <  2:286407/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGATTCCGCGATATTCAGGcac  >  2:354255/1‑251 (MQ=255)
                                                                                                                                                                                                                                                acacCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAACGCGGACGGTCCGGCGAAAGGCAAGCCGATGCTGAACACCGCCATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGcac  >  1:279625/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 cacCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGTCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGAc                                 <  1:388619/217‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                          
GCCATCTGGGTAAAGAAACCGACATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCAC  >  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: