Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,411,637 Δ1 bp 100% coding (453/651 nt) yfbT ← sugar phosphatas

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,411,6340C.98.0% 461.7 / 1.7 102coding (456/651 nt)yfbTsugar phosphatas
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (44/56);  total (45/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 115 total aligned reads displayed.

TGCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGCCC  >  NC_000913/2411535‑2411729
                                                                                                   |                                                                                               
tGCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGc                                                                                                 <  1:1967583/100‑1 (MQ=255)
tGCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGc                                                                                                 >  2:226460/1‑100 (MQ=255)
 gCCGGAGCGTTAACCGCAATGACCTGACAACCCGCCGCCAGTCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCACGcc                                                                                                >  1:436766/1‑100 (MQ=255)
   cGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGccca                                                                                              >  2:1403747/1‑99 (MQ=255)
   cGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGccca                                                                                              <  1:490511/100‑2 (MQ=255)
   cGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGccca                                                                                              <  1:935642/100‑2 (MQ=255)
    ggAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGcccag                                                                                             <  1:465251/100‑3 (MQ=255)
    ggAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGcccag                                                                                             >  2:2117039/1‑98 (MQ=255)
    ggAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGccca                                                                                              >  2:744690/1‑98 (MQ=255)
        cGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCagcagc                                                                                        <  1:653322/100‑1 (MQ=255)
        cGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCagcagc                                                                                        <  2:1128749/100‑1 (MQ=255)
         gTTCACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCCGCGGGAGCATCTTCCCCCACCACACACTCCTGCGGCGGAAGcccagccg                                                                                          >  2:671832/1‑93 (MQ=255)
         gTTAACCGCAATGACATGACACCCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCt                                                                                       >  2:394312/1‑100 (MQ=255)
         gTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCt                                                                                       >  1:487296/1‑100 (MQ=255)
         gTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCt                                                                                       >  2:1358787/1‑100 (MQ=255)
            aaCCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTgcg                                                                                    <  1:678273/100‑1 (MQ=255)
            aaCCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTgcg                                                                                    >  2:1410806/1‑100 (MQ=255)
                gCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGcc                                                                                <  2:2048774/100‑1 (MQ=255)
                gCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGcc                                                                                <  2:676260/100‑1 (MQ=255)
                gCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGcc                                                                                <  2:1449388/100‑1 (MQ=255)
                 cAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCt                                                                               <  1:2203702/100‑1 (MQ=255)
                    tgacatgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTaa                                                                             <  1:166469/100‑1 (MQ=255)
                    tgacatgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAAc                                                                            <  1:193559/100‑1 (MQ=255)
                     gacatgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCCGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACa                                                                           <  1:1747265/100‑1 (MQ=255)
                     gacatgacaACCCGCCGCAAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACa                                                                           >  1:820463/1‑100 (MQ=255)
                      acatgacaACCCTCCTTCAGTCTAGAAAGCACGCCAGCGGGAGTATCTTCCACCACCACACATTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCTTAACa                                                                           <  2:1029802/99‑1 (MQ=255)
                      acatgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAg                                                                          <  1:673018/100‑1 (MQ=255)
                       catgaccACCCACCGCAATGCCAGAAAGCACGCCAGCCCGAGCAACTCCCACCACCACACACTCCTGCGGCACAAG‑CCCAGCAGCTGCGCGCCt                                                                               >  2:2327854/1‑94 (MQ=255)
                       catgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGa                                                                         <  1:1169080/100‑1 (MQ=255)
                       catgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGa                                                                         <  2:2240759/100‑1 (MQ=255)
                         tgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATa                                                                       >  1:477771/1‑100 (MQ=255)
                         tgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATa                                                                       <  2:2275504/100‑1 (MQ=255)
                               cccGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATc                                                                 >  2:977353/1‑100 (MQ=255)
                               cccGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATc                                                                 >  2:2315433/1‑100 (MQ=255)
                                ccgccgCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCa                                                                >  1:2403798/1‑100 (MQ=255)
                                 cgccgcCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAg                                                               <  1:1216652/100‑1 (MQ=255)
                                 cgccgcCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAg                                                               >  1:1582366/1‑100 (MQ=255)
                                  gccgccAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAgg                                                              <  1:45548/100‑1 (MQ=255)
                                   ccgccAGGCCAGAAAGCACGCCAGCGTGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGt                                                             <  2:803471/100‑1 (MQ=255)
                                    cgccAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGtt                                                            >  1:1597788/1‑100 (MQ=255)
                                     gccaggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTc                                                           >  1:2268278/1‑100 (MQ=255)
                                     gccaggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTc                                                           >  2:578709/1‑100 (MQ=255)
                                      ccaggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCt                                                          >  2:1204406/1‑100 (MQ=255)
                                       caggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTg                                                         <  2:545490/100‑1 (MQ=255)
                                        aggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTgg                                                        >  2:450747/1‑100 (MQ=255)
                                         ggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGt                                                       >  2:767576/1‑100 (MQ=255)
                                          gccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGt                                                       <  2:2128334/99‑1 (MQ=255)
                                           ccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGttt                                                     >  2:2333762/1‑100 (MQ=255)
                                               aaggCACGCCAGCGGGAGCATCTTCCACCAGCACGCACTCCGGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCcg                                                 <  2:182543/97‑1 (MQ=255)
                                               aaaGCACGACCGCGGCAGCATCTTCCACCACCACACCCTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACTGATACGCATCAGGTTCTGGTTTTCcg                                                 <  2:2032516/100‑1 (MQ=255)
                                                aaGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCcgc                                                <  1:2022741/100‑1 (MQ=255)
                                                 aGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCcgc                                                <  2:2298651/99‑1 (MQ=255)
                                                   cACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCt                                             <  2:487291/100‑1 (MQ=255)
                                                   cACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCt                                             <  2:16453/100‑1 (MQ=255)
                                                   cACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCt                                             <  2:2474779/100‑1 (MQ=255)
                                                     cGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTc                                           >  1:2469749/1‑100 (MQ=255)
                                                      gCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCa                                          <  1:308360/100‑1 (MQ=255)
                                                       ccAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCAc                                         <  1:1358796/100‑1 (MQ=255)
                                                        cAGCGGGAGCATCTTCCACCACCACACACTCCTGTGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACt                                        <  1:725659/100‑1 (MQ=255)
                                                          gCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCg                                      <  2:2465721/100‑1 (MQ=255)
                                                            gggAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCt                                    <  1:849867/100‑1 (MQ=255)
                                                             ggAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTc                                   >  2:2262378/1‑100 (MQ=255)
                                                                  aTCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCgg                              <  2:637591/100‑1 (MQ=255)
                                                                   tCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATGCGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGt                             >  2:1146734/1‑100 (MQ=255)
                                                                   tCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGt                             >  1:190870/1‑100 (MQ=255)
                                                                   tCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGt                             >  1:388507/1‑100 (MQ=255)
                                                                   tCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTGCACTCGCTCAGCgg                              >  1:1507416/1‑99 (MQ=255)
                                                                      tCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTa                           >  2:1369120/1‑100 (MQ=255)
                                                                      tCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTAc                          >  1:2043998/1‑100 (MQ=255)
                                                                       ccaccaccaCACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACa                         <  2:1159794/100‑1 (MQ=255)
                                                                        caccaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACaa                        <  1:1099585/100‑1 (MQ=255)
                                                                          ctcttccgatcTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAac                      <  1:1243477/90‑1 (MQ=255)
                                                                          ccaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAac                      <  1:2418406/100‑1 (MQ=255)
                                                                          ccaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAac                      <  1:2383168/100‑1 (MQ=255)
                                                                          ccaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAac                      <  1:744667/100‑1 (MQ=255)
                                                                           caccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAaca                     <  1:1572429/100‑1 (MQ=255)
                                                                             ccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACAcc                   <  1:1583402/100‑1 (MQ=255)
                                                                             ccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACAcc                   <  2:1744105/100‑1 (MQ=255)
                                                                             ccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACAcc                   <  2:1422557/100‑1 (MQ=255)
                                                                              cacacaCTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACct                  >  1:2064264/1‑100 (MQ=255)
                                                                               acacacTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATCCGCACCAGGTTCTGCTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACctc                 >  2:1235691/1‑100 (MQ=255)
                                                                                cacacTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACctct                <  1:1891052/100‑1 (MQ=255)
                                                                                  nacTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTgg              >  2:1822379/2‑100 (MQ=255)
                                                                                  cacTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAGCAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTgg              >  1:1805137/1‑100 (MQ=255)
                                                                                  cacTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTgg              >  2:35910/1‑100 (MQ=255)
                                                                                      ccTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCg          <  1:1616757/100‑1 (MQ=255)
                                                                                      ccTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCg          >  1:1360422/1‑100 (MQ=255)
                                                                                      ccTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCg          >  1:89845/1‑100 (MQ=255)
                                                                                      ccTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCg          >  1:90642/1‑100 (MQ=255)
                                                                                       cTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAACCACCTCTGGTGCgg         >  2:2073711/1‑100 (MQ=255)
                                                                                       cTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCgg         <  1:662410/100‑1 (MQ=255)
                                                                                       cTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCgg         <  1:1590393/100‑1 (MQ=255)
                                                                                          cggcgCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGaa      <  1:1767574/100‑1 (MQ=255)
                                                                                          cggcgCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGaa      <  1:2069919/100‑1 (MQ=255)
                                                                                          cggcgCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGaa      >  2:1000292/1‑100 (MQ=255)
                                                                                           ggcgCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAg     <  1:35911/100‑1 (MQ=255)
                                                                                           ggcgCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACCCGCTCAGCGGGTACAAACACCTCTGGTGCGGGAAg     >  2:1342507/1‑100 (MQ=255)
                                                                                             cgcAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGcc   >  2:136859/1‑100 (MQ=255)
                                                                                              gcAAG‑CCCAGCAGCTGCGCGCCTAACAGTTACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGccc  >  2:2246672/1‑100 (MQ=255)
                                                                                              gcAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGccc  <  1:583337/100‑1 (MQ=255)
                                                                                                   |                                                                                               
TGCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGCCC  >  NC_000913/2411535‑2411729

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: ‑