Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 968,384 T→C S266P (TCT→CCT)  ycaH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913968,3840TC98.2% 857.5 / 1.5 271S266P (TCT→CCT) ycaHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (135/131);  total (137/134)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 274 total aligned reads displayed.

CGGGGATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGTCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCTTT  >  NC_000913/968295‑968473
                                                                                         |                                                                                         
cGGGGATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCtt                                                                                        <  1:4140020/93‑1 (MQ=255)
 ggggATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCttt                                                                                       >  1:1541012/1‑93 (MQ=255)
 ggggATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCttt                                                                                       >  1:9545891/1‑93 (MQ=255)
 ggggATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCttt                                                                                       >  2:4497010/1‑93 (MQ=255)
  gggATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTg                                                                                      <  1:1909149/93‑1 (MQ=255)
  gggATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTg                                                                                      >  1:7359620/1‑93 (MQ=255)
    gATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGaa                                                                                    >  2:8485689/1‑93 (MQ=255)
       tGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCa                                                                                 <  2:9229111/93‑1 (MQ=255)
       tGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCa                                                                                 <  1:8822403/93‑1 (MQ=255)
       tGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCa                                                                                 >  1:834917/1‑93 (MQ=255)
           cgtCCGCCGCGCTGTTTTGCCACGCTGAAGAGGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCg                                                                             <  1:753929/91‑1 (MQ=255)
           cATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCg                                                                             >  2:8501428/1‑93 (MQ=255)
            aTCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCgg                                                                            >  2:7955105/1‑93 (MQ=255)
            aTCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCgg                                                                            <  2:8686246/93‑1 (MQ=255)
             tCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACGACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGa                                                                           >  1:4962149/1‑93 (MQ=255)
             tCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGa                                                                           <  2:2056778/93‑1 (MQ=255)
             tCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGa                                                                           >  1:5915231/1‑93 (MQ=255)
              ccgccgCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGAt                                                                          >  2:8036969/1‑93 (MQ=255)
                 ccggGCTTTTTTGCCACGCTGAAGATGTGTGGCGTAAAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTc                                                                       <  2:8472957/93‑1 (MQ=255)
                 ccgcGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTc                                                                       >  2:2077758/1‑93 (MQ=255)
                 ccgcGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTc                                                                       <  2:1421905/93‑1 (MQ=255)
                   gcgcTTTTTTGCCACGCTGAAGATGTGTGGCGTACAAACGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAg                                                                     >  2:8892465/1‑93 (MQ=255)
                       ttttttGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCg                                                                 <  1:1702539/93‑1 (MQ=255)
                       ttttttGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCg                                                                 >  1:7769480/1‑93 (MQ=255)
                        tttttGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGt                                                                <  1:323619/93‑1 (MQ=255)
                        tttttGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGt                                                                >  1:6888811/1‑93 (MQ=255)
                         ttttGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCGGCCTTTGAACCATGCGGATGTCAGTGCGtt                                                               >  1:6582419/1‑93 (MQ=255)
                         ttttGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGtt                                                               <  2:6226296/93‑1 (MQ=255)
                         ttttGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGtt                                                               >  2:7649208/1‑93 (MQ=255)
                         ttttGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGtt                                                               <  2:3404353/93‑1 (MQ=255)
                            tGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGt                                                            >  2:8996005/1‑93 (MQ=255)
                            tGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGt                                                            >  1:567289/1‑93 (MQ=255)
                              ccACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTaa                                                          >  1:638879/1‑93 (MQ=255)
                                aCGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAgc                                                        >  1:7045638/1‑93 (MQ=255)
                                  gCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAgcgc                                                      >  1:952238/1‑93 (MQ=255)
                                   cTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGcc                                                     <  1:351862/93‑1 (MQ=255)
                                    tGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTCCGCGCCg                                                    >  2:1283042/1‑93 (MQ=255)
                                    tGAAGATGTGCGGCGTACAACCGGAAAACCGTGTACCGCTGGCCGATCATCAGTCTTTGAACCATGCGGATTTCAGCGCGCTGGCACGCGTCg                                                    >  1:3038089/1‑93 (MQ=255)
                                     gAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCgg                                                   <  1:386236/93‑1 (MQ=255)
                                       aGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGc                                                 <  2:7340265/93‑1 (MQ=255)
                                         aTGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCaa                                               >  1:8248836/1‑93 (MQ=255)
                                         aTGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCaa                                               <  2:1603589/93‑1 (MQ=255)
                                           gtgtgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAAc                                             <  1:1149722/93‑1 (MQ=255)
                                           gtgtgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAAc                                             <  1:6230036/93‑1 (MQ=255)
                                              tgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGTCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCt                                          <  2:4716522/93‑1 (MQ=255)
                                              tgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCt                                          <  1:92242/93‑1 (MQ=255)
                                              tgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCt                                          >  2:6546008/1‑93 (MQ=255)
                                              tgGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGAGGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCt                                          >  2:5112273/1‑93 (MQ=255)
                                               ggCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGTCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTg                                         >  1:972756/1‑93 (MQ=255)
                                               ggCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGTCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTg                                         <  2:6438423/93‑1 (MQ=255)
                                               ggCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTg                                         <  2:8197671/93‑1 (MQ=255)
                                               ggCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTg                                         >  1:5222710/1‑93 (MQ=255)
                                               ggCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTg                                         <  2:5265248/93‑1 (MQ=255)
                                                gCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTgg                                        <  2:8271807/93‑1 (MQ=255)
                                                  gTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTa                                      >  1:5589177/1‑93 (MQ=255)
                                                  gTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTa                                      >  2:2974205/1‑93 (MQ=255)
                                                  gTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTa                                      >  1:7679763/1‑93 (MQ=255)
                                                  gTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTa                                      <  1:8427040/93‑1 (MQ=255)
                                                   tACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTaa                                     <  1:2444202/93‑1 (MQ=255)
                                                   tACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTaa                                     <  2:4694939/93‑1 (MQ=255)
                                                    aCAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAAt                                    <  2:5563086/93‑1 (MQ=255)
                                                     cAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATg                                   <  2:7138815/93‑1 (MQ=255)
                                                     cAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATg                                   <  2:5128813/93‑1 (MQ=255)
                                                      aaCCGGAAAAATGTGTACCGCTGGCCTATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGa                                  <  2:2246772/93‑1 (MQ=255)
                                                      aaCCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGa                                  <  1:8557672/93‑1 (MQ=255)
                                                       aCCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGAc                                 <  2:5151368/93‑1 (MQ=255)
                                                       aCCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGAc                                 <  1:8624995/93‑1 (MQ=255)
                                                       aCCGGAAAAATGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGAc                                 <  2:1717810/93‑1 (MQ=255)
                                                             aaaaTGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGaaaa                           >  1:8237984/1‑93 (MQ=255)
                                                                aTGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCATACGCTGGTAATGACTGAACAAGa                        >  2:3189829/1‑93 (MQ=255)
                                                                aTGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGa                        <  2:5875809/93‑1 (MQ=255)
                                                                aTGTGTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGa                        <  1:4811651/93‑1 (MQ=255)
                                                                 tgtgTACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGAt                       >  2:3797180/1‑93 (MQ=255)
                                                                  gtgtACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATg                      >  2:9690593/1‑93 (MQ=255)
                                                                  gtgtACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATg                      <  1:7599570/93‑1 (MQ=255)
                                                                  gtgtACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATg                      <  1:5360029/93‑1 (MQ=255)
                                                                  gtgtACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATg                      <  1:3696670/93‑1 (MQ=255)
                                                                   tgtACCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGc                     <  1:4432895/93‑1 (MQ=255)
                                                                      aCCGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGt                  <  2:5917687/93‑1 (MQ=255)
                                                                        cGCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGa                >  1:9307564/1‑93 (MQ=255)
                                                                         gCTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGaa               <  2:3579496/93‑1 (MQ=255)
                                                                          cTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGaaa              <  1:614679/93‑1 (MQ=255)
                                                                          cTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGaaa              <  2:2641479/93‑1 (MQ=255)
                                                                          cTGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGaaa              <  2:7137244/93‑1 (MQ=255)
                                                                           tGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAAt             >  1:5153338/1‑93 (MQ=255)
                                                                           tGGCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAAt             >  1:9073239/1‑93 (MQ=255)
                                                                            ggCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATg            <  1:4108343/93‑1 (MQ=255)
                                                                            ggCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATg            >  2:4271455/1‑93 (MQ=255)
                                                                             gCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGc           <  2:3871670/93‑1 (MQ=255)
                                                                             gCCGATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGc           >  1:1268725/1‑93 (MQ=255)
                                                                                gATCATCGGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCgg        <  1:8941364/93‑1 (MQ=255)
                                                                                gATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCgg        >  1:801242/1‑93 (MQ=255)
                                                                                gATCATCAGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCgg        >  1:5708539/1‑93 (MQ=255)
                                                                                  tcatcaGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGc      >  2:5640560/1‑93 (MQ=255)
                                                                                  tcatcaGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGc      <  1:1013395/93‑1 (MQ=255)
                                                                                     tcaGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCtt   <  1:7888296/93‑1 (MQ=255)
                                                                                     tcaGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCtt   >  2:6239234/1‑93 (MQ=255)
                                                                                     tcaGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCtt   >  2:4022830/1‑93 (MQ=255)
                                                                                      caGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCttt  <  1:399606/93‑1 (MQ=255)
                                                                                      caGCCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCttt  <  1:2892503/93‑1 (MQ=255)
                                                                                         |                                                                                         
CGGGGATTGGGCATCCGCCGCGCTTTTTTGCCACGCTGAAGATGTGTGGCGTACAACCGGAAAAATGTGTACCGCTGGCCGATCATCAGTCTTTGAACCATGCGGATGTCAGTGCGTTGGTAAGCGCCGGGCAAACGCTGGTAATGACTGAAAAAGATGCGGTGAAATGCCGGGCCTTT  >  NC_000913/968295‑968473

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