Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 129,360 G→T intergenic (+24/+47) lpd → / ← yacH dihydrolipoyl dehydrogenase; E3 component of pyruvate and 2‑oxoglutarate dehydrogenases complexes; glycine cleavage system L protein; dihydrolipoamide dehydrogenase/DUF3300 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913129,3600GT97.1% 455.0 / NA 140intergenic (+24/+47)lpd/yacHdihydrolipoyl dehydrogenase; E3 component of pyruvate and 2‑oxoglutarate dehydrogenases complexes; glycine cleavage system L protein; dihydrolipoamide dehydrogenase/DUF3300 family protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  major base T (64/70);  minor base A (3/1);  total (69/71)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

CTGTGGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGGCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCATACTG  >  NC_000913/129263‑129456
                                                                                                  |                                                                                                
cTGTGGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCa                                                                                                <  1:640081/100‑1 (MQ=255)
cTGTGGGCCTGGCGGCAGAAGTGTGCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCa                                                                                                <  2:2773166/100‑1 (MQ=255)
 tgtgGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCaa                                                                                               >  2:3147805/1‑100 (MQ=255)
 tgtgGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCaa                                                                                               <  2:2690052/100‑1 (MQ=255)
  gtgGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAAt                                                                                              <  2:2653521/100‑1 (MQ=255)
   tgGGCCTGGCGGTAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAAtt                                                                                             >  2:3119581/1‑100 (MQ=255)
   tgGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAAtt                                                                                             >  2:1481692/1‑100 (MQ=255)
    gggCCTGGCGGCCGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTa                                                                                            <  1:2325271/100‑1 (MQ=255)
     ggCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCAGGACCATCCCTCAATTaa                                                                                           >  2:637925/1‑100 (MQ=255)
       ccTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTaaaa                                                                                         >  2:1493319/1‑100 (MQ=255)
        cTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTaaaaa                                                                                        >  1:1941545/1‑100 (MQ=255)
          ggcggcAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTGAATTgaaaa                                                                                        >  1:1736862/1‑98 (MQ=255)
          ggcggcAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAg                                                                                      >  2:3277073/1‑100 (MQ=255)
            cggcAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCg                                                                                    <  2:3068093/100‑1 (MQ=255)
             ggcAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCgg                                                                                   >  1:2191944/1‑100 (MQ=255)
              gcAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGc                                                                                  <  1:3252326/100‑1 (MQ=255)
              gcAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTATCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGc                                                                                  <  2:2448233/100‑1 (MQ=255)
               cAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAAGCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCt                                                                                 >  1:1850809/1‑100 (MQ=255)
                aGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTa                                                                                <  1:1674354/100‑1 (MQ=255)
                  aaGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAAc                                                                              >  2:1314749/1‑100 (MQ=255)
                   aGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAAcc                                                                             >  1:1094142/1‑100 (MQ=255)
                   aGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAAcc                                                                             >  2:2640035/1‑100 (MQ=255)
                   aGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAAcc                                                                             >  1:465472/1‑100 (MQ=255)
                   aGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAAcc                                                                             <  1:452132/100‑1 (MQ=255)
                   aGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAAcc                                                                             >  1:1473222/1‑100 (MQ=255)
                      gtTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcg                                                                          >  1:540131/1‑100 (MQ=255)
                      gtTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcg                                                                          >  2:1157710/1‑100 (MQ=255)
                      gtTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcg                                                                          <  1:895901/100‑1 (MQ=255)
                       ttCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAATTTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAc                                                                           <  1:624288/99‑1 (MQ=255)
                       ttCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGcgtggttagc                                                                         <  1:1189528/100‑10 (MQ=255)
                       ttCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcgc                                                                         >  2:2565538/1‑100 (MQ=255)
                       ttCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcgc                                                                         >  2:3468289/1‑100 (MQ=255)
                       ttCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcgc                                                                         >  2:1653968/1‑100 (MQ=255)
                       ttCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcgc                                                                         >  2:1625586/1‑100 (MQ=255)
                        tCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcgcc                                                                        >  2:2436748/1‑100 (MQ=255)
                         cGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcgccg                                                                       >  1:192370/1‑100 (MQ=255)
                          gAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCAcgccgc                                                                      >  2:2165732/1‑100 (MQ=255)
                           aaGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCt                                                                     <  2:2732524/100‑1 (MQ=255)
                            aGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGTCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCtt                                                                    >  1:2432556/1‑100 (MQ=255)
                            aGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCtt                                                                    <  2:3404080/100‑1 (MQ=255)
                              gTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCtttt                                                                  >  1:2863283/1‑100 (MQ=255)
                                 gCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCttttttt                                                               <  1:900551/100‑1 (MQ=255)
                                 gCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCttttttt                                                               <  2:1831431/100‑1 (MQ=255)
                                  cATTACCGACCTGCCGAGCCCGAAAGCGAAGAAGCAGGAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTa                                                              <  2:621834/100‑1 (MQ=255)
                                  cATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTa                                                              <  2:2625226/100‑1 (MQ=255)
                                  cATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTa                                                              <  2:3205519/100‑1 (MQ=255)
                                  cATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTa                                                              <  2:451004/100‑1 (MQ=255)
                                  cATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTa                                                              <  2:615810/100‑1 (MQ=255)
                                   aTTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTAc                                                             <  1:3178937/100‑1 (MQ=255)
                                    ttACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACg                                                            <  1:2785044/100‑1 (MQ=255)
                                      atcGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTc                                                          <  1:2101879/98‑1 (MQ=255)
                                        cGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTg                                                        <  2:2211377/100‑1 (MQ=255)
                                         gACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGc                                                       >  2:2438713/1‑100 (MQ=255)
                                              gCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAAttt                                                  >  2:2621287/1‑100 (MQ=255)
                                               ccGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGGCAATTAAAAAAGCGGCTAACAACGCCGCTTTTTTTACGTCTGCAAt                                                    >  1:87460/1‑97 (MQ=255)
                                                cGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTAc                                                >  2:150460/1‑100 (MQ=255)
                                                cGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTAc                                                >  1:2041859/1‑100 (MQ=255)
                                                    cccGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCttt                                            >  2:802102/1‑100 (MQ=255)
                                                     ccGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTc                                           >  2:2829125/1‑100 (MQ=255)
                                                     ccGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTc                                           <  2:2530021/100‑1 (MQ=255)
                                                     ccGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTc                                           >  1:1529541/1‑100 (MQ=255)
                                                      cGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTAccagatc                                          >  2:1878343/1‑95 (MQ=255)
                                                      cGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTcc                                          <  1:1242336/100‑1 (MQ=255)
                                                      cGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTcc                                          >  2:1005437/1‑100 (MQ=255)
                                                       gAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCa                                         <  2:3398613/100‑1 (MQ=255)
                                                        aaaGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAg                                        <  2:3023866/100‑1 (MQ=255)
                                                        aaaGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAg                                        <  1:2616985/100‑1 (MQ=255)
                                                        aaaGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAg                                        <  1:456614/100‑1 (MQ=255)
                                                          aGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtc                                      >  1:3265238/1‑100 (MQ=255)
                                                            cGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtctt                                    >  2:2308393/1‑100 (MQ=255)
                                                            cGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtctt                                    >  2:1301607/1‑100 (MQ=255)
                                                              aagaagaagTAA‑TTTTTTTTTTTTTTTTACATCCGACAATTAAAAAAGCAGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtcttct                                  >  1:712650/1‑100 (MQ=255)
                                                              aagaagaagTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtcttct                                  <  1:2882690/100‑1 (MQ=255)
                                                              aagaagaagTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtcttct                                  >  2:3191255/1‑100 (MQ=255)
                                                              aagaagaagTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGtcttct                                  <  2:2143043/100‑1 (MQ=255)
                                                                 aagaagTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGc                               <  1:802091/100‑1 (MQ=255)
                                                                 aagaagTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGc                               <  2:3425159/100‑1 (MQ=255)
                                                                  agaagTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCt                              >  2:740737/1‑100 (MQ=255)
                                                                     agTAA‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCa                           >  1:1603998/1‑100 (MQ=255)
                                                                       tAA‑TTTTACGTATGCCGGAACATCCGACAATTAAAAAAGAGGCTAACCACGCCACATTTATTACGTCTGCAATTGACATTTCCAGTCTACTTGTACCACg                         >  1:1411418/1‑100 (MQ=255)
                                                                         a‑TTTTTCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGtt                       >  1:3187936/1‑100 (MQ=255)
                                                                            ttttCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCa                     >  1:2895462/1‑100 (MQ=255)
                                                                             tttCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCag                    <  2:3004124/100‑1 (MQ=255)
                                                                             tttCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCag                    <  2:1993844/100‑1 (MQ=255)
                                                                              ttCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCaga                   <  1:822739/100‑1 (MQ=255)
                                                                               tCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCggag                  <  1:3220724/100‑1 (MQ=255)
                                                                               tCGTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCagag                  <  2:2322862/100‑1 (MQ=255)
                                                                                 gTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCagagag                <  2:540137/100‑1 (MQ=255)
                                                                                 gTTTGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCagagag                >  1:2709212/1‑100 (MQ=255)
                                                                                  tttGCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAgagaga               <  1:1314730/100‑1 (MQ=255)
                                                                                     gCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGt            >  1:15537/1‑100 (MQ=255)
                                                                                     gCCGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGt            >  2:600532/1‑100 (MQ=255)
                                                                                      ccGGAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGtt           >  1:1446534/1‑100 (MQ=255)
                                                                                        ggAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCg         >  1:1686985/1‑100 (MQ=255)
                                                                                        ggAACATCCGTCAATTAAAAAAGCGGCTAACCACGCCCCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCg         <  2:2076176/100‑1 (MQ=255)
                                                                                         gAACATCCGTCAATTAGAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCGTTCCATTCTTCTTGCTCCACGTTCAGAGAGACGTTCg         <  2:3318099/99‑1 (MQ=255)
                                                                                          aaCATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGc        <  2:1155576/99‑1 (MQ=255)
                                                                                          aaCATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCa       <  2:3324045/100‑1 (MQ=255)
                                                                                           aCATCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCAt      <  1:1036823/100‑1 (MQ=255)
                                                                                              tCCGTCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCATAct   <  2:177103/100‑1 (MQ=255)
                                                                                                  |                                                                                                
CTGTGGGCCTGGCGGCAGAAGTGTTCGAAGGTAGCATTACCGACCTGCCGAACCCGAAAGCGAAGAAGAAGTAA‑TTTTTCGTTTGCCGGAACATCCGGCAATTAAAAAAGCGGCTAACCACGCCGCTTTTTTTACGTCTGCAATTTACCTTTCCAGTCTTCTTGCTCCACGTTCAGAGAGACGTTCGCATACTG  >  NC_000913/129263‑129456

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