Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,911,433 A→G intergenic (‑137/+25) glmS ← / ← glmU L‑glutamine:D‑fructose‑6‑phosphate aminotransferase/N‑acetyl glucosamine‑1‑phosphate uridyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,911,4330AG84.2% 644.3 / 117.0 284intergenic (‑137/+25)glmS/glmUL‑glutamine:D‑fructose‑6‑phosphate aminotransferase/N‑acetyl glucosamine‑1‑phosphate uridyltransferase
Reads supporting (aligned to +/- strand):  ref base A (24/21);  new base G (122/117);  total (146/138)
Fisher's exact test for biased strand distribution p-value = 8.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.91e-02

Only 100 of 286 total aligned reads displayed.

GTTTTTTCGGGCGCCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCACACGGCTGATAGCTAAT  >  NC_000913/3911343‑3911524
                                                                                          |                                                                                           
gTTTTTTCGGGCGCCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGtg                                                                                           >  2:3793796/1‑93 (MQ=255)
 ttttttCGGGCGCCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCgtgt                                                                                          <  1:7317527/93‑1 (MQ=255)
 ttttttCGGGCGCCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCAtgt                                                                                          >  2:4541531/1‑93 (MQ=255)
  tttttCGGGCGCCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGtgtg                                                                                         >  2:6842217/1‑93 (MQ=255)
   ttttCGGGCGCCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGa                                                                                        >  1:6587720/1‑93 (MQ=255)
     ttCGGGCGCCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGAcc                                                                                      <  2:2642450/93‑1 (MQ=255)
      tCGGGCGCCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCg                                                                                     >  2:7158363/1‑93 (MQ=255)
        gggCGCCCCGAGCCTTGTAGAGAGTGGAGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCggg                                                                                   <  1:4920345/93‑1 (MQ=255)
         ggCGCCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGt                                                                                  <  1:2038059/93‑1 (MQ=255)
          gcgcCCCGAGCTTGGTAGAGATTGGGTTTATTTTTATAGTTACTGCTTGTGGGAGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTa                                                                                <  2:1226204/93‑1 (MQ=255)
          gcgcCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGtt                                                                                 >  2:1841972/1‑93 (MQ=255)
           cgcCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTa                                                                                <  1:6939639/93‑1 (MQ=255)
           cgcCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTa                                                                                <  2:1419573/93‑1 (MQ=255)
            gcCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCACGTGACCGGGTTAg                                                                               >  1:3294739/1‑93 (MQ=255)
             ccccGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGc                                                                              <  2:2928525/93‑1 (MQ=255)
                 gAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGc                                                                          <  2:6829989/93‑1 (MQ=255)
                  aGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGcc                                                                         >  1:3922971/1‑93 (MQ=255)
                  aGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGcc                                                                         >  2:2734859/1‑93 (MQ=255)
                   gCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCa                                                                        >  1:4986883/1‑93 (MQ=255)
                   gCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCa                                                                        >  1:3747626/1‑93 (MQ=255)
                     cTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGa                                                                      <  2:4860784/93‑1 (MQ=255)
                      ttGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGaa                                                                     >  1:6995380/1‑93 (MQ=255)
                      ttGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGaa                                                                     >  1:152944/1‑93 (MQ=255)
                      ttGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGaa                                                                     >  1:71531/1‑93 (MQ=255)
                       tGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAAt                                                                    <  1:6236240/93‑1 (MQ=255)
                       tGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAAt                                                                    <  2:4258134/93‑1 (MQ=255)
                         tAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCa                                                                  <  2:3998500/93‑1 (MQ=255)
                         tAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCa                                                                  <  1:4363780/93‑1 (MQ=255)
                          agagagTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCAc                                                                 <  1:2933166/93‑1 (MQ=255)
                           gagagTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACt                                                                >  1:471068/1‑93 (MQ=255)
                            agagTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACtt                                                               <  1:106783/93‑1 (MQ=255)
                             gagTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACttt                                                              >  1:5820358/1‑93 (MQ=255)
                             gagTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACttt                                                              <  1:5419437/93‑1 (MQ=255)
                              agTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACtttt                                                             >  2:356344/1‑93 (MQ=255)
                              agTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACtttt                                                             >  1:2612821/1‑93 (MQ=255)
                              agTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACtttt                                                             <  1:3022930/93‑1 (MQ=255)
                                tGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTc                                                           <  2:4098386/93‑1 (MQ=255)
                                 ggggTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCt                                                          >  1:1932354/1‑93 (MQ=255)
                                  gggTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCtt                                                         >  1:3167335/1‑93 (MQ=255)
                                  gggTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCtt                                                         <  1:646317/93‑1 (MQ=255)
                                  gggTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCtt                                                         <  2:119629/93‑1 (MQ=255)
                                     ttATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTAc                                                      <  2:3428538/93‑1 (MQ=255)
                                     ttATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTAc                                                      >  2:5795728/1‑93 (MQ=255)
                                      tATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTAcc                                                     >  2:7255260/1‑93 (MQ=255)
                                      tATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTAcc                                                     >  1:2799552/1‑93 (MQ=255)
                                       aTTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCg                                                    >  2:2224285/1‑93 (MQ=255)
                                         ttttATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGga                                                  >  2:6554195/1‑93 (MQ=255)
                                          tttATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGgac                                                 <  2:1783015/93‑1 (MQ=255)
                                          tttATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGgac                                                 <  1:5622602/93‑1 (MQ=255)
                                            tataGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCcgacga                                               <  2:4595527/93‑1 (MQ=255)
                                             ataGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGgacgac                                              <  2:5579325/93‑1 (MQ=255)
                                             ataGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGgacgac                                              <  2:5161062/93‑1 (MQ=255)
                                              taGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGacgacg                                             >  2:4956366/1‑93 (MQ=255)
                                              taGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGacgacg                                             >  1:5081433/1‑93 (MQ=255)
                                               aGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGc                                            >  1:3444797/1‑93 (MQ=255)
                                                gTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGcc                                           <  2:6130971/93‑1 (MQ=255)
                                                gTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGcc                                           <  1:7035537/93‑1 (MQ=255)
                                                 ttACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCa                                          <  1:847374/93‑1 (MQ=255)
                                                 ttACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCa                                          >  1:6303189/1‑93 (MQ=255)
                                                   aCTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGc                                        <  1:444323/93‑1 (MQ=255)
                                                   aCTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGc                                        >  2:147581/1‑93 (MQ=255)
                                                    cTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGcc                                       <  2:1141288/93‑1 (MQ=255)
                                                     tGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCct                                      >  1:6351303/1‑93 (MQ=255)
                                                       cTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCcttc                                    <  1:2185725/93‑1 (MQ=255)
                                                       cTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCcttc                                    <  1:2458561/93‑1 (MQ=255)
                                                       cTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCcttc                                    <  1:6755076/93‑1 (MQ=255)
                                                        ttGTGGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGTTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCcttct                                   >  2:6903596/1‑93 (MQ=255)
                                                         tgtgGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCcttctt                                  <  1:3112815/93‑1 (MQ=255)
                                                         tgtgGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCcttctt                                  <  2:2859874/93‑1 (MQ=255)
                                                          gtgGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCttcttt                                 >  1:5229395/1‑93 (MQ=255)
                                                          gtgGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCttcttt                                 >  2:6486288/1‑93 (MQ=255)
                                                           tgGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCttctttc                                <  1:5034745/93‑1 (MQ=255)
                                                           tgGGAGGGAGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCttctttc                                <  1:1965263/93‑1 (MQ=255)
                                                            gggagggaGATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCttctttct                               >  1:1153692/1‑93 (MQ=255)
                                                             ggagggagATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTg                              >  2:636427/1‑93 (MQ=255)
                                                                gggagATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGt                           >  2:5872282/1‑93 (MQ=255)
                                                                gggagATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGt                           >  1:7479799/1‑93 (MQ=255)
                                                                 ggagATTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTc                          >  2:5675145/1‑93 (MQ=255)
                                                                 ggagATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTc                          >  2:5995624/1‑93 (MQ=255)
                                                                    gaTTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGc                       <  1:1497624/93‑1 (MQ=255)
                                                                    gaTTATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGc                       >  1:6094187/1‑93 (MQ=255)
                                                                      ttATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCgg                     >  2:7494296/1‑93 (MQ=255)
                                                                      ttATGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCgg                     >  2:4711614/1‑93 (MQ=255)
                                                                        aTGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGca                   <  1:364996/93‑1 (MQ=255)
                                                                         tGTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGcac                  >  1:170895/1‑93 (MQ=255)
                                                                         tGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGcac                  >  2:5237186/1‑93 (MQ=255)
                                                                          gTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGcaca                 <  1:451849/93‑1 (MQ=255)
                                                                          gTTATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGcaca                 <  1:6842217/93‑1 (MQ=255)
                                                                           ttATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCacac                >  1:6126513/1‑93 (MQ=255)
                                                                           ttATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCacac                <  1:4023847/93‑1 (MQ=255)
                                                                           ttATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCacac                <  2:3268940/93‑1 (MQ=255)
                                                                           ttATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCacac                >  2:4317498/1‑93 (MQ=255)
                                                                            tATCTCCTCATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCACACg               <  2:6707162/93‑1 (MQ=255)
                                                                                 cctcATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCACACGGCTGa          >  1:7516107/1‑93 (MQ=255)
                                                                                 cctcATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCACACGGCTGa          <  2:1359793/93‑1 (MQ=255)
                                                                                  ctcATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCACACGGCTGAt         >  1:4083462/1‑93 (MQ=255)
                                                                                  ctcATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCACACGGCTGAt         <  2:7311315/93‑1 (MQ=255)
                                                                                   tcATCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCACACGGCTGATa        <  2:2503750/93‑1 (MQ=255)
                                                                                     aTCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCACACGGCTGATAGc      >  2:5558172/1‑93 (MQ=255)
                                                                                      tCCCGTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCACACGGCTGATAGCt     <  1:5398681/93‑1 (MQ=255)
                                                                                          |                                                                                           
GTTTTTTCGGGCGCCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCACACGGCTGATAGCTAAT  >  NC_000913/3911343‑3911524

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