Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,912,152 A→G V226A (GTG→GCG)  glmU ← N‑acetyl glucosamine‑1‑phosphate uridyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,912,1520AG86.6% 787.5 / 109.5 329V226A (GTG→GCG) glmUN‑acetyl glucosamine‑1‑phosphate uridyltransferase
Reads supporting (aligned to +/- strand):  ref base A (19/25);  new base G (132/152);  total (152/177)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

Only 100 of 337 total aligned reads displayed.

GATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTAT  >  NC_000913/3912060‑3912244
                                                                                            |                                                                                            
  tCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACg                                                                                            <  2:3568620/93‑1 (MQ=255)
   cgcgCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCgcgc                                                                                           <  2:1123563/93‑1 (MQ=255)
     cgcAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCt                                                                                         >  2:4984807/1‑93 (MQ=255)
      gcagcaTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCtt                                                                                        >  2:3423973/1‑93 (MQ=255)
       cagcaTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTc                                                                                       <  2:7352136/93‑1 (MQ=255)
       cagcaTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTc                                                                                       <  2:142373/93‑1 (MQ=255)
        agcaTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCt                                                                                      >  2:6331912/1‑93 (MQ=255)
         gcaTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTa                                                                                     <  2:824228/93‑1 (MQ=255)
         gcaTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTa                                                                                     >  1:4294991/1‑93 (MQ=255)
         gcaTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTa                                                                                     <  2:6075111/93‑1 (MQ=255)
           aTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACt                                                                                   >  2:1655248/1‑93 (MQ=255)
           aTAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACt                                                                                   <  1:5821915/93‑1 (MQ=255)
            tAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACtt                                                                                  <  2:6723540/93‑1 (MQ=255)
            tAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACtt                                                                                  >  2:2010959/1‑93 (MQ=255)
             aaCGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTc                                                                                 >  2:4392213/1‑93 (MQ=255)
              aCGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCg                                                                                <  2:7259743/93‑1 (MQ=255)
                gCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCt                                                                              <  1:318261/93‑1 (MQ=255)
                  cTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTa                                                                            >  1:802795/1‑93 (MQ=255)
                   tGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTaa                                                                           >  1:1614746/1‑93 (MQ=255)
                    gCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTaaa                                                                          <  1:4881382/93‑1 (MQ=255)
                    gCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTaaa                                                                          <  1:4223764/93‑1 (MQ=255)
                      tAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACg                                                                        <  2:911612/93‑1 (MQ=255)
                        aCAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGtt                                                                      <  1:4416052/93‑1 (MQ=255)
                        aCAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGtt                                                                      >  2:1758363/1‑93 (MQ=255)
                        aCAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGtt                                                                      >  1:5124091/1‑93 (MQ=255)
                         cagcagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTg                                                                     >  2:5102134/1‑93 (MQ=255)
                         cagcagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTg                                                                     >  2:4472393/1‑93 (MQ=255)
                          agcagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGc                                                                    >  1:2561483/1‑93 (MQ=255)
                          agcagTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGc                                                                    <  1:247392/93‑1 (MQ=255)
                             agTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGa                                                                 >  2:3912860/1‑93 (MQ=255)
                              gTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGAt                                                                >  1:1317522/1‑93 (MQ=255)
                              gTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGAt                                                                <  1:6195147/93‑1 (MQ=255)
                              gTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGAt                                                                >  2:5555051/1‑93 (MQ=255)
                               tttttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATg                                                               >  1:3997201/1‑93 (MQ=255)
                               tttttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATg                                                               <  2:1283576/93‑1 (MQ=255)
                                ttttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGa                                                              >  1:2935662/1‑93 (MQ=255)
                                ttttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGa                                                              <  2:3484885/93‑1 (MQ=255)
                                 tttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGaa                                                             <  2:3395765/93‑1 (MQ=255)
                                 tttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGaa                                                             <  1:2969726/93‑1 (MQ=255)
                                 tttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGaa                                                             <  2:5921169/93‑1 (MQ=255)
                                  ttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAAc                                                            <  2:2163347/93‑1 (MQ=255)
                                  ttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAAc                                                            >  1:3107857/1‑93 (MQ=255)
                                  ttCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCCCCTTCTACTTCGCTTAAACGTTGCGGATGAAc                                                            <  1:7425241/93‑1 (MQ=255)
                                   tCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACg                                                           <  1:5944301/93‑1 (MQ=255)
                                   tCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACg                                                           <  2:3768908/93‑1 (MQ=255)
                                     aGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGc                                                         <  1:2213414/93‑1 (MQ=255)
                                      gCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGcg                                                        >  2:659121/1‑93 (MQ=255)
                                       ccTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGcga                                                       <  1:4683613/93‑1 (MQ=255)
                                       ccTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGcga                                                       <  1:3534340/93‑1 (MQ=255)
                                        cTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGcgac                                                      <  2:5790878/93‑1 (MQ=255)
                                         tGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGcgacg                                                     >  1:5491871/1‑93 (MQ=255)
                                         tGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGGTTTTTCCCCCTTCTACTTCCCTTAAACGTTGCGGGTGAACGGcgacg                                                     >  1:5281711/1‑93 (MQ=255)
                                          gTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGcgacga                                                    <  1:1114769/93‑1 (MQ=255)
                                           ttCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGAt                                                   >  1:694615/1‑93 (MQ=255)
                                           ttCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGAt                                                   <  1:5177845/93‑1 (MQ=255)
                                            tCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGAtt                                                  <  2:4159226/93‑1 (MQ=255)
                                              ggACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTc                                                >  1:7162994/1‑93 (MQ=255)
                                               gACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCa                                               <  1:6705628/93‑1 (MQ=255)
                                                     tAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCt                                         >  1:4765295/1‑93 (MQ=255)
                                                       aaCACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTTAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTc                                       <  1:43594/93‑1 (MQ=255)
                                                       aaCACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTc                                       >  1:988103/1‑93 (MQ=255)
                                                       aaCACGCTCCAGACGGGAGAGCTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGTTGAACGGCGTCGATTTCACGCCCTTc                                       >  1:7396007/1‑93 (MQ=255)
                                                        acacGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTcc                                      <  1:1848583/93‑1 (MQ=255)
                                                        acacGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTcc                                      <  1:1550733/93‑1 (MQ=255)
                                                         cacGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCt                                     <  1:1682679/93‑1 (MQ=255)
                                                         cacGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCt                                     <  2:6512785/93‑1 (MQ=255)
                                                          acGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTg                                    <  2:2657969/93‑1 (MQ=255)
                                                          acGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTg                                    <  1:6061726/93‑1 (MQ=255)
                                                           cGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGa                                   <  2:5674317/93‑1 (MQ=255)
                                                           cGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGa                                   <  2:278499/93‑1 (MQ=255)
                                                           cGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGa                                   <  2:5699609/93‑1 (MQ=255)
                                                           cGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGa                                   <  1:1398780/93‑1 (MQ=255)
                                                            gCTCCAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGAt                                  >  1:5413335/1‑93 (MQ=255)
                                                               ccAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACg                               >  1:5609048/1‑93 (MQ=255)
                                                                cAGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGc                              <  2:1471688/93‑1 (MQ=255)
                                                                 aGACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGcc                             <  1:2769619/93‑1 (MQ=255)
                                                                  gACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCa                            <  2:7212728/93‑1 (MQ=255)
                                                                  gACGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCa                            >  2:2886053/1‑93 (MQ=255)
                                                                    cGGGAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAgc                          >  2:540440/1‑93 (MQ=255)
                                                                      ggAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAgcgc                        <  2:4233162/93‑1 (MQ=255)
                                                                      ggAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAgcgc                        >  2:3284953/1‑93 (MQ=255)
                                                                      ggAGAGTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAgcgc                        >  2:3057440/1‑93 (MQ=255)
                                                                         gagTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCaat                     >  2:2526408/1‑93 (MQ=255)
                                                                          agTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCaata                    >  2:5143982/1‑93 (MQ=255)
                                                                           gTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCaataa                   >  1:4033554/1‑93 (MQ=255)
                                                                           gTTGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCaataa                   <  1:270299/93‑1 (MQ=255)
                                                                             tGCAGGCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATg                 >  1:3605160/1‑93 (MQ=255)
                                                                               cAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTc               <  2:6773378/93‑1 (MQ=255)
                                                                                 ggCGGTTATTCGCTCCTTCTACTTCGCTTAAGCGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCgg             <  1:6451196/93‑1 (MQ=255)
                                                                                 ggCGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCgg             <  2:2956173/93‑1 (MQ=255)
                                                                                   cGGTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTg           >  2:1865819/1‑93 (MQ=255)
                                                                                    ggTTATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGa          <  1:2136696/93‑1 (MQ=255)
                                                                                    ggTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGa          <  2:3993203/93‑1 (MQ=255)
                                                                                      ttATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATg        <  2:4579549/93‑1 (MQ=255)
                                                                                      ttATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATg        <  1:4585665/93‑1 (MQ=255)
                                                                                       tATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATgt       <  1:4990165/93‑1 (MQ=255)
                                                                                       tATTCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATgt       <  1:3941306/93‑1 (MQ=255)
                                                                                         ttCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATgtag     >  2:3651067/1‑93 (MQ=255)
                                                                                         ttCGCGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATgtag     >  1:6846020/1‑93 (MQ=255)
                                                                                           cACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATgtagta   <  2:6884912/93‑1 (MQ=255)
                                                                                            |                                                                                            
GATCGCGCAGCATAACGCCTGCTAACAGCAGTTTTTCAGCCTGTTCGGACTGATAAACACGCTCCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTAT  >  NC_000913/3912060‑3912244

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