Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 955,785 +A coding (71/1761 nt) ycaO ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913955,7821.A98.7% 821.5 / 1.6 225coding (74/1761 nt)ycaOhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base A (102/120);  total (104/121)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

ACAGTGCGCACTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACTTTAG  >  NC_000913/955690‑955869
                                                                                              |                                                                                       
 caGTGCGCACTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGa                                                                                        >  2:2793238/1‑93 (MQ=255)
  aGTGCGCACTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGaa                                                                                       <  2:157104/93‑1 (MQ=255)
   gTGCGCACTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGaaa                                                                                      >  2:6784404/1‑93 (MQ=255)
       gcACTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGaaaagtt                                                                                  >  2:5318173/1‑89 (MQ=255)
         aCTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAgtttt                                                                                <  2:7279115/93‑1 (MQ=255)
         aCTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAgtttt                                                                                >  2:1669439/1‑93 (MQ=255)
          ctctTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGttttt                                                                                <  1:7137986/93‑1 (MQ=255)
           tctTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGGAAAGTTTTTg                                                                               >  2:3229221/1‑93 (MQ=255)
           tctTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTg                                                                               <  1:3507421/93‑1 (MQ=255)
            ctTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGc                                                                              <  2:6188204/93‑1 (MQ=255)
             tttGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCt                                                                             <  1:2194392/93‑1 (MQ=255)
             tttGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCt                                                                             <  2:7597456/93‑1 (MQ=255)
             tttGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCt                                                                             <  1:1202053/93‑1 (MQ=255)
              ttGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTg                                                                            >  1:6531320/1‑93 (MQ=255)
              ttGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTg                                                                            >  2:4659797/1‑93 (MQ=255)
               tGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTgg                                                                           >  2:4741704/1‑93 (MQ=255)
                gTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGa                                                                          >  1:2917572/1‑93 (MQ=255)
                gTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGa                                                                          >  2:2566881/1‑93 (MQ=255)
                  cgcgAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAag                                                                        <  1:7597695/93‑1 (MQ=255)
                   gctaaTATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAagc                                                                       <  1:1286948/90‑1 (MQ=255)
                      aaTATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAagcgag                                                                    >  2:767698/1‑93 (MQ=255)
                       atatGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAagcgagc                                                                   <  1:5437138/93‑1 (MQ=255)
                        tatGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAagcgagcg                                                                  >  1:4092411/1‑93 (MQ=255)
                         atGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAgcgagcga                                                                 <  2:3137234/93‑1 (MQ=255)
                          tGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAgcgagcgag                                                                >  1:3876923/1‑92 (MQ=255)
                          tGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGAt                                                                <  1:188583/93‑1 (MQ=255)
                           gTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATg                                                               <  2:1153992/93‑1 (MQ=255)
                           gTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATg                                                               <  2:6439666/93‑1 (MQ=255)
                           gTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGGTTTTGCTGGAAGCGAGCGCTg                                                               >  2:6666058/1‑93 (MQ=255)
                            tACAGACCAGGCGTTAGGCACGGGCTTCAGACAGGAGGCCCCTTCAATCTGACAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCATCCGATgg                                                              >  2:5401885/1‑93 (MQ=255)
                             aCAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGa                                                             <  2:3491374/93‑1 (MQ=255)
                              cagaccagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGaa                                                            <  1:7777022/93‑1 (MQ=255)
                              cagaccagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGaa                                                            <  2:4866149/93‑1 (MQ=255)
                               agaccagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAAt                                                           <  2:2126129/93‑1 (MQ=255)
                               agaccagacGTTAGGCACGGGATTCAACCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAAt                                                           >  1:4908061/1‑93 (MQ=255)
                                gaccagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATc                                                          <  2:7488099/93‑1 (MQ=255)
                                 accagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCt                                                         <  2:4962944/93‑1 (MQ=255)
                                  ccagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCtt                                                        >  1:5927350/1‑93 (MQ=255)
                                   cagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTc                                                       <  2:2877229/93‑1 (MQ=255)
                                        gTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCagag                                                  <  1:3732688/93‑1 (MQ=255)
                                         ttAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGc                                                 <  1:6298864/93‑1 (MQ=255)
                                         ttAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGc                                                 <  2:7839658/93‑1 (MQ=255)
                                         ttAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGc                                                 >  2:7177790/1‑93 (MQ=255)
                                         ttAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGc                                                 >  2:1053957/1‑93 (MQ=255)
                                         ttAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGCAATCTTCCAGAGc                                                 >  2:3996573/1‑93 (MQ=255)
                                            ggCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCgga                                              >  2:285100/1‑92 (MQ=255)
                                            ggCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGc                                              <  1:2206225/93‑1 (MQ=255)
                                            ggCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGc                                              >  2:231418/1‑93 (MQ=255)
                                            ggCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGc                                              >  1:6793757/1‑93 (MQ=255)
                                             gCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCa                                             <  2:3883637/93‑1 (MQ=255)
                                              cACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCAt                                            >  1:5134792/1‑93 (MQ=255)
                                               aCGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATc                                           >  2:4566140/1‑93 (MQ=255)
                                               aCGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATc                                           <  1:6657542/93‑1 (MQ=255)
                                               aCGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATc                                           <  1:978156/93‑1 (MQ=255)
                                                cGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCt                                          <  1:7276186/93‑1 (MQ=255)
                                                cGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCt                                          >  1:867264/1‑93 (MQ=255)
                                                 gggATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCtt                                         <  2:3601219/93‑1 (MQ=255)
                                                  ggATTCAGCCAGGAGGCCTCTNCNATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCttt                                        <  1:2601633/93‑1 (MQ=255)
                                                   gATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTg                                       <  1:1227705/93‑1 (MQ=255)
                                                   gATTCAGCCAGAAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTg                                       >  2:7800251/1‑93 (MQ=255)
                                                     ttCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGcc                                     <  1:3659700/93‑1 (MQ=255)
                                                     ttCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGcc                                     <  1:4536192/93‑1 (MQ=255)
                                                         gCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCgggg                                 >  1:1775154/1‑93 (MQ=255)
                                                           cAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGAt                               <  1:1068344/93‑1 (MQ=255)
                                                           cAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGAt                               <  1:3103289/93‑1 (MQ=255)
                                                           cAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGAt                               >  2:5908854/1‑93 (MQ=255)
                                                           cAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGAt                               <  1:480541/93‑1 (MQ=255)
                                                              gaggCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATaaa                            <  2:877386/93‑1 (MQ=255)
                                                                 gCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGt                         >  1:7499115/1‑93 (MQ=255)
                                                                   ctctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGttt                       <  2:6399493/93‑1 (MQ=255)
                                                                   ctctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGttt                       <  2:5266557/93‑1 (MQ=255)
                                                                   ctctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGttt                       <  2:7908936/93‑1 (MQ=255)
                                                                    tctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTg                      >  1:7910002/1‑93 (MQ=255)
                                                                    tctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTg                      >  1:5325895/1‑93 (MQ=255)
                                                                    tctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTg                      >  1:6059738/1‑93 (MQ=255)
                                                                     ctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGc                     >  2:951577/1‑93 (MQ=255)
                                                                     ctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGc                     <  1:5787377/93‑1 (MQ=255)
                                                                         aaTCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCa                 <  1:1424675/93‑1 (MQ=255)
                                                                          aTCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCAt                <  1:5170022/93‑1 (MQ=255)
                                                                          aTCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCAt                >  1:3978062/1‑93 (MQ=255)
                                                                           tCTGAAAGCCGGGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATg               >  2:3695105/1‑93 (MQ=255)
                                                                           tCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATg               <  2:4206150/93‑1 (MQ=255)
                                                                           tCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATg               <  1:3632311/93‑1 (MQ=255)
                                                                           tCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATg               >  2:2723346/1‑93 (MQ=255)
                                                                             tGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGaa             <  1:2752765/93‑1 (MQ=255)
                                                                             tGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGaa             >  2:6355719/1‑93 (MQ=255)
                                                                              gAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGaaa            >  1:5310616/1‑93 (MQ=255)
                                                                              gAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGaaa            >  2:4424965/1‑93 (MQ=255)
                                                                               aaaGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGaaaa           <  1:3892218/93‑1 (MQ=255)
                                                                                  gCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCa        <  1:3327511/93‑1 (MQ=255)
                                                                                  gCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCa        <  1:6700581/93‑1 (MQ=255)
                                                                                     gAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACtt     >  2:7714107/1‑93 (MQ=255)
                                                                                     gAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACtt     <  1:3198005/93‑1 (MQ=255)
                                                                                     gAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACtt     >  1:2511237/1‑93 (MQ=255)
                                                                                      aGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACttt    >  1:6097025/1‑93 (MQ=255)
                                                                                       ggTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACTTTa   <  1:4567578/93‑1 (MQ=255)
                                                                                       ggTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACTTTa   <  1:7685237/93‑1 (MQ=255)
                                                                                        gTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACTTTAg  <  1:4539944/93‑1 (MQ=255)
                                                                                        gTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACTTTAg  >  1:3727535/1‑93 (MQ=255)
                                                                                        gTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACTTTAg  >  2:5682734/1‑93 (MQ=255)
                                                                                              |                                                                                       
ACAGTGCGCACTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACTTTAG  >  NC_000913/955690‑955869

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