Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,729,414 Δ1 bp 15.8% intergenic (‑43/+29) xylB ← / ← xylA xylulokinase/D‑xylose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,729,4130A.15.8% 542.5 / 96.3 184intergenic (‑42/+30)xylB/xylAxylulokinase/D‑xylose isomerase
Reads supporting (aligned to +/- strand):  ref base A (82/73);  new base . (20/9);  total (102/82)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 187 total aligned reads displayed.

TCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGA  >  NC_000913/3729316‑3729512
                                                                                                  |                                                                                                   
tCGTTGAGCAAAATAACTTTTACGACCGAGGTGCCAAGATCTATCCCGATATACATATCGATAAATACTAAAAAAAATGCCCGG‑TATCGCTACCGATa                                                                                                     <  1:3931117/98‑1 (MQ=255)
  gTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCg                                                                                                 >  1:869572/1‑100 (MQ=255)
    tGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCggg                                                                                               <  1:3011935/100‑1 (MQ=255)
    tGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCggg                                                                                               <  2:4357134/100‑1 (MQ=255)
     gAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGc                                                                                              >  2:939765/1‑100 (MQ=255)
     gAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGc                                                                                              >  2:1328974/1‑100 (MQ=255)
      aGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGcc                                                                                             >  1:2316351/1‑100 (MQ=255)
      aGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGcc                                                                                             >  1:3208697/1‑100 (MQ=255)
       gCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGGTTATCGGTAGCGAT‑ACCGGGcat                                                                                            <  2:3685349/100‑3 (MQ=37)
           aaTAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACgg                                                                                        <  1:1275612/100‑1 (MQ=255)
               aCTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACt                                                                                     >  1:665857/1‑99 (MQ=255)
                cTTTTACGCCCGAGGTGCCAAGATCTATCCCGCTATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTg                                                                                    <  2:1913970/99‑1 (MQ=255)
                cTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGc                                                                                   >  2:3928947/1‑100 (MQ=255)
                 ttttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGca                                                                                  <  1:189327/100‑1 (MQ=255)
                  tttACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGcac                                                                                 <  1:1839701/100‑1 (MQ=255)
                      cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                                             >  1:2588937/1‑100 (MQ=255)
                      cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                                             >  2:1143102/1‑100 (MQ=255)
                      cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                                             >  2:4011116/1‑100 (MQ=255)
                      cGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGtt                                                                             <  1:1933565/100‑1 (MQ=255)
                       gCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTa                                                                            <  1:2804504/100‑1 (MQ=255)
                         ccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAg                                                                           >  2:3604996/1‑99 (MQ=255)
                         ccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                                          <  1:4314844/100‑1 (MQ=255)
                         ccGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                                          >  1:1148236/1‑100 (MQ=255)
                          cGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGcc                                                                         <  2:1462020/100‑1 (MQ=255)
                          cGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGc                                                                          <  1:3054689/99‑1 (MQ=255)
                               tGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGtttt                                                                    <  2:4093656/100‑3 (MQ=255)
                               tGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAt                                                                    >  1:1907478/1‑100 (MQ=255)
                                gCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAtt                                                                   >  2:704978/1‑100 (MQ=255)
                                gCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAtt                                                                   >  2:258506/1‑100 (MQ=255)
                                gCCAAGATCTATCCCGATATACATATCGATCGTTACTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTAt                                                                    <  1:49239/99‑1 (MQ=255)
                                  cAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGTCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTg                                                                 >  2:320436/1‑100 (MQ=255)
                                    aGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTc                                                               <  2:2347115/100‑1 (MQ=255)
                                    aGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTc                                                               >  2:4162946/1‑100 (MQ=255)
                                    aGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACAGACTGCACAGTTAGCCGTTATTTGt                                                                >  1:4086801/1‑99 (MQ=255)
                                      aTCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGa                                                             <  1:4060153/100‑1 (MQ=255)
                                       tCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGaa                                                            <  2:4207955/100‑1 (MQ=255)
                                        cTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAAc                                                           <  1:156067/100‑1 (MQ=255)
                                           tCCCGATATACATATCGATCGTTCATTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGa                                                        >  1:2155736/1‑100 (MQ=255)
                                             ccGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATa                                                      >  2:2536930/1‑100 (MQ=255)
                                               gATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTTCACAGTTAGCCGTTATTTGTCGAACAGATAAt                                                    >  2:2857498/1‑100 (MQ=255)
                                                  ataCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGt                                                 <  1:2010562/100‑1 (MQ=255)
                                                   taCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGtt                                                >  1:789805/1‑100 (MQ=255)
                                                   taCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGATACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGtt                                                >  1:3522131/1‑100 (MQ=255)
                                                   naCATATCGATCGTTCCTTACAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCATTATTTGTCGAACAGATAATgg                                                  >  1:1624053/2‑98 (MQ=255)
                                                    aCATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGttt                                               >  1:1162824/1‑100 (MQ=255)
                                                     cATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTa                                              >  1:2505580/1‑100 (MQ=255)
                                                     cATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTa                                              >  1:3855714/1‑100 (MQ=255)
                                                      atatCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTAc                                             <  2:311203/100‑1 (MQ=255)
                                                       tatCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGATTccc                                            >  2:2060783/1‑100 (MQ=255)
                                                        atcgatcgTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCa                                           <  1:612262/100‑1 (MQ=255)
                                                          cgagcgTTCCTTAAAAAAATGCCCGG‑TATCGGTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTTGAACAGATAATGGTTTACCAGa                                         <  1:2558439/100‑1 (MQ=255)
                                                             tcgTTCCTTAAAAAAATGCCAGG‑TATCGATACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGAttt                                      >  1:529394/1‑100 (MQ=255)
                                                              cgTTCCATAAAAAAATGCCCGG‑TATCGCAACCGATAACCGGGCCAACGGACTGCACAGTTATCCGTTATTTGTCGAACAGATAATGGTTTACCAGAtttt                                     <  1:1387656/100‑1 (MQ=255)
                                                               gTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTc                                    >  2:3437896/1‑100 (MQ=255)
                                                                ttCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTcc                                   <  2:2613630/100‑1 (MQ=255)
                                                                ttCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTcc                                   <  1:3501053/100‑1 (MQ=255)
                                                                  ccTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAATGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAg                                 <  2:1614736/100‑1 (MQ=255)
                                                                  ccTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAg                                 >  1:630957/1‑100 (MQ=255)
                                                                   cTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgt                                >  1:2279886/1‑100 (MQ=255)
                                                                   cTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgt                                >  2:2666931/1‑100 (MQ=255)
                                                                   cTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgt                                >  2:3024253/1‑100 (MQ=255)
                                                                   cTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgt                                >  1:3664347/1‑100 (MQ=255)
                                                                    ttAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgtt                               >  1:3402665/1‑100 (MQ=255)
                                                                     tAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTTCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgttg                              <  1:2843335/100‑1 (MQ=255)
                                                                     tAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgttg                              >  2:1873925/1‑100 (MQ=255)
                                                                      aaaaaaaTGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAgttgt                             <  1:3164076/100‑1 (MQ=255)
                                                                        ccgttggCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTCGTCGAACAGCTACTGGTTTACCAGATTTTCCAGTTGTTc                           >  1:4338477/7‑100 (MQ=21)
                                                                            aTGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTgg                       <  1:3560413/100‑1 (MQ=255)
                                                                             ggCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGc                      >  2:3570659/2‑100 (MQ=38)
                                                                              gCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCg                     >  2:1261059/1‑100 (MQ=38)
                                                                              gCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCg                     >  1:3865094/1‑100 (MQ=38)
                                                                              gCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCg                     >  2:3502265/1‑100 (MQ=38)
                                                                              gCCCGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCg                     >  1:2997611/1‑100 (MQ=38)
                                                                               cccGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa                    >  2:826655/1‑100 (MQ=38)
                                                                               cccGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa                    >  2:898526/1‑100 (MQ=38)
                                                                               cccGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa                    >  2:2762663/1‑100 (MQ=38)
                                                                               cccGGTTATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGa                    >  2:2239613/1‑100 (MQ=38)
                                                                               cccGGTTATCGGTACCGAT‑ACCGGGCCAACGGACTGCACAGTTTGCCGTTATTTGTCGAACAGATAATGGCTTACCAGGTTTTCCACTTGTTACTGGCg                     >  1:3305594/1‑99 (MQ=255)
                                                                                ccggttATCGGTAGCGAT‑ACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAc                   >  2:2804526/6‑100 (MQ=255)
                                                                                ccGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAc                   >  2:110063/1‑100 (MQ=255)
                                                                                ccGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAc                   <  2:3141431/100‑1 (MQ=255)
                                                                                cagg‑TATCGCTACAGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAc                   >  2:2966454/3‑100 (MQ=255)
                                                                                 cGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGAcc                  >  2:3278728/1‑100 (MQ=255)
                                                                                   g‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCAc                <  1:2390101/100‑1 (MQ=255)
                                                                                       tCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCActct             >  2:3034767/1‑100 (MQ=255)
                                                                                        cGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCGGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTg            >  2:2448543/1‑100 (MQ=255)
                                                                                          cTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGAt          <  1:3819874/100‑1 (MQ=255)
                                                                                          cTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGAt          <  1:704977/100‑1 (MQ=255)
                                                                                           tACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATg         <  2:3738597/100‑1 (MQ=255)
                                                                                           tACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATg         <  1:2463129/100‑1 (MQ=255)
                                                                                            aCCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGc        <  2:3313401/100‑1 (MQ=255)
                                                                                            aCCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGc        <  1:3455741/100‑1 (MQ=255)
                                                                                            aCCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGc        <  1:3811380/100‑1 (MQ=255)
                                                                                            aCCGATAACCGGGCCAACGGACTACACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATg         >  1:3800695/1‑99 (MQ=255)
                                                                                             ccGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGAt          >  2:2740514/1‑97 (MQ=255)
                                                                                             ccGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCa       >  1:2478045/1‑100 (MQ=255)
                                                                                              cGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCAc      >  2:1189645/1‑100 (MQ=255)
                                                                                               gATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCAcc     >  2:105697/1‑100 (MQ=255)
                                                                                                aTAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCg    <  2:1785758/100‑1 (MQ=255)
                                                                                                 tAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCgg   <  2:2425011/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TCGTTGAGCAAAATAACTTTTACGCCCGAGGTGCCAAGATCTATCCCGATATACATATCGATCGTTCCTTAAAAAAATGCCCGG‑TATCGCTACCGATAACCGGGCCAACGGACTGCACAGTTAGCCGTTATTTGTCGAACAGATAATGGTTTACCAGATTTTCCAGTTGTTCCTGGCGACCACTCTGATGCACCGGA  >  NC_000913/3729316‑3729512

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