Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,088,657 T→C I431V (ATC→GTC)  frvB ← PTS system, fructose‑like enzyme IIBC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,088,6570TC99.5% 685.1 / ‑5.3 218I431V (ATC→GTC) frvBPTS system, fructose‑like enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (120/97);  total (120/98)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 234 total aligned reads displayed.

TAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCG  >  NC_000913/4088567‑4088746
                                                                                          |                                                                                         
tAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGacac                                                                                         >  1:4273820/1‑93 (MQ=255)
tAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGacac                                                                                         >  1:1381187/1‑93 (MQ=255)
 aatgaatgAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACAcc                                                                                        >  2:4094387/1‑93 (MQ=255)
  atgaatgaAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACAccc                                                                                       <  2:722940/93‑1 (MQ=255)
   tgaatgaaATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCa                                                                                      >  2:1447215/1‑93 (MQ=255)
    gaatgaaATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCaa                                                                                     >  2:3186232/1‑93 (MQ=255)
    gaatgaaATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCaa                                                                                     <  2:420086/93‑1 (MQ=255)
     aatgaaATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCaaa                                                                                    <  2:2583932/93‑1 (MQ=255)
     aatgaaATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCaaa                                                                                    >  1:174854/1‑93 (MQ=255)
     aatgaaATGCTCAGCGCGAGGCCAATCGCCAGCCCCAGATAAAACGAACCCCCCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCaaa                                                                                    <  2:866308/93‑1 (MQ=255)
      atgaaATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAAc                                                                                   >  2:2829839/1‑93 (MQ=255)
      atgaaATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAAc                                                                                   >  1:1905643/1‑93 (MQ=255)
       tgaaATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAAcc                                                                                  <  2:1590540/93‑1 (MQ=255)
       tgaaATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAAcc                                                                                  >  1:458843/1‑93 (MQ=255)
       tgaaATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAAcc                                                                                  >  1:76171/1‑93 (MQ=255)
        gaaATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCa                                                                                 >  1:642284/1‑93 (MQ=255)
        gaaATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCa                                                                                 >  1:4190984/1‑93 (MQ=255)
         aaaTTTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAg                                                                                <  1:562042/93‑1 (MQ=255)
         aaaTGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAg                                                                                >  1:2382090/1‑93 (MQ=255)
         aaaTGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAg                                                                                >  2:2000303/1‑93 (MQ=255)
           aTGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAggg                                                                              >  2:3260823/1‑92 (MQ=255)
           aTGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGt                                                                              >  1:1004885/1‑93 (MQ=255)
           aTGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGt                                                                              >  1:2962696/1‑93 (MQ=255)
             gTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGCCAAAGACACCCAAACCAGGTGc                                                                            >  2:3032378/1‑93 (MQ=255)
              ttCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGcc                                                                           >  2:3041963/1‑93 (MQ=255)
              ttCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGcc                                                                           >  2:822915/1‑93 (MQ=255)
                cAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAg                                                                         >  1:2076436/1‑93 (MQ=255)
                 aGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGa                                                                        <  2:2132063/93‑1 (MQ=255)
                  gcgcGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGAc                                                                       <  1:2842283/93‑1 (MQ=255)
                    gcgAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGt                                                                     <  2:345588/93‑1 (MQ=255)
                     cgAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGtt                                                                    >  2:4644363/1‑93 (MQ=255)
                     cgAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGtt                                                                    <  2:3300976/93‑1 (MQ=255)
                     cgAGACCAATCGCAAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGtt                                                                    <  1:1088053/93‑1 (MQ=255)
                      gAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGttt                                                                   <  2:360796/93‑1 (MQ=255)
                           cAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTc                                                              <  1:4482130/93‑1 (MQ=255)
                            aaTCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTcc                                                             <  2:2260607/93‑1 (MQ=255)
                            aaTCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTcc                                                             >  1:447780/1‑93 (MQ=255)
                              tCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGa                                                           <  2:447847/93‑1 (MQ=255)
                                gCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGaaa                                                         <  1:889179/93‑1 (MQ=255)
                                cccAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGCTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGaaa                                                         <  1:3302153/92‑1 (MQ=255)
                                 ccAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAg                                                        >  2:2153287/1‑93 (MQ=255)
                                 ccAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAg                                                        <  2:4274119/93‑1 (MQ=255)
                                  cAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGc                                                       <  2:3491083/93‑1 (MQ=255)
                                   aGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACGAGGTGCCAGACGTTTGATTCCGAAACCg                                                      >  1:1514383/1‑93 (MQ=255)
                                    gCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCCGACGTTTGATTCCGAAAGCGa                                                     >  1:2825394/1‑93 (MQ=255)
                                    gCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGa                                                     >  1:3290545/1‑93 (MQ=255)
                                    gCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGa                                                     <  1:2433144/93‑1 (MQ=255)
                                     cACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGAt                                                    <  1:2785751/93‑1 (MQ=255)
                                       ccAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCa                                                  <  2:3348173/93‑1 (MQ=255)
                                        cAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACGCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCaa                                                 >  1:1715519/1‑92 (MQ=255)
                                          gATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCa                                               >  2:2547338/1‑93 (MQ=255)
                                           aTAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAg                                              >  2:2698709/1‑93 (MQ=255)
                                           aTAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAg                                              >  2:1969483/1‑93 (MQ=255)
                                            tAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGc                                             >  2:603611/1‑93 (MQ=255)
                                                aCGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGc                                         <  1:304529/93‑1 (MQ=255)
                                                   aaCCCACCGGCGACATCAGGCCGAGCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCgg                                      >  1:3644503/1‑93 (MQ=255)
                                                    aCCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGATTGATTCCGAAAGCGATCACCAGCACGCCGGt                                     >  1:1877119/1‑93 (MQ=255)
                                                      ccACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTaa                                   <  2:2491800/93‑1 (MQ=255)
                                                       cACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGTTCACCAGCACGCCGGTAAt                                  <  2:4564407/93‑1 (MQ=255)
                                                       cACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAAt                                  <  2:2741582/93‑1 (MQ=255)
                                                       cACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAAt                                  <  1:4022811/93‑1 (MQ=255)
                                                       cACCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAAt                                  <  1:3218758/93‑1 (MQ=255)
                                                        aCCGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATg                                 >  2:3628011/1‑93 (MQ=255)
                                                         ccGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATgc                                >  2:1137323/1‑93 (MQ=255)
                                                         ccGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATgc                                <  1:1175888/93‑1 (MQ=255)
                                                          cGGCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATAACCAGCACGCCGGTAATgcc                               >  2:193493/1‑93 (MQ=255)
                                                            gCGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATgccgc                             >  1:984401/1‑93 (MQ=255)
                                                             cGACATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATgccgcc                            >  2:1232745/1‑93 (MQ=255)
                                                               aCATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGccgcccg                          >  1:79324/1‑93 (MQ=255)
                                                                cATCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGccgcccgc                         >  1:3711077/1‑93 (MQ=255)
                                                                 aTCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCcgcccgcc                        >  1:1039517/1‑93 (MQ=255)
                                                                 aTCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCcgcccgcc                        >  1:2618013/1‑93 (MQ=255)
                                                                  tCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCa                       <  1:632262/93‑1 (MQ=255)
                                                                  tCAGGCCGATCAGCGGGTCAAAGACACCCAAACCAGGCGCCAGACGTTTGATTCCGAAAGCGATCACCACCACGACGGTAATGCCGCCCGCCa                       >  2:520287/1‑93 (MQ=255)
                                                                     ggCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCg                    <  1:4344730/93‑1 (MQ=255)
                                                                     ggCCGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCg                    >  2:4358644/1‑93 (MQ=255)
                                                                       ccGATCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTa                  >  2:3516600/1‑93 (MQ=255)
                                                                         gaccAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTAtt                <  2:1816013/90‑1 (MQ=255)
                                                                         gATCAGGGGGTCAAAGCCCCCAAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGACCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTAtt                <  2:3951727/93‑1 (MQ=255)
                                                                          aTCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTg               >  2:1679303/1‑93 (MQ=255)
                                                                           tCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTgg              >  1:183158/1‑93 (MQ=255)
                                                                           tCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTgg              >  2:3294347/1‑93 (MQ=255)
                                                                           tCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTgg              >  1:4639071/1‑93 (MQ=255)
                                                                               cgagTCAAAGACACCCAACCCGGGTCCAGGACGTTAAATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGt          <  2:1645940/90‑1 (MQ=255)
                                                                               cGGGTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGt          <  1:452046/93‑1 (MQ=255)
                                                                                  gTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTaca       >  2:2771970/1‑91 (MQ=255)
                                                                                  gTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAat       <  2:2040775/93‑1 (MQ=255)
                                                                                  gTCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAat       >  1:2123315/1‑93 (MQ=255)
                                                                                   tCAAAGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAatc      <  2:2366284/93‑1 (MQ=255)
                                                                                     aaaGACACCCAAAGCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCCCGCCGGTAATGCCCCCCGCCAGCGTATTGGCGGTAatcat    <  2:3827050/93‑1 (MQ=255)
                                                                                     aaaGACACCCAAACCAGGTGCCAGACGTTTGATTACGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAatcat    <  1:1650998/93‑1 (MQ=255)
                                                                                      aaGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAatcatc   >  1:2614252/1‑93 (MQ=255)
                                                                                      aaGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAatcatc   <  1:2230055/93‑1 (MQ=255)
                                                                                      aaGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAatcatc   <  1:3518338/93‑1 (MQ=255)
                                                                                      aaGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAatcatc   <  1:4438359/93‑1 (MQ=255)
                                                                                      aaGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAatcatc   >  1:742812/1‑93 (MQ=255)
                                                                                      aaGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAatcatc   <  2:1444873/93‑1 (MQ=255)
                                                                                      aaGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAatcatc   <  2:1154001/93‑1 (MQ=255)
                                                                                       aGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCg  >  1:2319898/1‑93 (MQ=255)
                                                                                       aGACACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCg  <  1:1849877/93‑1 (MQ=255)
                                                                                          |                                                                                         
TAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCG  >  NC_000913/4088567‑4088746

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