Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,807,923 (T)6→5 100% coding (128/861 nt) yniA → fructosamine kinase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,807,9180T.83.3% 238.4 / 56.4 84coding (123/861 nt)yniAfructosamine kinase family protein
Reads supporting (aligned to +/- strand):  ref base T (10/4);  new base . (41/29);  total (51/33)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCACCGT  >  NC_000913/1807820‑1808013
                                                                                                  |                                                                                               
ttGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGAtt                                                                                                <  2:184268/100‑1 (MQ=255)
 tGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGAttt                                                                                               >  1:517020/1‑100 (MQ=255)
 tGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGAtt                                                                                                >  2:806536/1‑99 (MQ=255)
  gagcgagcAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGAtttt                                                                                              <  2:1538315/100‑1 (MQ=255)
   agcgagcAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGAttttt                                                                                             >  1:844985/1‑100 (MQ=255)
    gcgagcAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGAtttttc                                                                                            >  2:2368003/1‑99 (MQ=255)
     cgagcAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGAtttttcg                                                                                           <  1:1459081/100‑3 (MQ=255)
     cgagcAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTc                                                                                           <  1:997286/100‑1 (MQ=255)
      gagcAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAATTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGAtttttcgt                                                                                          <  1:831400/100‑4 (MQ=255)
      gagcAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGAtttttcgt                                                                                          <  2:1690833/100‑4 (MQ=255)
      gagcAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGAtttttcgt                                                                                          >  1:1855149/1‑97 (MQ=255)
       agcAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGAtttttcgtc                                                                                         <  1:1422496/100‑5 (MQ=255)
         cAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCAATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCaa                                                                                      >  1:1883489/1‑100 (MQ=255)
          aGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCaaa                                                                                     >  1:1393448/1‑100 (MQ=255)
          aGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCaaa                                                                                     >  1:1108538/1‑100 (MQ=255)
           gTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAAt                                                                                    >  2:1285492/1‑100 (MQ=255)
             tAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAAtgt                                                                                  >  2:1027293/1‑100 (MQ=255)
              aGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAAtgt                                                                                  >  1:383063/1‑100 (MQ=255)
              aGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAAtgtg                                                                                 <  2:1595958/100‑1 (MQ=255)
               ggTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGtga                                                                                >  2:1014448/1‑100 (MQ=255)
                gTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTGTTCGTc                                                                                        >  1:2197341/1‑92 (MQ=255)
                 tGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGtgatg                                                                              >  2:1733656/1‑100 (MQ=255)
                  gAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGtgatga                                                                             >  2:2568559/1‑100 (MQ=255)
                     ggCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAGAtcggaagagca                                                                          >  2:2288499/1‑90 (MQ=255)
                     ggCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAg                                                                          >  1:1009166/1‑100 (MQ=255)
                     ggCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAg                                                                          >  1:2090196/1‑100 (MQ=255)
                      gCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGaaa                                                                           >  1:737573/1‑98 (MQ=255)
                      gCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAATgg                                                                         >  1:274091/1‑100 (MQ=255)
                        gAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAggg                                                                        >  2:1262300/1‑100 (MQ=255)
                        gAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGa                                                                       <  2:1073076/100‑1 (MQ=255)
                          aaTCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAAc                                                                     >  2:2030315/1‑100 (MQ=255)
                             cGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTg                                                                   >  1:1700173/1‑100 (MQ=255)
                             cGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGc                                                                  >  2:61753/1‑100 (MQ=255)
                                   gCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCg                                                            >  2:987385/1‑100 (MQ=255)
                                       aaTGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGtt                                                         <  2:383065/100‑1 (MQ=255)
                                       aaTGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGttt                                                        <  2:25313/100‑1 (MQ=255)
                                         tGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGtttt                                                       >  2:1768574/1‑100 (MQ=255)
                                           aaCTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAcc                                                    >  1:2011161/1‑100 (MQ=255)
                                           aaCTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAcc                                                    <  2:358235/100‑1 (MQ=255)
                                            aCTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccg                                                   >  1:1861972/1‑100 (MQ=255)
                                             cTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccgc                                                  >  1:76513/1‑100 (MQ=255)
                                              tGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGCAACTGCTTCCCGGTTTTAccgc                                                  <  1:1503927/100‑1 (MQ=255)
                                              tGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccgcc                                                 >  2:1789659/1‑100 (MQ=255)
                                              tGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccgcc                                                 >  1:41885/1‑100 (MQ=255)
                                               gCCTGGCGGAGAAGTCCATTCTGCATGGTTTTTTCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccgccg                                                >  1:1564985/1‑100 (MQ=255)
                                               gCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGGCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccgccg                                                >  2:1544181/1‑100 (MQ=255)
                                                 cTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGa                                               >  2:1556250/1‑100 (MQ=255)
                                                 cTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGaa                                              >  2:249545/1‑100 (MQ=255)
                                                   ggCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAg                                             >  2:173631/1‑100 (MQ=255)
                                                   ggCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGc                                            <  1:1265479/100‑1 (MQ=255)
                                                   ggCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGc                                            >  2:18181/1‑100 (MQ=255)
                                                   ggCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGc                                            <  2:1257169/100‑1 (MQ=255)
                                                   ggCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGc                                            <  2:561037/100‑1 (MQ=255)
                                                    gCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGcc                                           >  1:1520298/1‑100 (MQ=255)
                                                    gCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGcc                                           <  2:655101/100‑1 (MQ=255)
                                                     cGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCg                                          >  2:1607798/1‑100 (MQ=255)
                                                     cGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCg                                          >  2:1854185/1‑100 (MQ=255)
                                                       gagaAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGa                                         <  1:2259988/100‑1 (MQ=255)
                                                       cttccgatcttGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGAc                                        <  1:485886/90‑1 (MQ=255)
                                                        agaAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGAcc                                       <  2:1730227/100‑1 (MQ=255)
                                                           aGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAAc                                    <  1:717590/100‑1 (MQ=255)
                                                            gTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACt                                   <  2:1855159/100‑1 (MQ=255)
                                                             tCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTg                                  <  2:416457/100‑1 (MQ=255)
                                                             tCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTg                                  <  1:1607789/100‑1 (MQ=255)
                                                             tCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTg                                  <  1:1482339/100‑1 (MQ=255)
                                                              ccATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTgg                                 >  2:234421/1‑100 (MQ=255)
                                                                aTGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTTTTTTTTTTTATGAAAGGGAACTGCTTCCCAGTTTTACCGCCGAAGCCGACCAACTGGAg                               >  1:1783398/1‑100 (MQ=255)
                                                                aTGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAg                               <  2:779166/100‑1 (MQ=255)
                                                                 tGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGAcagatcggaag                              >  2:485888/1‑90 (MQ=255)
                                                                    cGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTAc                           <  2:919686/100‑1 (MQ=255)
                                                                    cGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTAc                           <  2:921642/100‑1 (MQ=255)
                                                                    cGCATGGCATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTAc                           <  1:1778036/100‑1 (MQ=255)
                                                                      cATGGCATTTGCGCTATGCAGGACATGATTTTCCGTCAATTGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTg                         <  1:1818194/100‑1 (MQ=255)
                                                                           cATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTcgcg                    <  2:1948288/100‑1 (MQ=255)
                                                                           cATTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTcgcg                    >  1:2183264/1‑100 (MQ=255)
                                                                            aTTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgt                   <  2:2276515/100‑1 (MQ=255)
                                                                            aTTTGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgt                   <  1:1758404/100‑1 (MQ=255)
                                                                              ttGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgtag                 <  2:679138/100‑1 (MQ=255)
                                                                               tGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgtag                 >  1:2264447/1‑100 (MQ=255)
                                                                               tGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgtagt                <  2:2394375/100‑1 (MQ=255)
                                                                               tGCGCTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgtagt                >  2:2416008/1‑100 (MQ=255)
                                                                                gcgcTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgtagta               >  1:744930/1‑100 (MQ=255)
                                                                                  gcTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTaaa             >  2:1794805/1‑100 (MQ=255)
                                                                                   cTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTaaaa            >  2:884429/1‑100 (MQ=255)
                                                                                   cTATGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTaaaa            >  2:883830/1‑100 (MQ=255)
                                                                                     aTGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAAc           >  2:2339465/1‑100 (MQ=255)
                                                                                      tGCAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCg         <  1:589625/100‑1 (MQ=255)
                                                                                        cAGGACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgtagatcggaagag       >  1:1080533/1‑90 (MQ=255)
                                                                                          tgACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCTGTTTTTACCGCGGCAGAAGACAAACAGGCGTTACTGTCGCGTAGTAAAACCg         >  1:980414/2‑96 (MQ=255)
                                                                                          ggACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCAc     <  2:2307599/100‑1 (MQ=255)
                                                                                          ggACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCAc     <  1:1829960/100‑1 (MQ=255)
                                                                                           gACATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCAcc    <  1:2446535/100‑1 (MQ=255)
                                                                                             cATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCACCGt  >  2:427414/1‑100 (MQ=255)
                                                                                             cATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCACCGt  >  2:233484/1‑100 (MQ=255)
                                                                                             cATGATTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCACCGt  >  2:1039680/1‑100 (MQ=255)
                                                                                                  |                                                                                               
TTGAGCGAGCAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCACCGT  >  NC_000913/1807820‑1808013

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