Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 73,833 A→G I230T (ATC→ACC)  yabK ← putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091373,8330AG99.6% 753.0 / NA 234I230T (ATC→ACC) yabKputative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (114/119);  total (114/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

AACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCA  >  NC_000913/73743‑73922
                                                                                          |                                                                                         
aaCAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGtt                                                                                         <  1:4386097/93‑1 (MQ=255)
aaCAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGtt                                                                                         <  1:618965/93‑1 (MQ=255)
aaCAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGtt                                                                                         >  1:2525927/1‑93 (MQ=255)
 aCAGCACCAACCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGAATCGTAGTCGTAACTCAGCGCCTGATAGGTTg                                                                                        <  1:6307524/93‑1 (MQ=255)
   aGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGcc                                                                                      >  2:954732/1‑93 (MQ=255)
   aGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGcc                                                                                      <  1:2749113/93‑1 (MQ=255)
   aGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGcc                                                                                      <  2:372849/93‑1 (MQ=255)
    gCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCGGATAGGATGCCa                                                                                     >  1:469356/1‑93 (MQ=255)
      cccAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGc                                                                                   <  1:692595/92‑1 (MQ=255)
      aCCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGc                                                                                   >  2:2234191/1‑93 (MQ=255)
      aCCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGc                                                                                   >  2:1980930/1‑93 (MQ=255)
       ccAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCt                                                                                  >  1:3069700/1‑93 (MQ=255)
        cAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTc                                                                                 >  1:1989171/1‑93 (MQ=255)
        cAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTc                                                                                 <  2:3236211/93‑1 (MQ=255)
         aGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCg                                                                                <  2:3286362/93‑1 (MQ=255)
         aGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCg                                                                                >  1:4777582/1‑93 (MQ=255)
            ccGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCGGCGCCTGATAGGTTGCCAGCTCGATa                                                                             <  1:5892878/93‑1 (MQ=255)
            ccGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATa                                                                             <  2:3679616/93‑1 (MQ=255)
            ccGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATa                                                                             >  2:5378069/1‑93 (MQ=255)
            ccGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATa                                                                             >  1:321889/1‑93 (MQ=255)
             cGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAg                                                                            <  1:317651/93‑1 (MQ=255)
              gAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGt                                                                           <  1:1484325/93‑1 (MQ=255)
               aGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTg                                                                          <  1:6131198/93‑1 (MQ=255)
                ggCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTgg                                                                         <  1:221547/93‑1 (MQ=255)
                ggCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATACGTTGCCAGCTCGATAGTgg                                                                         >  1:4806032/1‑93 (MQ=255)
                 gcagcaCACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGt                                                                        <  1:1899642/93‑1 (MQ=255)
                  cagcaCACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTc                                                                       >  2:6113739/1‑93 (MQ=255)
                    gcaCACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGc                                                                     >  1:3799477/1‑93 (MQ=255)
                      acacCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCt                                                                   <  1:850565/93‑1 (MQ=255)
                      acacCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCt                                                                   >  1:3994238/1‑93 (MQ=255)
                       cacCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTg                                                                  >  2:5099926/1‑93 (MQ=255)
                       cacCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTg                                                                  >  2:2325701/1‑93 (MQ=255)
                         ccATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCg                                                                <  1:1749026/93‑1 (MQ=255)
                          cATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCgg                                                               >  2:2437964/1‑93 (MQ=255)
                          cATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCgg                                                               >  2:4682332/1‑93 (MQ=255)
                               ggAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGc                                                          >  1:5286431/1‑93 (MQ=255)
                                     gcgcCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCa                                                    <  1:3713947/93‑1 (MQ=255)
                                     gcgcCAGCATTGCCGCGCGGGCAGGATCGAAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCa                                                    >  2:1155130/1‑93 (MQ=255)
                                      cgcCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAg                                                   <  1:467054/93‑1 (MQ=255)
                                       gcCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGc                                                  >  2:4257629/1‑93 (MQ=255)
                                         cAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGa                                                <  2:1289862/93‑1 (MQ=255)
                                           gCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATa                                              >  2:6228413/1‑93 (MQ=255)
                                            cATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAg                                             <  1:3251315/93‑1 (MQ=255)
                                              ttGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCa                                           >  1:983390/1‑93 (MQ=255)
                                              ttGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCa                                           <  2:2554618/93‑1 (MQ=255)
                                              ttGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCa                                           >  2:2752394/1‑93 (MQ=255)
                                              ttGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCa                                           >  2:2923803/1‑93 (MQ=255)
                                              ttGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCa                                           >  1:5046626/1‑93 (MQ=255)
                                                gCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACg                                         <  2:5241447/93‑1 (MQ=255)
                                                  cgcgcgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGt                                       <  1:1715820/93‑1 (MQ=255)
                                                  cgcgcgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGt                                       <  2:822111/93‑1 (MQ=255)
                                                   gcgcgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTg                                      >  2:3819453/1‑93 (MQ=255)
                                                   gcgcgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTg                                      <  1:4268908/93‑1 (MQ=255)
                                                   gcgcgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTg                                      <  2:4289206/93‑1 (MQ=255)
                                                     gcgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGc                                    <  1:4692253/93‑1 (MQ=255)
                                                       gggCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGa                                  <  2:1057907/93‑1 (MQ=255)
                                                       gggCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGa                                  >  1:5644486/1‑93 (MQ=255)
                                                       gggCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGa                                  >  1:3094148/1‑93 (MQ=255)
                                                        ggCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGaa                                 >  1:1267962/1‑93 (MQ=255)
                                                          cAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGc                               >  2:1552996/1‑93 (MQ=255)
                                                           aGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCt                              <  2:2041357/93‑1 (MQ=255)
                                                           aGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCt                              >  2:2832969/1‑93 (MQ=255)
                                                           aGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCt                              >  1:4967776/1‑93 (MQ=255)
                                                            ggATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTg                             >  2:1702451/1‑93 (MQ=255)
                                                            ggATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTg                             <  1:445728/93‑1 (MQ=255)
                                                             gATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTgg                            >  2:5544844/1‑93 (MQ=255)
                                                             gATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTgg                            >  2:3017264/1‑93 (MQ=255)
                                                              aTCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGc                           >  1:3468953/1‑93 (MQ=255)
                                                              aTCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGc                           <  1:1676170/93‑1 (MQ=255)
                                                               tCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCg                          >  2:6273909/1‑93 (MQ=255)
                                                               tCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCg                          >  2:3944859/1‑93 (MQ=255)
                                                                   aGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAAc                      >  2:2181145/1‑93 (MQ=255)
                                                                    gTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACa                     >  1:769783/1‑93 (MQ=255)
                                                                    atCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACa                     <  1:1399645/92‑1 (MQ=255)
                                                                     tCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACag                    >  2:72431/1‑93 (MQ=255)
                                                                     tCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACag                    >  2:2147174/1‑93 (MQ=255)
                                                                       gTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACagag                  >  2:3502795/1‑93 (MQ=255)
                                                                        tAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGc                 <  1:4177203/93‑1 (MQ=255)
                                                                         aaCTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCa                >  1:4523051/1‑93 (MQ=255)
                                                                         aaCTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCa                >  2:1655396/1‑93 (MQ=255)
                                                                          aCTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCAt               <  2:3833601/93‑1 (MQ=255)
                                                                            tCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATaa             >  1:3619965/1‑93 (MQ=255)
                                                                             cAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATaaa            <  1:2116723/93‑1 (MQ=255)
                                                                             cAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCG‑CCCCCAGCGATAGCACCGGGGCGAAGCTGGGGAAACAGAGGATaaaa           >  2:4988490/1‑92 (MQ=255)
                                                                                 gcgtGATAGGTTGCCAGATCGATAGTGGTCGCCTGCGGACCGCCAGCCAGCGATAGCACGTTGGCGAGGCTGGCGAAACAGAGCATAAAGATa        <  2:3996089/90‑1 (MQ=255)
                                                                                 gcCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATa        >  1:2918195/1‑93 (MQ=255)
                                                                                 gcCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATa        <  1:493185/93‑1 (MQ=255)
                                                                                 gcCTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATa        <  1:5192383/93‑1 (MQ=255)
                                                                                  ccTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATaa       <  2:1890421/93‑1 (MQ=255)
                                                                                  ccTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATaa       <  2:5674690/93‑1 (MQ=255)
                                                                                  ccTGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCAAGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATaa       >  1:3917111/1‑93 (MQ=255)
                                                                                   cttaTAGGTTCCAAGCTGTATAGTGGTTGCCTGCAGACCGCTCTCAAGCGATAGCACGGTGGCGAGGCCGGCGAAACAGAGCATAAAGATAAg      <  2:6016312/90‑1 (MQ=255)
                                                                                    tGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTTGCGAAGCTGGCGAAACAGAGCATAAAGATAAga     >  2:4405534/1‑92 (MQ=255)
                                                                                    tGATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAgc     >  1:5781995/1‑93 (MQ=255)
                                                                                     gATAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAgcg    >  1:5929259/1‑93 (MQ=255)
                                                                                      aTAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAgcgc   <  1:249563/93‑1 (MQ=255)
                                                                                       tAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCgca  <  2:5929816/93‑1 (MQ=255)
                                                                                       tAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCgca  >  2:669439/1‑93 (MQ=255)
                                                                                       tAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCgca  >  1:2770178/1‑93 (MQ=255)
                                                                                       tAGGTTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCgca  >  2:3099683/1‑93 (MQ=255)
                                                                                          |                                                                                         
AACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCA  >  NC_000913/73743‑73922

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