Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,419,837 A→G intergenic (‑212/‑143) ulaG ← / → ulaA L‑ascorbate‑6‑phosphate lactonase/L‑ascorbate specific PTS enzyme IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,419,8370AG99.8% 1845.3 / ‑5.7 431intergenic (‑212/‑143)ulaG/ulaAL‑ascorbate‑6‑phosphate lactonase/L‑ascorbate specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (215/215);  total (215/216)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 431 total aligned reads displayed.

AGTTGTTGTTGATTTTGAGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCACTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTAGAGTGAGTGCACAACATT  >  NC_000913/4419739‑4419934
                                                                                                  |                                                                                                 
aGTTGTTGTTGATTTTGAGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGct                                                                                                 <  1:8209080/101‑1 (MQ=255)
aGTTGTTGTTGATTTTGAGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGct                                                                                                 <  2:8208984/101‑1 (MQ=255)
aGTTGTTGTTGATTTTGAGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGct                                                                                                 >  2:8209511/1‑101 (MQ=255)
  ttgttgttgATTTTGAGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGctct                                                                                               >  1:8209516/1‑101 (MQ=255)
    gttgttgATTTTGAGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTAt                                                                                             <  1:8209305/101‑1 (MQ=255)
     ttgttgATTTTGAGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATc                                                                                            >  2:8209524/1‑101 (MQ=255)
     ttgttgATTTTGAGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATc                                                                                            >  1:8209522/1‑101 (MQ=255)
       gttgATTTTGAGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCGCGCTCTATCCa                                                                                          >  2:8209528/1‑101 (MQ=255)
         tgATTTTGAGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCAtg                                                                                        <  1:8209044/101‑1 (MQ=255)
          gATTTTGAGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCAtgt                                                                                       >  1:8209531/1‑101 (MQ=255)
              ttGAGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGaa                                                                                   >  2:8209537/1‑101 (MQ=255)
                 aGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAatta                                                                                <  1:8208531/101‑1 (MQ=255)
                  gTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAattat                                                                               >  2:8209543/1‑101 (MQ=255)
                   tATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAattatt                                                                              >  1:8209546/1‑101 (MQ=255)
                     tGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTg                                                                            <  2:8209310/101‑1 (MQ=255)
                     tGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTg                                                                            >  2:8209552/1‑101 (MQ=255)
                      ggAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGc                                                                           <  1:8208623/101‑1 (MQ=255)
                       gAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCg                                                                          >  1:8209554/1‑101 (MQ=255)
                        aaaGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCgg                                                                         >  2:8209556/1‑101 (MQ=255)
                        aaaGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCgg                                                                         <  2:8208695/101‑1 (MQ=255)
                          aGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGt                                                                       >  1:8209562/1‑101 (MQ=255)
                           gATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTc                                                                      >  2:8209566/1‑101 (MQ=255)
                           gATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTc                                                                      <  2:8208770/101‑1 (MQ=255)
                            aTTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTcg                                                                     >  1:8209570/1‑101 (MQ=255)
                             tttAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTcgc                                                                    >  1:8209575/1‑101 (MQ=255)
                                aaTGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTc                                                                 <  1:8209137/101‑1 (MQ=255)
                                aaTGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTc                                                                 >  2:8209582/1‑101 (MQ=255)
                                aaTGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTc                                                                 >  1:8209581/1‑101 (MQ=255)
                                 aTGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTca                                                                <  1:8208717/101‑1 (MQ=255)
                                 aTGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTca                                                                <  1:8208783/101‑1 (MQ=255)
                                   ggAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTcaca                                                              >  1:8209590/1‑101 (MQ=255)
                                   ggAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTcaca                                                              <  1:8209121/101‑1 (MQ=255)
                                    gAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACAt                                                             >  2:8209594/1‑101 (MQ=255)
                                    gAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACAt                                                             <  2:8208728/101‑1 (MQ=255)
                                     aaTGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACAtt                                                            <  2:8208902/101‑1 (MQ=255)
                                     aaTGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACAtt                                                            <  2:8208916/101‑1 (MQ=255)
                                      aTGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACAttt                                                           <  1:8208967/101‑1 (MQ=255)
                                       tgtgTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTa                                                          <  2:8209207/101‑1 (MQ=255)
                                       tgtgTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTa                                                          <  2:8208935/101‑1 (MQ=255)
                                       tgtgTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTa                                                          >  1:8209598/1‑101 (MQ=255)
                                        gtgtAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTaa                                                         >  2:8209601/1‑101 (MQ=255)
                                         tgtAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAAt                                                        <  2:8208826/101‑1 (MQ=255)
                                         tgtAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAAt                                                        <  2:8208994/101‑1 (MQ=255)
                                          gtAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATc                                                       >  1:8209607/1‑101 (MQ=255)
                                            aaTTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCat                                                     <  2:8208643/101‑1 (MQ=255)
                                             aTTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCata                                                    >  2:8209611/1‑101 (MQ=255)
                                               tCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCataaa                                                  >  1:8209615/1‑101 (MQ=255)
                                                cAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCataaat                                                 <  1:8208795/101‑1 (MQ=255)
                                                 aaTTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCataaata                                                <  1:8208319/101‑1 (MQ=255)
                                                 aaTTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCataaata                                                >  2:8209620/1‑101 (MQ=255)
                                                 aaTTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCataaata                                                <  1:8208917/101‑1 (MQ=255)
                                                  aTTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCataaataa                                               >  1:8209622/1‑101 (MQ=255)
                                                   ttAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCGCATTTAATCATAaataat                                              >  1:8209626/1‑101 (MQ=255)
                                                     aaCTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACAttt                                                           >  2:8209631/1‑86 (MQ=255)
                                                        tAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCttgt                                         >  2:8209639/1‑101 (MQ=255)
                                                        tAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCttgt                                         >  1:8209640/1‑101 (MQ=255)
                                                        tAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCttgt                                         <  1:8208950/101‑1 (MQ=255)
                                                         aaaTGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCttgtt                                        >  2:8209647/1‑101 (MQ=255)
                                                         aaaTGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCttgtt                                        >  1:8209651/1‑101 (MQ=255)
                                                         aaaTGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCttgtt                                        <  2:8208976/101‑1 (MQ=255)
                                                         aaaTGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCttgtt                                        <  2:8208944/101‑1 (MQ=255)
                                                          aaTGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCttgttg                                       >  2:8209652/1‑101 (MQ=255)
                                                           aTGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTtgttgt                                      <  1:8208878/101‑1 (MQ=255)
                                                            tGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTtgtg                                     <  2:8209053/101‑1 (MQ=255)
                                                              aaTTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCat                                                     >  1:8209660/1‑83 (MQ=255)
                                                              aaTTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGAt                                   <  2:8209236/101‑1 (MQ=255)
                                                               aTTTAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGAtt                                  >  2:8209661/1‑101 (MQ=255)
                                                                tttAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTa                                 <  1:8209175/101‑1 (MQ=255)
                                                                 ttAAATGGATAATTGTTTCGTTGTGTGANTCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTAc                                >  2:8209668/1‑101 (MQ=255)
                                                                 ttAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTAc                                <  1:8209001/101‑1 (MQ=255)
                                                                 ttAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTAc                                <  1:8208929/101‑1 (MQ=255)
                                                                  tAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACt                               >  2:8209672/1‑101 (MQ=255)
                                                                  tAAATGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACt                               >  1:8209674/1‑101 (MQ=255)
                                                                   aaaTGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACtt                              >  2:8209679/1‑101 (MQ=255)
                                                                      tGGATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTg                           >  1:8209685/1‑101 (MQ=255)
                                                                       ggATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGa                          <  1:8208889/101‑1 (MQ=255)
                                                                       ggATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGa                          >  1:8209689/1‑101 (MQ=255)
                                                                        gATAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGaa                         <  2:8209086/101‑1 (MQ=255)
                                                                         aTAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGaaa                        <  2:8209114/101‑1 (MQ=255)
                                                                         aTAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGaaa                        >  2:8209691/1‑101 (MQ=255)
                                                                          tAATTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGaaaa                       <  1:8208823/101‑1 (MQ=255)
                                                                           aaTTGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAAt                      <  1:8209212/101‑1 (MQ=255)
                                                                             ttGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTa                    >  2:8209698/1‑101 (MQ=255)
                                                                              tGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTag                   <  2:8208476/101‑1 (MQ=255)
                                                                              tGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTag                   >  2:8209703/1‑101 (MQ=255)
                                                                              tGTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTag                   >  1:8209705/1‑101 (MQ=255)
                                                                               gTTTCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCttgt                                         <  1:8209426/78‑1 (MQ=255)
                                                                                tttCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTagag                 >  2:8209709/1‑101 (MQ=255)
                                                                                tttCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTagag                 <  2:8209149/101‑1 (MQ=255)
                                                                                tttCGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTagag                 <  2:8208872/101‑1 (MQ=255)
                                                                                   cGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGaaaa                       >  2:8209717/1‑92 (MQ=255)
                                                                                   cGTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTAgagtga              >  1:8209716/1‑101 (MQ=255)
                                                                                    gTTGTGTGAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTAgagtgag             <  1:8208853/101‑1 (MQ=255)
                                                                                        tgtgAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTAGAGTGAGTGca         <  2:8209266/101‑1 (MQ=255)
                                                                                         gtgAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTAGAGTGAGTGcac        >  1:8209723/1‑101 (MQ=255)
                                                                                           gAATCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTAGAGTGAGTGCacaa      <  1:8209157/101‑1 (MQ=255)
                                                                                              tCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTAGAGTGAGTGCACAACAt   <  1:8209247/101‑1 (MQ=255)
                                                                                              tCCCGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTAGAGTGAGTGCACAACAt   <  1:8209068/101‑1 (MQ=255)
                                                                                               cccGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTAGAGTGAGTGCACAACAtt  <  2:8209028/101‑1 (MQ=255)
                                                                                               cccGCTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTAGAGTGAGTGCACAACAtt  <  1:8209029/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
AGTTGTTGTTGATTTTGAGTATGGAAAGATTTAATGGAATGTGTAATTCAATTAACTAAATGAATTTAAATGGATAATTGTTTCGTTGTGTGAATCCCACTCTATCCATGTGGAATTATTTGCGGGTCGCGTCACATTTAATCATAAATAATCTTGTTGTGATTACTTTTGAAAATTAGAGTGAGTGCACAACATT  >  NC_000913/4419739‑4419934

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: ‑