Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,524,976 G→A S12F (TCC→TTC)  yjhH ← KpLE2 phage‑like element; putative 2‑dehydro‑3‑deoxy‑D‑pentonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,524,9760GA99.8% 1886.2 / ‑4.7 434S12F (TCC→TTC) yjhHKpLE2 phage‑like element; putative 2‑dehydro‑3‑deoxy‑D‑pentonate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (227/205);  total (227/207)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 434 total aligned reads displayed.

CCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGGATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCAGCCCATATGATATTT  >  NC_000913/4524878‑4525073
                                                                                                  |                                                                                                 
ccAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAAt                                                                                                 <  2:8791815/101‑1 (MQ=255)
ccAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAAt                                                                                                 >  2:8792621/1‑101 (MQ=255)
 cAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATa                                                                                                <  2:8791882/101‑1 (MQ=255)
 cAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATa                                                                                                >  2:8792625/1‑101 (MQ=255)
  aGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATAc                                                                                               >  1:8792628/1‑101 (MQ=255)
  aGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATAc                                                                                               <  2:8792220/101‑1 (MQ=255)
   gATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATAcc                                                                                              >  2:8792633/1‑101 (MQ=255)
    aTAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCg                                                                                             <  2:8791993/101‑1 (MQ=255)
      aaaaCAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGt                                                                                           <  1:8791450/101‑1 (MQ=255)
      aaaaCAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGt                                                                                           >  1:8792639/1‑101 (MQ=255)
       aaaCAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTg                                                                                          >  2:8792640/1‑101 (MQ=255)
           aGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGaat                                                                                      <  1:8792017/101‑1 (MQ=255)
            gCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGaata                                                                                     >  2:8792648/1‑101 (MQ=255)
             cccGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGaataa                                                                                    <  1:8791819/101‑1 (MQ=255)
             cccGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGaataa                                                                                    >  1:8792652/1‑101 (MQ=255)
              ccGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGaataat                                                                                   <  2:8792296/101‑1 (MQ=255)
               cGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATg                                                                                  >  2:8792655/1‑101 (MQ=255)
               cGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATg                                                                                  <  2:8791957/101‑1 (MQ=255)
                gTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCgg                                                                                            >  1:8792659/1‑90 (MQ=255)
                gTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGc                                                                                 <  2:8792154/101‑1 (MQ=255)
                  cGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCg                                                                               >  1:8792664/1‑101 (MQ=255)
                   gACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGc                                                                              >  2:8792671/1‑101 (MQ=255)
                    aCCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCt                                                                             >  1:8792676/1‑101 (MQ=255)
                     ccccTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCt                                                                             >  2:8792679/1‑100 (MQ=255)
                      cccTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGaataat                                                                                   <  2:8791929/93‑1 (MQ=255)
                         tttATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAAtt                                                                        <  1:8791878/101‑1 (MQ=255)
                          ttattaATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAAttt                                                                       >  2:8792686/1‑101 (MQ=255)
                          ttattaATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAAttt                                                                       >  1:8792687/1‑101 (MQ=255)
                           tattaATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAAtttt                                                                      <  2:8792258/101‑1 (MQ=255)
                           tattaATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAAtttt                                                                      >  2:8792691/1‑101 (MQ=255)
                            attaATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAAttttt                                                                     <  2:8792037/101‑1 (MQ=255)
                              taATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTc                                                                   >  2:8792697/1‑101 (MQ=255)
                              taATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTc                                                                   <  1:8791992/101‑1 (MQ=255)
                               aaTCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCa                                                                  >  1:8792700/1‑101 (MQ=255)
                                aTCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCAt                                                                 <  2:8792016/101‑1 (MQ=255)
                                 tCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCAtt                                                                >  2:8792709/1‑101 (MQ=255)
                                 tCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCAtt                                                                >  1:8792708/1‑101 (MQ=255)
                                 tCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCAtt                                                                <  2:8791975/101‑1 (MQ=255)
                                    ggAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTAcc                                                             >  2:8792716/1‑101 (MQ=255)
                                     gAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCt                                                            >  1:8792719/1‑101 (MQ=255)
                                      aaGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTc                                                           <  1:8791649/101‑1 (MQ=255)
                                       aGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCa                                                          <  2:8792056/101‑1 (MQ=255)
                                       aGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCa                                                          <  2:8792113/101‑1 (MQ=255)
                                       aGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCa                                                          <  1:8791774/101‑1 (MQ=255)
                                        gTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAg                                                         >  2:8792725/1‑101 (MQ=255)
                                         tCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGt                                                        <  1:8791855/101‑1 (MQ=255)
                                         tCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGt                                                        <  1:8792295/101‑1 (MQ=255)
                                           ggCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATa                                                                                                <  1:8791932/59‑1 (MQ=255)
                                           ggCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGttt                                                      <  1:8792040/101‑1 (MQ=255)
                                           ggCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGttt                                                      >  1:8792729/1‑101 (MQ=255)
                                            gCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTc                                                     <  1:8792079/101‑1 (MQ=255)
                                            gCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTc                                                     >  2:8792732/1‑101 (MQ=255)
                                              aaCTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCt                                                   <  2:8792100/101‑1 (MQ=255)
                                              aaCTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCt                                                   >  1:8792737/1‑101 (MQ=255)
                                               aCTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCtt                                                  <  1:8792157/101‑1 (MQ=255)
                                                cTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCttt                                                 <  1:8792190/101‑1 (MQ=255)
                                                cTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCttt                                                 >  2:8792741/1‑101 (MQ=255)
                                                 ttCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCtttt                                                <  1:8792030/101‑1 (MQ=255)
                                                   cgcgCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCg                                              >  1:8792747/1‑101 (MQ=255)
                                                   cgcgCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCg                                              >  2:8792748/1‑101 (MQ=255)
                                                     cgcATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTa                                            <  1:8791978/101‑1 (MQ=255)
                                                     cgcATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTa                                            <  1:8792103/101‑1 (MQ=255)
                                                      gcATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACTTCAGTTTCCTTTTCGTaa                                           >  2:8792762/1‑101 (MQ=255)
                                                      gcATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTaa                                           >  1:8792761/1‑101 (MQ=255)
                                                      gcATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTaa                                           >  2:8792756/1‑101 (MQ=255)
                                                      gcATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTaa                                           <  2:8792338/101‑1 (MQ=255)
                                                       cATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAg                                          <  2:8792049/101‑1 (MQ=255)
                                                        aTTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGt                                         >  1:8792770/1‑101 (MQ=255)
                                                         ttGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTg                                        >  2:8792772/1‑101 (MQ=255)
                                                           gCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCt                                      >  2:8792779/1‑101 (MQ=255)
                                                             cTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCa                                    >  1:8792784/1‑101 (MQ=255)
                                                              tttttATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCAt                                   >  2:8792785/1‑101 (MQ=255)
                                                                tttATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCAttt                                 >  2:8792791/1‑101 (MQ=255)
                                                                 ttATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCAtttt                                >  1:8792793/1‑101 (MQ=255)
                                                                  tATCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTc                               >  2:8792802/1‑101 (MQ=255)
                                                                    tCAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGt                             >  1:8792806/1‑101 (MQ=255)
                                                                     cAAGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGtt                            <  2:8791904/101‑1 (MQ=255)
                                                                      aaGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGttt                           <  2:8792401/101‑1 (MQ=255)
                                                                      aaGGGTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGttt                           >  2:8792808/1‑101 (MQ=255)
                                                                       agtgTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTc                          >  1:8792813/4‑101 (MQ=255)
                                                                         ggTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTg                        >  2:8792816/1‑101 (MQ=255)
                                                                         ggTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTg                        <  2:8792077/101‑1 (MQ=255)
                                                                          gTTCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGt                       >  1:8792820/1‑101 (MQ=255)
                                                                           ttCCGTCACGATGAAACGTGCTGGATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTTATTTTCGTTTCTGTa                      <  1:8792225/101‑1 (MQ=255)
                                                                           ttCCGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTa                      <  1:8792063/101‑1 (MQ=255)
                                                                              cGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTAt                     <  1:8792199/99‑1 (MQ=255)
                                                                              cGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATcc                   >  2:8792828/1‑101 (MQ=255)
                                                                              cGTCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATcc                   >  1:8792829/1‑101 (MQ=255)
                                                                                tCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCa                 <  2:8792124/101‑1 (MQ=255)
                                                                                tCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCa                 <  1:8792345/101‑1 (MQ=255)
                                                                                tCACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCa                 <  1:8792252/101‑1 (MQ=255)
                                                                                 cACGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCAg                >  1:8792834/1‑101 (MQ=255)
                                                                                   cGATGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCAGcc              >  2:8792839/1‑101 (MQ=255)
                                                                                     aTGAAACGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCAGCCCa            >  1:8792841/1‑101 (MQ=255)
                                                                                        aaaCGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCAGCCCatat         <  1:8792127/101‑1 (MQ=255)
                                                                                         aaCGTGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCAGCCCATATg        <  2:8792182/101‑1 (MQ=255)
                                                                                             tGCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCAGCCCATATGatat    >  1:8792848/1‑101 (MQ=255)
                                                                                              gCTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCAGCCCATATGATAtt   <  1:8792414/101‑1 (MQ=255)
                                                                                               cTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCAGCCCATATGATAttt  >  1:8792854/1‑101 (MQ=255)
                                                                                               cTGAATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCAGCCCATATGATAttt  <  1:8792091/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
CCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGATGAAACGTGCTGGATACCGGTGGAATAATGCCGCTGAATTTTTTCATTACCTCAGTTTCCTTTTCGTAAGTGCTCATTTTCGTTTCTGTATCCCAGCCCATATGATATTT  >  NC_000913/4524878‑4525073

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑