Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 343,467 A→C 100% H195P (CAC→CCC)  yahK → broad specificity NADPH‑dependent aldehyde reductase, Zn‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913343,4670AC53.9% 33.1 / 110.3 87H195P (CAC→CCC) yahKbroad specificity NADPH‑dependent aldehyde reductase, Zn‑containing
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base C (23/24);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

TATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAAT  >  NC_000913/343373‑343564
                                                                                              |                                                                                                 
tATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGa                                                                                              <  2:1107464/100‑1 (MQ=255)
  ttCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATg                                                                                            <  1:1190348/100‑1 (MQ=255)
   ttgccgcTACGTCTCTGGCAGGCCGGGCCGGCTAAAATAGTGCGCGTGGTCGTCATCGGCGGTCTTGGACATTTTGGGATTAAGCTGGCCCAcgcg                                                                                               <  2:2096941/94‑1 (MQ=255)
    cgccgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATggg                                                                                          >  1:491749/1‑100 (MQ=255)
    cgccgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATggg                                                                                          >  1:968387/1‑100 (MQ=255)
    cgccgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCAGCCTGGGACACATGGGGATTAAGCTGGCCCCcgc                                                                                                >  2:397341/1‑94 (MQ=255)
      ccgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATggggg                                                                                        <  1:1055311/100‑1 (MQ=255)
          tACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGcaca                                                                                    >  1:901326/1‑100 (MQ=255)
          tACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGcaca                                                                                    >  2:1779475/1‑100 (MQ=255)
           aCGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACAt                                                                                   <  1:828157/100‑1 (MQ=255)
               cACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACAtgtg                                                                                >  1:1881736/1‑99 (MQ=255)
                 cTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATgtggtg                                                                             >  2:2016575/1‑100 (MQ=255)
                       ggtcgggccgggTAAAAAAGTGGGCGGGGTCGGCAGCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCAttt                                                                       <  1:537999/97‑1 (MQ=255)
                       ggccgggccgGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCAttt                                                                       >  1:1809746/1‑100 (MQ=255)
                       ggccgggccgGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCAttt                                                                       >  1:751711/1‑100 (MQ=255)
                        gccgggccggGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTa                                                                      >  2:1411596/1‑100 (MQ=255)
                         ccgggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTAc                                                                     <  2:1039359/100‑1 (MQ=255)
                          cgggccgggTCCAAAAGTGGGAGTGGTCGGTCTAGGCGGTGTGGGCCATATGGGGGTTAAGCTGGCACCAGAGATGGGGGCACATGTGGTGGAATTTAca                                                                    >  2:478806/1‑99 (MQ=255)
                            ggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTAcc                                                                    <  1:879837/98‑1 (MQ=255)
                            ggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCAc                                                                  <  2:1721602/100‑1 (MQ=255)
                            ggccgggTAAAAAAGTGGGCGTGGTCGGCATCGGCAGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATAGGGGCACATGTGGTGGCATTTa                                                                      <  2:717670/96‑1 (MQ=255)
                                 ggTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTg                                                             >  2:1919290/1‑100 (MQ=255)
                                  gTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGa                                                            <  2:81517/100‑1 (MQ=255)
                                   tAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAg                                                           <  1:983570/100‑1 (MQ=255)
                                     aaaaaGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCACCCGCGCTGGGGGCACATATGGTGGCATTTACCACTTCTGAGGc                                                         >  2:2113449/1‑100 (MQ=255)
                                      aaaaGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCa                                                        <  1:728874/100‑1 (MQ=255)
                                      aaaaGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCa                                                        >  1:1646648/1‑100 (MQ=255)
                                         aGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCaaaa                                                     >  1:35554/1‑100 (MQ=255)
                                             ggCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAAcgcg                                                 <  2:1588475/100‑1 (MQ=255)
                                             ggCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAAcgcg                                                 >  1:1397903/1‑100 (MQ=255)
                                              gCGTGGTCGGCATCGACGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTCGCATTTACCCCTTCTGAGACAAAAcgc                                                  >  2:610557/1‑98 (MQ=255)
                                               cGTGGTCGGCATCGGCTGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGaa                                               <  2:1400174/100‑1 (MQ=255)
                                                  ggTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCg                                            <  2:1412347/100‑1 (MQ=255)
                                                  ggTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCg                                            <  1:407620/100‑1 (MQ=255)
                                                    tCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGc                                          >  2:35015/1‑100 (MQ=255)
                                                    tCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGc                                          <  1:2291118/100‑1 (MQ=255)
                                                     cGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCa                                         <  2:1904196/100‑1 (MQ=255)
                                                     cGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCa                                         >  1:1683527/1‑100 (MQ=255)
                                                          tCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGc                                    >  2:1166608/1‑100 (MQ=255)
                                                          tCGGCGGTCTGGGACATAAAAGAAATAAGCAGGAACACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCaaaacc                                    >  2:2086742/1‑98 (MQ=255)
                                                            ggcggTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGccc                                  <  1:1835701/100‑1 (MQ=255)
                                                            ggcggTCTGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGcc                                   >  2:1560011/1‑99 (MQ=255)
                                                            ggcggTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGccc                                  >  1:1581390/1‑100 (MQ=255)
                                                              cggTCTGCGACATATGGGGATGATGGCGGCCCCCGCGTTGGGGGCATATGTTGTGGCATTCACCACTGCTGAGGCAAAGCTCGAAGCGGCAAAAGccccg                                <  2:1098069/100‑3 (MQ=255)
                                                               ggTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGTCAAAACGCGAAGCGGCAAAAGCCCTg                                <  1:2029003/99‑1 (MQ=255)
                                                               ggTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTgg                               >  2:796036/1‑100 (MQ=255)
                                                                gTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTggg                              >  1:1381626/1‑100 (MQ=255)
                                                                 tCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTgggg                             <  1:595571/100‑1 (MQ=255)
                                                                  cTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTggg                              >  2:946559/1‑98 (MQ=255)
                                                                   tgtgACATATGGGGATTAAGCTAGCCAACGGAATGGGGGCACATGTGGTGACATATCCCACTTCTGAAGCAAAACGCGAAGCTGCAATAGCCCTGGGGGc                           <  2:2247924/97‑1 (MQ=255)
                                                                   tGGGACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGc                           >  2:729566/1‑100 (MQ=255)
                                                                    gggACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGcc                          >  1:1801421/1‑100 (MQ=255)
                                                                    gggACATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGcc                          >  1:186761/1‑100 (MQ=255)
                                                                      gACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGAGGCCGa                        >  1:937603/1‑100 (MQ=255)
                                                                        cATATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGACGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATg                      >  1:1652376/1‑100 (MQ=255)
                                                                        cATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGACATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATg                      >  2:1643389/1‑100 (MQ=255)
                                                                         atatGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGa                     <  1:1338030/100‑1 (MQ=255)
                                                                          tatGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGaa                    >  1:1386681/1‑100 (MQ=255)
                                                                          tatGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGaa                    >  2:277746/1‑100 (MQ=255)
                                                                          tatGGGGATTAAGCTGGCACCAGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGAtaa                     >  2:973670/1‑97 (MQ=255)
                                                                            tGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgt                  >  1:1444331/1‑100 (MQ=255)
                                                                            tGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgt                  >  1:563406/1‑100 (MQ=255)
                                                                            tGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgt                  <  2:627640/100‑1 (MQ=255)
                                                                              gggATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttg                <  2:2063845/100‑1 (MQ=255)
                                                                              gggATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttg                <  1:919493/100‑1 (MQ=255)
                                                                              gggATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttg                <  2:1381612/100‑1 (MQ=255)
                                                                              gggATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttg                >  2:2149093/1‑100 (MQ=255)
                                                                               ggATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttgt               >  2:310741/1‑100 (MQ=255)
                                                                               ggATTAAGCTGGCCCACGCGATGGGGGAACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttgt               >  2:311052/1‑100 (MQ=255)
                                                                                gATTAAGCTGGCCCACGGGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttgtt              <  2:696187/100‑1 (MQ=255)
                                                                                gATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttgtt              <  1:2105301/100‑1 (MQ=255)
                                                                                gATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttgtt              <  1:2105331/100‑1 (MQ=255)
                                                                                gATTAAGCGGGCCCCCGTGATGGGGGCACATGTGGTGTCATTTAAAAAAAATAAAAAAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttgtt              <  2:187445/100‑1 (MQ=255)
                                                                                 aTTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTa             <  2:1987706/100‑1 (MQ=255)
                                                                                  ttAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttgtt              <  1:493673/98‑1 (MQ=255)
                                                                                   tAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAAc           <  2:35553/100‑1 (MQ=255)
                                                                                   tAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAAc           <  2:1397881/100‑1 (MQ=255)
                                                                                   tAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAAc           <  1:929622/100‑1 (MQ=255)
                                                                                   tAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAAc           <  2:161175/100‑1 (MQ=255)
                                                                                    aaGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACt          <  1:1166619/100‑1 (MQ=255)
                                                                                    aaGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACt          <  1:472936/100‑1 (MQ=255)
                                                                                     aGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACt          >  1:2086229/1‑99 (MQ=255)
                                                                                         ggCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACg      >  1:607524/1‑99 (MQ=255)
                                                                                          gctcACGCGGTGGCGGCACATGTGGTGGCATTTACCACTTCTGAGGAGAAACGCGAAGCGGCAAAAGCCTTGGGGGCCGATTAAGTTGTTAACTCACga     <  2:642289/96‑2 (MQ=255)
                                                                                          gCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCa    <  1:1422017/100‑1 (MQ=255)
                                                                                           ccccCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCaa   >  1:609290/1‑100 (MQ=255)
                                                                                           ccccCGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCaa   >  1:1757842/1‑100 (MQ=255)
                                                                                           cccACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCaa   >  1:860279/1‑100 (MQ=255)
                                                                                            ccACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAAt  <  1:1738358/100‑1 (MQ=255)
                                                                                              |                                                                                                 
TATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAAT  >  NC_000913/343373‑343564

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