Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,490,864 T→C Q19R (CAG→CGG)  idnO ← 5‑keto‑D‑gluconate 5‑reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,490,8640TC98.4% 1021.8 / 2.6 307Q19R (CAG→CGG) idnO5‑keto‑D‑gluconate 5‑reductase
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base C (128/174);  total (132/175)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 336 total aligned reads displayed.

GCAAGTTCTGCGCGTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCTGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTATCTCTTAGAAAACAA  >  NC_000913/4490774‑4490953
                                                                                          |                                                                                         
gCAAGTTCTGCGCGTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCcgt                                                                                         <  1:5767586/93‑1 (MQ=255)
gCAAGTTCTGCGCGTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCcgt                                                                                         >  1:2943434/1‑93 (MQ=255)
 cAAGTTCTGCGCGTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCcgtg                                                                                        <  1:4774540/93‑1 (MQ=255)
 cAAGTTCTGCGCGTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCcgtg                                                                                        >  1:5201932/1‑93 (MQ=255)
  aaGTTCTGCGCGTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGc                                                                                       <  2:862620/93‑1 (MQ=255)
   aGTTCTGCGCGTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCt                                                                                      <  2:5124826/93‑1 (MQ=255)
    gTTCTGCGCGTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCACGTGCTg                                                                                     >  1:7434938/1‑93 (MQ=255)
       cTGCGCGTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAAc                                                                                  >  2:3708505/1‑93 (MQ=255)
        tGCGCGTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAAcc                                                                                 <  2:7500552/93‑1 (MQ=255)
        tGCGCGTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAAcc                                                                                 <  1:2734901/93‑1 (MQ=255)
         gcgcGTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCg                                                                                >  2:156526/1‑93 (MQ=255)
          cgcgTTCGGCAGTAATATCATTCATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAACAGCCACTGCCCCGTGCTGAACCgg                                                                               >  2:5204688/1‑93 (MQ=255)
            cgTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTa                                                                             <  2:85114/93‑1 (MQ=255)
                cGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCa                                                                         >  2:3230966/1‑93 (MQ=255)
                 ggCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCaa                                                                        >  1:6254002/1‑93 (MQ=255)
                  gCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCTGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAg                                                                       >  1:8648339/1‑93 (MQ=255)
                  gCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAg                                                                       <  2:4172662/93‑1 (MQ=255)
                   cAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGa                                                                      <  2:5288945/93‑1 (MQ=255)
                   cAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGa                                                                      >  2:4792722/1‑93 (MQ=255)
                    aGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGat                                                                     >  2:1434123/1‑93 (MQ=255)
                     gTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGata                                                                    >  1:3078651/1‑93 (MQ=255)
                     gTAATATCATGAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGata                                                                    <  1:2515506/93‑1 (MQ=255)
                      tAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGatat                                                                   <  2:2562794/93‑1 (MQ=255)
                       aaTATGATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATAtt                                                                  >  2:8743352/1‑93 (MQ=255)
                       aaTATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATAtt                                                                  >  2:6974462/1‑93 (MQ=255)
                        atatCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATAttt                                                                 <  2:2726205/93‑1 (MQ=255)
                        atatCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATAttt                                                                 <  2:6040818/93‑1 (MQ=255)
                          atCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATAttttt                                                               <  2:3403246/93‑1 (MQ=255)
                          atCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATAttttt                                                               <  1:5685830/93‑1 (MQ=255)
                           tCATTAATAATTATTTGTGCGCCATATTTACCCAGGTCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATAtttttt                                                              >  2:763996/1‑93 (MQ=255)
                           tCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATAtttttt                                                              >  2:4862342/1‑93 (MQ=255)
                            cATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTc                                                             <  1:3530222/93‑1 (MQ=255)
                             aTTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTTCTTAACCGGTAATCAAGATATTTTTTcc                                                            <  2:1109682/93‑1 (MQ=255)
                             aTTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTcc                                                            <  1:3372781/93‑1 (MQ=255)
                              ttAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCt                                                           <  1:2171091/93‑1 (MQ=255)
                                 ataatTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGcc                                                        <  2:303246/93‑1 (MQ=255)
                                  taatTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCa                                                       <  1:8140998/93‑1 (MQ=255)
                                  taatTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCa                                                       >  2:6267175/1‑93 (MQ=255)
                                  taatTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCa                                                       <  1:375521/93‑1 (MQ=255)
                                   aatTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAg                                                      <  1:7572406/93‑1 (MQ=255)
                                    attattTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGt                                                     <  2:1914089/93‑1 (MQ=255)
                                    attattTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGt                                                     >  2:1305466/1‑93 (MQ=255)
                                     ttattTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCTGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTg                                                    <  1:1460710/93‑1 (MQ=255)
                                       ttttGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGaa                                                  <  2:3604369/92‑1 (MQ=255)
                                          tgtgCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAAt                                               <  2:5769881/93‑1 (MQ=255)
                                           gtgCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATa                                              <  1:1206652/93‑1 (MQ=255)
                                           gtgCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATa                                              <  1:4531838/93‑1 (MQ=255)
                                           gtgCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATa                                              <  1:1061522/93‑1 (MQ=255)
                                           gtgCGCAATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATa                                              <  2:7850880/93‑1 (MQ=255)
                                             gcgcCATATTTTCCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGa                                            <  1:7048599/93‑1 (MQ=255)
                                             gcgcCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGa                                            <  2:1526171/93‑1 (MQ=255)
                                             gcgcAATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGa                                            <  2:4356244/93‑1 (MQ=255)
                                              ggcCACATTTACCAAGGCGGGTTGCCAGTAAAAACCCAATGCCCCGTTCTGAACCGGTAATAAAGATTTTTTTTCCTGCCAGTGAAAATAGAt                                           <  1:4937467/92‑1 (MQ=255)
                                              cgcCATATTTACCCAGTCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGAt                                           <  2:2290356/93‑1 (MQ=255)
                                              cgcCATATTTACCCAGCCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGAt                                           <  2:5491770/93‑1 (MQ=255)
                                               gcCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATc                                          <  2:4966666/93‑1 (MQ=255)
                                                 cATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGt                                        <  1:4317213/93‑1 (MQ=255)
                                                 cATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGt                                        <  1:5535617/93‑1 (MQ=255)
                                                   tatTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTc                                      <  1:3593311/93‑1 (MQ=255)
                                                    atTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCa                                     <  1:440154/93‑1 (MQ=255)
                                                      ttACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCAtt                                   <  2:3769362/93‑1 (MQ=255)
                                                      ttACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCAtt                                   >  2:8150734/1‑93 (MQ=255)
                                                       tACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTg                                  <  1:2801235/93‑1 (MQ=255)
                                                       tACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTg                                  <  1:5314642/93‑1 (MQ=255)
                                                         cccAGGCCGGCTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGta                                >  2:6154972/1‑93 (MQ=255)
                                                            aGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGtatta                             >  2:3117376/1‑93 (MQ=255)
                                                            aGGCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGtatta                             <  1:2659568/93‑1 (MQ=255)
                                                             ggCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGtattat                            >  2:5670347/1‑93 (MQ=255)
                                                             ggCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGtattat                            >  2:5928976/1‑93 (MQ=255)
                                                             ggCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGtattat                            >  1:8481854/1‑93 (MQ=255)
                                                              gCCGGTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTattatt                           >  1:3817616/1‑93 (MQ=255)
                                                                 ggTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCt                        >  1:4013316/1‑93 (MQ=255)
                                                                 ggTTGCCAGTAAAAAGCAAATGCCCCGTGCTGAACCGGTAATAAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCt                        <  2:7346683/93‑1 (MQ=255)
                                                                  gTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGTCAGTGAAAATTGATCGTTCATTGTattatttttt                       >  2:4668713/1‑89 (MQ=255)
                                                                  gTTGCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCtt                       <  1:1918297/93‑1 (MQ=255)
                                                                     gCCAGTAAAAAGCCAATTCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCtttat                    <  1:872602/93‑1 (MQ=255)
                                                                     gCCAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCtttat                    <  2:6406843/93‑1 (MQ=255)
                                                                      ccAGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCtttatt                   <  1:6499652/93‑1 (MQ=255)
                                                                        aGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCtttattta                 >  1:7964865/1‑93 (MQ=255)
                                                                        aGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCtttattta                 <  1:8261862/93‑1 (MQ=255)
                                                                        aGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCtttattta                 >  2:8039111/1‑93 (MQ=255)
                                                                        aGTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCtttattta                 <  1:6165996/93‑1 (MQ=255)
                                                                         gTAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTttatttat                <  2:1772898/93‑1 (MQ=255)
                                                                         atAAAAAGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTttatttat                <  1:7639850/92‑1 (MQ=255)
                                                                           aaaaaTCCAATGCCCCGTGCTGACCCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTAtct              <  2:4113686/93‑1 (MQ=255)
                                                                           aaaaaGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTAtct              <  2:6615704/93‑1 (MQ=255)
                                                                           aaaaaGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTAtct              >  1:2312113/1‑93 (MQ=255)
                                                                           aaaaaGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTAtct              <  2:4000083/93‑1 (MQ=255)
                                                                            aaaaGCGAATGCCCAGTCCTGAACCGGGAATCAAGATATTTTTTGCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTAtctc             >  2:6481785/1‑93 (MQ=255)
                                                                             aaaGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTATctct            <  2:1079602/93‑1 (MQ=255)
                                                                              aaGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTATCTCtt           <  2:7085078/93‑1 (MQ=255)
                                                                               aGCCAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTATCTCTTa          <  1:8132346/93‑1 (MQ=255)
                                                                               aGCCAATGCCCAGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTATCTCTTa          <  1:1830248/93‑1 (MQ=255)
                                                                                gCCAATGCCCCGTGCTGAACCGGTGATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTATCTCTTAg         <  1:7813768/93‑1 (MQ=255)
                                                                                 ccAATGGCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTATCTCTTAGa        <  1:1692067/93‑1 (MQ=255)
                                                                                  cAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTATCTCTTAGaa       <  2:7300492/93‑1 (MQ=255)
                                                                                  cAATGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTATCTCTTAGaa       <  1:6704692/93‑1 (MQ=255)
                                                                                     tGCCCTGTGATGAACCGGTAATCAAGACATTTTTTCCTGCAAGTGAAAATAGATCGGTCATTGTCTTCTTCCTTTATTTATCTCTTAGaagaa    >  1:689088/1‑90 (MQ=255)
                                                                                     tGCCCCGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTATCTCTTAGAAAAc    <  2:211597/93‑1 (MQ=255)
                                                                                       ccccGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGACAATAGATCGTTCATTGTATTATTCCTTTATTTATCTCTTAGAAAAcca  >  1:50983/1‑91 (MQ=255)
                                                                                          |                                                                                         
GCAAGTTCTGCGCGTTCGGCAGTAATATCATTAATAATTATTTGTGCGCCATATTTACCCAGGCCGGTTGCCAGTAAAAAGCCAATGCCCTGTGCTGAACCGGTAATCAAGATATTTTTTCCTGCCAGTGAAAATAGATCGTTCATTGTATTATTCCTTTATTTATCTCTTAGAAAACAA  >  NC_000913/4490774‑4490953

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