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)
         cAGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCAATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCaa                                                                                      >  1:1883489/1‑100 (MQ=255)
          aGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCaaa                                                                                     >  1:1393448/1‑100 (MQ=255)
          aGTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCaaa                                                                                     >  1:1108538/1‑100 (MQ=255)
           gTTAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAAt                                                                                    >  2:1285492/1‑100 (MQ=255)
             tAGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAAtgt                                                                                  >  2:1027293/1‑100 (MQ=255)
              aGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAAtgt                                                                                  >  1:383063/1‑100 (MQ=255)
              aGGTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAAtgtg                                                                                 <  2:1595958/100‑1 (MQ=255)
               ggTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGtga                                                                                >  2:1014448/1‑100 (MQ=255)
                gTGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTGTTCGTc                                                                                        >  1:2197341/1‑92 (MQ=255)
                 tGAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGtgatg                                                                              >  2:1733656/1‑100 (MQ=255)
                  gAAGGCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGtgatga                                                                             >  2:2568559/1‑100 (MQ=255)
                     ggCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAGAtcggaagagca                                                                          >  2:2288499/1‑90 (MQ=255)
                     ggCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAg                                                                          >  1:1009166/1‑100 (MQ=255)
                     ggCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAg                                                                          >  1:2090196/1‑100 (MQ=255)
                      gCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGaaa                                                                           >  1:737573/1‑98 (MQ=255)
                      gCGAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAATgg                                                                         >  1:274091/1‑100 (MQ=255)
                        gAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAggg                                                                        >  2:1262300/1‑100 (MQ=255)
                        gAAATCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGa                                                                       <  2:1073076/100‑1 (MQ=255)
                          aaTCGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAAc                                                                     >  2:2030315/1‑100 (MQ=255)
                             cGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTg                                                                   >  1:1700173/1‑100 (MQ=255)
                             cGAACTGCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGc                                                                  >  2:61753/1‑100 (MQ=255)
                                   gCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCg                                                            >  2:987385/1‑100 (MQ=255)
                                       aaTGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGtt                                                         <  2:383065/100‑1 (MQ=255)
                                       aaTGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGttt                                                        <  2:25313/100‑1 (MQ=255)
                                         tGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGtttt                                                       >  2:1768574/1‑100 (MQ=255)
                                           aaCTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAcc                                                    >  1:2011161/1‑100 (MQ=255)
                                           aaCTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAcc                                                    <  2:358235/100‑1 (MQ=255)
                                            aCTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccg                                                   >  1:1861972/1‑100 (MQ=255)
                                             cTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccgc                                                  >  1:76513/1‑100 (MQ=255)
                                              tGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGCAACTGCTTCCCGGTTTTAccgc                                                  <  1:1503927/100‑1 (MQ=255)
                                              tGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccgcc                                                 >  2:1789659/1‑100 (MQ=255)
                                              tGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccgcc                                                 >  1:41885/1‑100 (MQ=255)
                                               gCCTGGCGGAGAAGTCCATTCTGCATGGTTTTTTCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccgccg                                                >  1:1564985/1‑100 (MQ=255)
                                               gCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGGCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccgccg                                                >  2:1544181/1‑100 (MQ=255)
                                                 cTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGa                                               >  2:1556250/1‑100 (MQ=255)
                                                 cTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGaa                                              >  2:249545/1‑100 (MQ=255)
                                                   ggCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAg                                             >  2:173631/1‑100 (MQ=255)
                                                   ggCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGc                                            <  1:1265479/100‑1 (MQ=255)
                                                   ggCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGc                                            >  2:18181/1‑100 (MQ=255)
                                                   ggCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGc                                            <  2:1257169/100‑1 (MQ=255)
                                                   ggCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGc                                            <  2:561037/100‑1 (MQ=255)
                                                    gCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGcc                                           >  1:1520298/1‑100 (MQ=255)
                                                    gCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGcc                                           <  2:655101/100‑1 (MQ=255)
                                                     cGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCg                                          >  2:1607798/1‑100 (MQ=255)
                                                     cGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCg                                          >  2:1854185/1‑100 (MQ=255)
                                                       gagaAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGa                                         <  1:2259988/100‑1 (MQ=255)
                                                       cttccgatcttGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGAc                                        <  1:485886/90‑1 (MQ=255)
                                                        agaAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGAcc                                       <  2:1730227/100‑1 (MQ=255)
                                                           aGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAAc                                    <  1:717590/100‑1 (MQ=255)
                                                            gTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACt                                   <  2:1855159/100‑1 (MQ=255)
                                                             tCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTg                                  <  2:416457/100‑1 (MQ=255)
                                                             tCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTg                                  <  1:1607789/100‑1 (MQ=255)
                                                             tCCATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTg                                  <  1:1482339/100‑1 (MQ=255)
                                                              ccATGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTgg                                 >  2:234421/1‑100 (MQ=255)
                                                                aTGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTT‑TTTTTTTTTTATGAAAGGGAACTGCTTCCCAGTTTTACCGCCGAAGCCGACCAACTGGAg                               >  1:1783398/1‑100 (MQ=255)
                                                                aTGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAg                               <  2:779166/100‑1 (MQ=255)
                                                                 tGCCGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGAcagatcggaag                              >  2:485888/1‑90 (MQ=255)
                                                                    cGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTAc                           <  2:919686/100‑1 (MQ=255)
                                                                    cGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTAc                           <  2:921642/100‑1 (MQ=255)
                                                                    cGCATGGCATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTAc                           <  1:1778036/100‑1 (MQ=255)
                                                                      cATGGCATTTGCGCTATGCAGGACATGA‑TTTTCCGTCAATTGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTg                         <  1:1818194/100‑1 (MQ=255)
                                                                           cATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTcgcg                    <  2:1948288/100‑1 (MQ=255)
                                                                           cATTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTcgcg                    >  1:2183264/1‑100 (MQ=255)
                                                                            aTTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgt                   <  2:2276515/100‑1 (MQ=255)
                                                                            aTTTGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgt                   <  1:1758404/100‑1 (MQ=255)
                                                                              ttGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgtag                 <  2:679138/100‑1 (MQ=255)
                                                                               tGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgtag                 >  1:2264447/1‑100 (MQ=255)
                                                                               tGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgtagt                <  2:2394375/100‑1 (MQ=255)
                                                                               tGCGCTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgtagt                >  2:2416008/1‑100 (MQ=255)
                                                                                gcgcTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgtagta               >  1:744930/1‑100 (MQ=255)
                                                                                  gcTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTaaa             >  2:1794805/1‑100 (MQ=255)
                                                                                   cTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTaaaa            >  2:884429/1‑100 (MQ=255)
                                                                                   cTATGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTaaaa            >  2:883830/1‑100 (MQ=255)
                                                                                     aTGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAAc           >  2:2339465/1‑100 (MQ=255)
                                                                                      tGCAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCg         <  1:589625/100‑1 (MQ=255)
                                                                                        cAGGACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCgtagatcggaagag       >  1:1080533/1‑90 (MQ=255)
                                                                                          tgACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCTGTTTTTACCGCGGCAGAAGACAAACAGGCGTTACTGTCGCGTAGTAAAACCg         >  1:980414/2‑96 (MQ=255)
                                                                                          ggACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCAc     <  2:2307599/100‑1 (MQ=255)
                                                                                          ggACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCAc     <  1:1829960/100‑1 (MQ=255)
                                                                                           gACATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCAcc    <  1:2446535/100‑1 (MQ=255)
                                                                                             cATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCACCGt  >  2:427414/1‑100 (MQ=255)
                                                                                             cATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCACCGt  >  2:233484/1‑100 (MQ=255)
                                                                                             cATGA‑TTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAACTGGAGTTACTGTCGCGTAGTAAAACCGTCACCGt  >  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: ‑