Predicted mutation
evidence seq id position mutation annotation gene description
RA AC_000091 1,132,560 1 bp→TT intergenic (‑120/‑35) flgA ← / → flgB assembly protein for flagellar basal‑body periplasmic P ring/flagellar component of cell‑proximal portion of basal‑body rod

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AC_0000911,132,5591.T96.2% 411.5 / 9.1 105intergenic (‑119/‑36)flgA/flgBassembly protein for flagellar basal‑body periplasmic P ring/flagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base T (64/37);  total (66/39)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-02
*AC_0000911,132,5600GT100.0% 428.1 / NA 105intergenic (‑120/‑35)flgA/flgBassembly protein for flagellar basal‑body periplasmic P ring/flagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (66/39);  total (66/39)

Only 100 of 114 total aligned reads displayed.

GCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCG  >  AC_000091/1132483‑1132650
                                                                              |||                                                                                         
gCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGa                                                                                 >  2:2290859/1‑90 (MQ=255)
gCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGa                                                                                 >  2:407344/1‑90 (MQ=255)
 cAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGaa                                                                                >  1:1021629/1‑90 (MQ=255)
  aaTGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAAt                                                                               >  1:1517904/1‑90 (MQ=255)
  aaTGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAAt                                                                               <  1:692877/90‑1 (MQ=255)
   aTGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAAtt                                                                              >  1:3140192/1‑90 (MQ=255)
    tGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAAttt                                                                             >  2:2018367/1‑90 (MQ=255)
      aCCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTg                                                                           >  2:465062/1‑90 (MQ=255)
      aCCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTg                                                                           >  1:1561301/1‑90 (MQ=255)
       ccATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGa                                                                          >  2:67963/1‑90 (MQ=255)
       ccATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGa                                                                          <  2:1131840/90‑1 (MQ=255)
       ccATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGa                                                                          >  2:1451291/1‑90 (MQ=255)
        cATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGAt                                                                         >  1:1524271/1‑90 (MQ=255)
         aTCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATa                                                                        <  1:2548616/90‑1 (MQ=255)
          tCAACGGCATAAAGAGCGACCCATTTCGCGTTTATTCCGCGGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATAc                                                                       <  2:1972452/90‑1 (MQ=255)
           cAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATAcc                                                                      >  1:1227028/1‑90 (MQ=255)
           cAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATAcc                                                                      >  2:1290224/1‑90 (MQ=255)
           cAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATAcc                                                                      >  1:2580463/1‑90 (MQ=255)
            aaCGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCt                                                                     >  2:1020824/1‑90 (MQ=255)
            aaCGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCt                                                                     >  2:2450936/1‑90 (MQ=255)
             aCGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTg                                                                    <  1:2828383/90‑1 (MQ=255)
              cGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGc                                                                   >  1:17128/1‑90 (MQ=255)
               ggCATAAATAGCGCCCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCg                                                                  <  2:660350/90‑1 (MQ=255)
               ggCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCg                                                                  >  2:2762570/1‑90 (MQ=255)
                 cATAAATAGCGCCCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCgga                                                                <  2:361219/90‑1 (MQ=255)
                 cATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCgga                                                                >  1:1043893/1‑90 (MQ=255)
                   tAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCggagg                                                              >  2:2920973/1‑90 (MQ=255)
                   tAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCggagg                                                              <  2:102622/90‑1 (MQ=255)
                    aaaTAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCggagga                                                             >  2:825936/1‑90 (MQ=255)
                     aaTAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGgaggag                                                            >  2:869530/1‑90 (MQ=255)
                      aTAGTGACCTATTTTTCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGgaga                                                           <  2:1102740/90‑1 (MQ=255)
                      aTAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGgaga                                                           >  1:2524552/1‑90 (MQ=255)
                      aTAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGgaga                                                           <  2:3173992/90‑1 (MQ=255)
                       tAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGat                                                          >  1:3103953/1‑90 (MQ=255)
                       tAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGat                                                          <  1:33941/90‑1 (MQ=255)
                         gcgcCCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGatat                                                        <  2:948152/90‑1 (MQ=255)
                          cGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATg                                                       >  2:2779141/1‑90 (MQ=255)
                          cGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATg                                                       <  1:2759405/90‑1 (MQ=255)
                           gACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGc                                                      >  1:1665574/1‑90 (MQ=255)
                           gACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGc                                                      <  2:2919519/90‑1 (MQ=255)
                           gACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGc                                                      <  1:208373/90‑1 (MQ=255)
                            aCCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCt                                                     <  1:346161/90‑1 (MQ=255)
                            aCCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCt                                                     <  2:1231665/90‑1 (MQ=255)
                            aCCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCt                                                     <  1:1913531/90‑1 (MQ=255)
                              ccATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCg                                                   >  2:227848/1‑90 (MQ=255)
                                aTTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGTATTTTGATACCTGCGGAGGAGATATGCTCGAt                                                 >  2:2394769/1‑90 (MQ=255)
                                aTTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGAt                                                 <  2:748750/90‑1 (MQ=255)
                                aTTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGAt                                                 >  1:295292/1‑90 (MQ=255)
                                aTTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGAt                                                 >  1:1174424/1‑90 (MQ=255)
                                aTTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGAt                                                 <  1:3112669/90‑1 (MQ=255)
                                aTTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGAt                                                 >  1:218672/1‑90 (MQ=255)
                                  tttGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATaa                                               <  1:493745/90‑1 (MQ=255)
                                  tttGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATaa                                               <  1:1210431/90‑1 (MQ=255)
                                   ttGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAg                                              >  1:2535260/1‑90 (MQ=255)
                                   ttGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAg                                              >  1:1509305/1‑90 (MQ=255)
                                    tGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGc                                             >  2:2240309/1‑90 (MQ=255)
                                     gCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCt                                            <  1:1284210/90‑1 (MQ=255)
                                      cGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTc                                           <  2:789422/90‑1 (MQ=255)
                                       gTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCg                                          >  2:3062367/1‑90 (MQ=255)
                                       gTTTATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCg                                          >  2:914294/1‑90 (MQ=255)
                                         ttATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGAc                                        <  1:501120/90‑1 (MQ=255)
                                          tATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGAcg                                       >  2:2238248/1‑90 (MQ=255)
                                          tATTCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGAcg                                       <  2:2093064/90‑1 (MQ=255)
                                            ttCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGAcgcc                                     >  1:2083364/1‑90 (MQ=255)
                                            ttCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGAcgcc                                     <  2:104551/90‑1 (MQ=255)
                                            ttCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGAcgcc                                     >  1:1132933/1‑90 (MQ=255)
                                            ttCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGAcgcc                                     >  2:845820/1‑90 (MQ=255)
                                             tCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGAcgccg                                    <  1:630569/90‑1 (MQ=255)
                                             tCCGCCGATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGAcgccg                                    >  1:3214877/1‑90 (MQ=255)
                                               cgccgATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACgccgcc                                  <  1:457649/90‑1 (MQ=255)
                                               cgccgATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACgccgcc                                  >  2:2338746/1‑90 (MQ=255)
                                                 ccgATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCtt                                >  1:2575764/1‑90 (MQ=255)
                                                  cgATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGANTTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCTGCCTTa                               <  2:1244767/90‑1 (MQ=255)
                                                  cgATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTa                               >  2:2068522/1‑90 (MQ=255)
                                                  cgATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTa                               >  1:2898262/1‑90 (MQ=255)
                                                  cgATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTa                               >  1:3221805/1‑90 (MQ=255)
                                                  cgATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTa                               >  1:473256/1‑90 (MQ=255)
                                                   gATAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTAc                              >  1:927349/1‑90 (MQ=255)
                                                    aTAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACg                             <  2:2546401/90‑1 (MQ=255)
                                                    aTAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACg                             >  1:483558/1‑90 (MQ=255)
                                                    aTAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACg                             <  2:2908180/90‑1 (MQ=255)
                                                     tAACGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGt                            >  1:1419782/1‑90 (MQ=255)
                                                      aaCGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGtt                           >  2:1721471/1‑90 (MQ=255)
                                                       aCGCGCGCGTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGttt                          >  2:3316088/1‑90 (MQ=255)
                                                        cgcgcgcgTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGtttt                         >  2:1539748/1‑90 (MQ=255)
                                                        cgcgcgcgTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGtttt                         <  2:3159655/90‑1 (MQ=255)
                                                        cgcgcgcgTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGtttt                         >  1:1119595/1‑90 (MQ=255)
                                                          cgcgcgTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTca                       >  2:955694/1‑90 (MQ=255)
                                                            cgcgTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTcaac                     >  1:1769385/1‑90 (MQ=255)
                                                            cgcgTAAAGGCATCAAATTATGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTcaac                     >  2:943571/1‑90 (MQ=255)
                                                              cgTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTcaacaa                   >  2:2206572/1‑90 (MQ=255)
                                                              cgTAAAGGCATCAAATTCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTcaacaa                   >  1:1614250/1‑90 (MQ=255)
                                                                            aattcTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATc      <  2:713441/86‑1 (MQ=255)
                                                                            aattcTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATc      <  1:2405940/86‑1 (MQ=255)
                                                                             attcTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCt     >  2:537131/4‑90 (MQ=255)
                                                                             attcTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCt     <  1:2003993/87‑1 (MQ=255)
                                                                                cTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTgcg  >  1:1755751/1‑90 (MQ=255)
                                                                                cTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTgcg  >  1:571774/1‑90 (MQ=255)
                                                                                cTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTgcg  <  2:80613/90‑1 (MQ=255)
                                                                                cTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTgcg  <  1:3251629/90‑1 (MQ=255)
                                                                              |||                                                                                         
GCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATATGCTCGATAAGCTCGACGCCGCCTTACGTTTTCAACAAGAGGCGCTCAATCTGCG  >  AC_000091/1132483‑1132650

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