Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,061,235 A→C Y24S (TAT→TCT)  yihM → putative sugar phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,061,2350AC99.4% 547.1 / NA 160Y24S (TAT→TCT) yihMputative sugar phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (77/81);  total (79/81)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

ACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGT  >  NC_000913/4061138‑4061331
                                                                                                 |                                                                                                
 cgatctcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGc                                                                                               <  2:483131/95‑1 (MQ=255)
  cTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGcc                                                                                              >  1:2286732/1‑100 (MQ=255)
  cTCAACATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGcc                                                                                              >  2:3379032/1‑100 (MQ=255)
   tcatcaTTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCt                                                                                             >  1:15167/1‑100 (MQ=255)
    catcatTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCtt                                                                                            >  2:1317473/1‑100 (MQ=255)
       catTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCa                                                                                         >  1:296732/1‑100 (MQ=255)
       catTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCa                                                                                         >  1:3489295/1‑100 (MQ=255)
       catTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCa                                                                                         <  1:1356685/100‑1 (MQ=255)
       catTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCa                                                                                         <  2:1508510/100‑1 (MQ=255)
        atTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCAt                                                                                        <  1:1427223/100‑1 (MQ=255)
        atTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCAt                                                                                        <  1:2061113/100‑1 (MQ=255)
        atTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCAt                                                                                        <  1:1065685/100‑1 (MQ=255)
        atTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCAt                                                                                        >  2:345289/1‑100 (MQ=255)
           aaTCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTa                                                                                     >  1:695680/1‑100 (MQ=255)
            aTCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTaa                                                                                    >  1:984256/1‑100 (MQ=255)
            aTCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTaa                                                                                    >  2:2312976/1‑100 (MQ=255)
              cTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAat                                                                                  <  2:1781140/100‑1 (MQ=255)
              cTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAat                                                                                  <  2:3998232/100‑1 (MQ=255)
               tGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAata                                                                                 >  1:1941195/1‑100 (MQ=255)
               tGTAGGTGTTTTATCGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAata                                                                                 <  2:3653850/100‑1 (MQ=255)
               tGTAGGTGGTTTATGGTGACAATAAATAACGCAAGAAAGATTCAACAACGTGTCGACACTCTTCCTCTTTATTTACACGCTTCTGCCCTTCATTTAAata                                                                                 >  2:3122342/1‑100 (MQ=255)
                 tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATg                                                                               >  1:2619588/1‑100 (MQ=255)
                 tAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATg                                                                               >  2:276012/1‑100 (MQ=255)
                   ggTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCg                                                                             >  1:2750780/1‑100 (MQ=255)
                   ggTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCg                                                                             >  1:1187923/1‑100 (MQ=255)
                      nttttATGGTGACAATAAATAACGCAAGAAAGATTATACAACGTGTCGACACTCTTCCTCTTTATTTACCTGCTTATGCCTTTCATATAACTATGCgccc                                                                          >  1:3245139/2‑97 (MQ=255)
                      gtTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCt                                                                          >  1:425779/1‑100 (MQ=255)
                        tttATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTgg                                                                        <  1:1069396/100‑1 (MQ=255)
                        tttATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTgg                                                                        >  2:1840012/1‑100 (MQ=255)
                             ggTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAaga                                                                   >  1:2320866/1‑100 (MQ=255)
                               tGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAgt                                                                 <  1:684292/100‑1 (MQ=255)
                               tGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAgt                                                                 >  2:876727/1‑100 (MQ=255)
                                gACAATAAATAACGCAATAAAAATTCTACAACGTGCTGACACTCTTCCTCTTTCTGAACATGCTTCTCCCTTTCACTTTTATATGCGGCTGGAAagagag                                                                <  2:616878/100‑3 (MQ=255)
                                gACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAgtg                                                                >  2:384051/1‑100 (MQ=255)
                                gACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAgtg                                                                <  2:399525/100‑1 (MQ=255)
                                gACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAgtg                                                                <  2:2339013/100‑1 (MQ=255)
                                gACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCCTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAgtg                                                                >  1:336315/1‑100 (MQ=255)
                                 aCAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAgtgt                                                               <  1:3829283/100‑1 (MQ=255)
                                 aCAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAgtgt                                                               <  2:2709055/100‑1 (MQ=255)
                                  cAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGtt                                                              >  2:366206/1‑100 (MQ=255)
                                  cAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGtt                                                              >  1:361852/1‑100 (MQ=255)
                                    ataaataaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGc                                                            >  2:1148588/1‑100 (MQ=255)
                                     taaataaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGcc                                                           <  1:1507943/100‑1 (MQ=255)
                                     taaataaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGcc                                                           <  1:839887/100‑1 (MQ=255)
                                      aaataaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCct                                                          >  2:2433643/1‑100 (MQ=255)
                                       aataaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctg                                                         >  2:3283333/1‑100 (MQ=255)
                                        ataaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctgc                                                        <  2:2355781/100‑1 (MQ=255)
                                          aaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctgctg                                                      >  1:1634804/1‑100 (MQ=255)
                                          aaCGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctgctg                                                      <  1:1453132/100‑1 (MQ=255)
                                             gCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGAtt                                                   <  2:3682346/100‑1 (MQ=255)
                                              cacgaaagaTTGTACAACCTGTCGACGCTCTTCCTCTTTATTTACATGATTCTACCCTTCATTTAAATATGCGGCTGGCGAGAGTGGTGCCGGCGGCttt                                                  >  2:966091/4‑100 (MQ=255)
                                                 gaaagaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACt                                               >  2:2353290/1‑100 (MQ=255)
                                                 gaaagaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACt                                               >  2:559821/1‑100 (MQ=255)
                                                 gaaagaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACt                                               <  2:2686074/100‑1 (MQ=255)
                                                 gaaagaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACt                                               >  1:3087330/1‑100 (MQ=255)
                                                 gaaagaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACt                                               >  1:2272373/1‑100 (MQ=255)
                                                  aaagaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTCATTTACtt                                              >  2:120221/1‑100 (MQ=255)
                                                   aagaTTCTACAACGTGTCGACCCTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTg                                             <  1:3201021/100‑1 (MQ=255)
                                                   aagaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTg                                             >  2:65672/1‑100 (MQ=255)
                                                   aagaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTg                                             <  1:4053866/100‑1 (MQ=255)
                                                     gaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGat                                           <  1:1915721/100‑1 (MQ=255)
                                                     gaTTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGat                                           <  1:2432497/100‑1 (MQ=255)
                                                        tCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATc                                        >  1:3655747/1‑100 (MQ=255)
                                                        tCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATc                                        >  1:3096536/1‑100 (MQ=255)
                                                        tCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCC‑TTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCg                                       >  1:3592600/1‑100 (MQ=255)
                                                         cTACAACGTGTCGACACTCTTCCTCCTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCg                                       <  1:2235703/100‑1 (MQ=255)
                                                           aCAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCa                                     <  1:378046/100‑1 (MQ=255)
                                                            cAACGTTTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCaa                                    <  2:3586628/100‑1 (MQ=255)
                                                             aaCGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAg                                   >  2:2935267/1‑100 (MQ=255)
                                                              aCGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGt                                  >  2:3531708/1‑100 (MQ=255)
                                                              aCGTGTCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGt                                  >  2:1296629/1‑100 (MQ=255)
                                                               cGTGTCTACACGCGGCCTCTGTAGTTACATTCTGCGGCCTTGCATTTAAATATTCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCTCAAGTg                                 <  1:1840620/100‑1 (MQ=255)
                                                                 tgtCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGaa                               >  2:3838394/1‑100 (MQ=255)
                                                                  gtCGACACTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGaaa                              >  2:309160/1‑100 (MQ=255)
                                                                       cacTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAaataat                         <  1:901674/100‑1 (MQ=255)
                                                                       cacTCTTCCCCTTTATTTACATGCTTCTTCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCACATGTGAAaataat                         >  2:1893326/1‑100 (MQ=255)
                                                                        acTCTTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATc                        <  2:2779917/100‑1 (MQ=255)
                                                                         ctctTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCt                       <  2:1259937/100‑1 (MQ=255)
                                                                           ctTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAATAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTAc                     <  2:3734350/100‑1 (MQ=255)
                                                                           ctTCCTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTAc                     >  1:2869033/1‑100 (MQ=255)
                                                                              ccTCTTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTg                  >  1:3821022/1‑100 (MQ=255)
                                                                                tctTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGc                <  1:3833853/100‑1 (MQ=255)
                                                                                tctTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGc                >  2:1744766/1‑100 (MQ=255)
                                                                                 ctTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAGGAGTGTGGCCTGCTGATTTACTTGATATCGCAAGTGAAAGTAATCTACGTGGCg               <  2:3740864/100‑1 (MQ=255)
                                                                                 ctTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGcn               <  1:1688977/100‑2 (MQ=255)
                                                                                 ctTTATTTACATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCg               <  2:1252447/100‑1 (MQ=255)
                                                                                 ctTGATTTACATGCTACTGCCTTTCATTTAAATATGCGGCTGGAAGGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCg               <  2:19189/100‑1 (MQ=255)
                                                                                  tttatttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGt              >  2:2522259/1‑100 (MQ=255)
                                                                                   ttatttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTc             <  2:2414277/100‑1 (MQ=255)
                                                                                     atttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCaa           <  2:2235377/100‑1 (MQ=255)
                                                                                      tttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCaaa          >  2:2228456/1‑100 (MQ=255)
                                                                                       ttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCaaaa         <  1:3942551/100‑1 (MQ=255)
                                                                                       ttaCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCaaaa         <  2:14954/100‑1 (MQ=255)
                                                                                        taCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAAt        <  1:2535779/100‑1 (MQ=255)
                                                                                         aCATGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATc       >  2:1510894/1‑100 (MQ=255)
                                                                                           aTGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCa     <  2:2181205/100‑1 (MQ=255)
                                                                                           aTGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCa     <  1:3466164/100‑1 (MQ=255)
                                                                                           aTGCTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCa     <  1:1138661/100‑1 (MQ=255)
                                                                                           aTGCTTCTGCCTTTCATTTAAAGATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTCGCGTCAAACTCCa     <  2:1235356/100‑1 (MQ=255)
                                                                                              cTTCTGCCTTTCATTTAAATATGCGGCTGGAAAGGGTGTTGCCTGAGGATTTACTTGATATCTCAAGTGAAAATCATCTACGTGGCTTCAAAATCCgtgt  >  2:3782939/1‑100 (MQ=255)
                                                                                                 |                                                                                                
ACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATATCGCAAGTGAAAATAATCTACGTGGCGTCAAAATCCATGT  >  NC_000913/4061138‑4061331

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