Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,838,044 +G coding (498/1458 nt) ascF → PTS system enzyme IIBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,838,0431.G99.3% 893.9 / ‑1.9 284coding (497/1458 nt)ascFPTS system enzyme IIBC
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (145/137);  total (147/137)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 308 total aligned reads displayed.

GATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGA  >  NC_000913/2837951‑2838132
                                                                                             |                                                                                          
gaTGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgc                                                                                             >  1:1536322/1‑93 (MQ=255)
  tGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgcgg                                                                                          >  1:7433067/1‑92 (MQ=255)
  tGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgcgg                                                                                          >  2:5736246/1‑92 (MQ=255)
  tGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgcgg                                                                                          >  2:6843689/1‑92 (MQ=255)
     gcggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgcggcat                                                                                       <  2:7261558/93‑5 (MQ=255)
       ggcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCt                                                                                     <  2:814536/93‑1 (MQ=255)
        gcgTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTg                                                                                    <  1:6780536/93‑1 (MQ=255)
          gTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGcc                                                                                  >  2:1510136/1‑93 (MQ=255)
           tGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcg                                                                                 >  1:2609718/1‑93 (MQ=255)
           tGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcg                                                                                 >  1:1393176/1‑93 (MQ=255)
             cTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgcc                                                                               >  2:6547978/1‑92 (MQ=255)
              tGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgcca                                                                              >  1:5478821/1‑91 (MQ=255)
                aCAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgccatc                                                                            >  2:3037460/1‑89 (MQ=255)
                aCAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCcgccatc                                                                            <  1:1161103/93‑5 (MQ=255)
                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCa                                                                          <  2:5243863/93‑1 (MQ=255)
                 cAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCa                                                                          >  2:6468134/1‑93 (MQ=255)
                  aaaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCaa                                                                         >  1:2014227/1‑93 (MQ=255)
                   aaaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCaaa                                                                        >  2:7706669/1‑93 (MQ=255)
                    aaGGATCGCCGACCTTAACCATTCTGAATGTTATTGGTGACGGTTCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAt                                                                       <  1:1982729/93‑1 (MQ=255)
                    aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAt                                                                       >  2:5635217/1‑93 (MQ=255)
                    aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAt                                                                       >  1:998661/1‑93 (MQ=255)
                    aaGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAt                                                                       <  1:7035886/93‑1 (MQ=255)
                            ccGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAAcc                                                               >  2:2519221/1‑93 (MQ=255)
                            ccGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAAcc                                                               <  1:3214121/93‑1 (MQ=255)
                              gACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCaa                                                             >  1:7631040/1‑93 (MQ=255)
                                ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAata                                                           <  1:4153702/93‑1 (MQ=255)
                                ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAata                                                           >  2:5465033/1‑93 (MQ=255)
                                ccTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAata                                                           <  1:7294130/93‑1 (MQ=255)
                                 cTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAatat                                                          <  1:4338670/93‑1 (MQ=255)
                                 cTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAatat                                                          >  2:4620136/1‑93 (MQ=255)
                                 cTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCGGCCGCCATCAAATTTAAAACCAatat                                                          <  2:1856515/93‑1 (MQ=255)
                                  ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATg                                                         <  1:4245527/93‑1 (MQ=255)
                                  ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATg                                                         >  1:2801438/1‑93 (MQ=255)
                                  ttAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATg                                                         >  2:338414/1‑93 (MQ=255)
                                   tAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGt                                                        <  2:2243634/93‑1 (MQ=255)
                                   tAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGt                                                        <  2:2882166/93‑1 (MQ=255)
                                   tAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGt                                                        <  2:3273202/93‑1 (MQ=255)
                                    aaCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTc                                                       <  1:1349175/93‑1 (MQ=255)
                                     aCCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCg                                                      >  1:2044375/1‑93 (MQ=255)
                                      ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGc                                                     <  2:3900780/93‑1 (MQ=255)
                                      ccATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGc                                                     <  1:6468135/93‑1 (MQ=255)
                                       cATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCt                                                    >  1:5683492/1‑93 (MQ=255)
                                        aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTg                                                   >  2:4472772/1‑93 (MQ=255)
                                        aTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTg                                                   <  2:6903103/93‑1 (MQ=255)
                                         ttCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTgg                                                  >  2:455282/1‑93 (MQ=255)
                                          tCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGc                                                 <  2:4071315/93‑1 (MQ=255)
                                          tCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGc                                                 <  2:4899889/93‑1 (MQ=255)
                                          tCTGAATGGGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGc                                                 <  1:3503975/93‑1 (MQ=255)
                                            tGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGa                                               <  1:3330569/93‑1 (MQ=255)
                                            tGAATGTGATGGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGa                                               <  2:7462596/93‑1 (MQ=255)
                                             gAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGAt                                              <  2:1072905/93‑1 (MQ=255)
                                              aaTGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGAtt                                             >  1:2582036/1‑93 (MQ=255)
                                                   gATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                        <  1:3784256/93‑1 (MQ=255)
                                                   gATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                        <  2:2144348/93‑1 (MQ=255)
                                                   gATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                        <  2:4860040/93‑1 (MQ=255)
                                                   gATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCAt                                        >  1:3593201/1‑93 (MQ=255)
                                                    aTTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATc                                       >  2:1720413/1‑93 (MQ=255)
                                                      tGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATcgc                                     >  1:4707666/1‑93 (MQ=255)
                                                         tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCTGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCggg                                  <  2:6055463/93‑1 (MQ=255)
                                                         tGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCggg                                  <  1:7355054/93‑1 (MQ=255)
                                                          gACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgt                                 <  1:2539987/93‑1 (MQ=255)
                                                          gACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgt                                 <  2:5980527/93‑1 (MQ=255)
                                                            cGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGgtgt                               >  2:1115331/1‑93 (MQ=255)
                                                             ggTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGtgtg                              <  1:5990443/93‑1 (MQ=255)
                                                             ggTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGACATCGCGGGtgtg                              >  1:4535765/1‑93 (MQ=255)
                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                             >  2:6700650/1‑93 (MQ=255)
                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                             >  2:7891932/1‑93 (MQ=255)
                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                             >  1:5315765/1‑93 (MQ=255)
                                                              gTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGc                             >  1:1215474/1‑93 (MQ=255)
                                                               tGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCt                            <  1:7813441/93‑1 (MQ=255)
                                                               tGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCt                            <  1:5078733/93‑1 (MQ=255)
                                                                 cTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTgg                          <  2:6070335/93‑1 (MQ=255)
                                                                   ttcttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTa                        >  2:4356123/1‑93 (MQ=255)
                                                                    tcttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTAc                       <  1:6192994/93‑1 (MQ=255)
                                                                    tctactccCTGCCGCTGATGCTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTAc                       <  2:3744605/93‑1 (MQ=255)
                                                                     cttcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACa                      >  2:4974840/1‑93 (MQ=255)
                                                                       tcttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATc                    >  2:4320668/1‑93 (MQ=255)
                                                                       tctccCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAACCCAATATGTCGCTGCCGATTGCCATCGCGGGTGTGCTGGTACATc                    <  2:7583665/93‑1 (MQ=255)
                                                                        cttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATcc                   >  1:2351194/1‑93 (MQ=255)
                                                                         ttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCg                  >  1:5609273/1‑93 (MQ=255)
                                                                         ttcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGGGTGCTGGTACATCCg                  >  1:683938/1‑93 (MQ=255)
                                                                          tcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa                 <  2:1311064/93‑1 (MQ=255)
                                                                          tcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa                 >  1:2128054/1‑93 (MQ=255)
                                                                          tcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa                 >  1:2938014/1‑93 (MQ=255)
                                                                          tcCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGa                 <  2:1932388/93‑1 (MQ=255)
                                                                            cTGCCGCTGATGGTCGCGTCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATGGCCATCGCGGGTGTGCTGGTACATCCGAGt               <  1:1125031/93‑1 (MQ=255)
                                                                             tGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGtt              >  2:3605742/1‑93 (MQ=255)
                                                                              gCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGttt             >  1:5799421/1‑93 (MQ=255)
                                                                                cGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa           >  1:2479887/1‑93 (MQ=255)
                                                                                cGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa           >  1:7572568/1‑93 (MQ=255)
                                                                                cGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTa           >  1:3993655/1‑93 (MQ=255)
                                                                                 gCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAt          <  2:6243746/93‑1 (MQ=255)
                                                                                 gCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAt          <  1:3131780/93‑1 (MQ=255)
                                                                                  cTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTAtt         >  1:694052/1‑93 (MQ=255)
                                                                                   tGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTg        >  1:4930811/1‑93 (MQ=255)
                                                                                    gATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGa       >  2:2970888/1‑93 (MQ=255)
                                                                                     aTGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGaa      <  2:3538122/93‑1 (MQ=255)
                                                                                      tGGTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAAc     >  1:1017310/1‑93 (MQ=255)
                                                                                       ggTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACt    >  1:1850509/1‑93 (MQ=255)
                                                                                        gTCGCGGCATCTGCCGCCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTg   >  1:1747229/1‑93 (MQ=255)
                                                                                             |                                                                                          
GATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGCATCTGCCGCGCATCAAATTTAAAACCAATATGTCGCTGGCGATTGCCATCGCGGGTGTGCTGGTACATCCGAGTTTTATTGAACTGA  >  NC_000913/2837951‑2838132

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