Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 446,355 T→C N192S (AAC→AGC)  cyoE ← protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913446,3550TC99.7% 1140.0 / ‑5.3 352N192S (AAC→AGC) cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (183/168);  total (184/168)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 363 total aligned reads displayed.

ATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCT  >  NC_000913/446265‑446444
                                                                                          |                                                                                         
 tCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTg                                                                                        >  2:9872646/1‑93 (MQ=255)
 tCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTg                                                                                        >  1:6709751/1‑93 (MQ=255)
 tCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTg                                                                                        >  1:6827363/1‑93 (MQ=255)
  cAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgc                                                                                       <  2:2135032/93‑1 (MQ=255)
     cGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgccgc                                                                                    >  2:7228588/1‑93 (MQ=255)
       tGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCt                                                                                  >  1:6402712/1‑93 (MQ=255)
        ggCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTg                                                                                 <  2:5097405/93‑1 (MQ=255)
        ggCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTg                                                                                 <  1:4297542/93‑1 (MQ=255)
         gCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTgg                                                                                >  1:2380623/1‑93 (MQ=255)
           aaCGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTa                                                                              >  1:8439014/1‑93 (MQ=255)
           aaCGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTa                                                                              >  2:8815518/1‑93 (MQ=255)
             cGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAAt                                                                            >  2:6798826/1‑93 (MQ=255)
             cGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAAt                                                                            >  2:5471285/1‑93 (MQ=255)
             cGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAAt                                                                            <  2:3508678/93‑1 (MQ=255)
             cGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAAt                                                                            <  2:9702718/93‑1 (MQ=255)
                cAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCt                                                                         >  1:1879065/1‑93 (MQ=255)
                 aaaGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCtt                                                                        >  1:205507/1‑93 (MQ=255)
                  aaGGCGATGATATACAGCGTGATGTGATTCTTCTCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTa                                                                       >  2:4263914/1‑93 (MQ=255)
                  aaGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTa                                                                       <  2:5854493/93‑1 (MQ=255)
                  aaGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTa                                                                       <  2:5645757/93‑1 (MQ=255)
                  aaGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTa                                                                       <  1:7805458/93‑1 (MQ=255)
                    ggCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTaaa                                                                     >  1:1137965/1‑93 (MQ=255)
                      cGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGc                                                                   <  1:1641905/93‑1 (MQ=255)
                      cGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGc                                                                   <  1:3850992/93‑1 (MQ=255)
                       gatgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCg                                                                  >  1:2425789/1‑93 (MQ=255)
                         tgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGa                                                                >  2:5791087/1‑93 (MQ=255)
                         tgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGa                                                                >  1:7703093/1‑93 (MQ=255)
                         tgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGa                                                                >  1:3983466/1‑93 (MQ=255)
                          gatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGaa                                                               >  2:171051/1‑93 (MQ=255)
                           atatACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGaaa                                                              <  2:1946479/93‑1 (MQ=255)
                            tataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGaaaa                                                             >  1:9756616/1‑93 (MQ=255)
                            tataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGaaaa                                                             >  2:7586017/1‑93 (MQ=255)
                            tataCAGCGGAATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGaaaa                                                             <  1:1358834/93‑1 (MQ=255)
                            aataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGaaaa                                                             <  2:9581214/92‑1 (MQ=255)
                             ataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAAt                                                            >  1:4780473/1‑93 (MQ=255)
                             ataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAAt                                                            >  2:6214776/1‑93 (MQ=255)
                             ataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAAt                                                            <  2:3605684/93‑1 (MQ=255)
                              taCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATg                                                           >  1:5506820/1‑93 (MQ=255)
                              taCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATg                                                           >  2:7911418/1‑93 (MQ=255)
                               aCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATgg                                                          <  1:795201/93‑1 (MQ=255)
                                  gCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGa                                                       <  2:6967844/93‑1 (MQ=255)
                                  gCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGa                                                       <  1:884100/93‑1 (MQ=255)
                                    gtgatgtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATg                                                     >  2:8579849/1‑93 (MQ=255)
                                      gatgtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGc                                                   >  1:4410400/1‑93 (MQ=255)
                                       atgtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCa                                                  <  1:8311551/93‑1 (MQ=255)
                                        tgttatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCAt                                                 <  1:5245848/93‑1 (MQ=255)
                                        tgtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCAt                                                 >  2:7101307/1‑93 (MQ=255)
                                        tgtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGAGATCGCAt                                                 >  1:3508282/1‑93 (MQ=255)
                                        agtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCAt                                                 <  1:5773105/92‑1 (MQ=255)
                                         gtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATa                                                >  2:2650686/1‑93 (MQ=255)
                                          tgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAg                                               >  2:4223979/1‑93 (MQ=255)
                                           gatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAgg                                              <  2:2292895/93‑1 (MQ=255)
                                            atTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGa                                             >  2:1101290/1‑93 (MQ=255)
                                             ttcttcGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAg                                            <  2:8178265/93‑1 (MQ=255)
                                             ttcttcGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAg                                            <  1:2269223/93‑1 (MQ=255)
                                             ttcttcGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAg                                            <  1:698992/93‑1 (MQ=255)
                                                 tcGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAg                                        <  2:639951/93‑1 (MQ=255)
                                                  cGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAgg                                       <  1:5899485/93‑1 (MQ=255)
                                                   gCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGc                                      <  2:4614113/93‑1 (MQ=255)
                                                   gCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGc                                      >  1:6355739/1‑93 (MQ=255)
                                                     cACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCAt                                    <  1:9176400/93‑1 (MQ=255)
                                                      aCCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATc                                   >  1:1534134/1‑93 (MQ=255)
                                                         gAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGc                                >  1:9448409/1‑93 (MQ=255)
                                                          aaaTGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGcc                               >  2:7726921/1‑93 (MQ=255)
                                                          aaaTGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGcc                               >  2:8448107/1‑93 (MQ=255)
                                                          aaaTGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGcc                               >  2:1392504/1‑93 (MQ=255)
                                                          aaaTGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGcc                               >  1:9869416/1‑93 (MQ=255)
                                                           aaTGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGccc                              >  1:7980414/1‑92 (MQ=255)
                                                            aTGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCcac                             >  1:1008206/1‑93 (MQ=255)
                                                             tGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCcaca                            >  1:1981296/1‑93 (MQ=255)
                                                             tGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCcaca                            >  1:2184914/1‑93 (MQ=255)
                                                              gCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAg                           >  2:6117284/1‑93 (MQ=255)
                                                               ccTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGGGAGGCATCTGCCACAgg                          >  2:934442/1‑93 (MQ=255)
                                                                 ttttACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCt                        <  2:4677403/93‑1 (MQ=255)
                                                                 ttttACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCt                        >  1:6154220/1‑93 (MQ=255)
                                                                 ttttACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCt                        >  2:1027759/1‑93 (MQ=255)
                                                                  tttACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTg                       >  1:5047502/1‑93 (MQ=255)
                                                                   ttACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGa                      <  2:3917116/93‑1 (MQ=255)
                                                                    tACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGaa                     <  2:4861696/93‑1 (MQ=255)
                                                                      ccACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGa                   <  2:5337823/93‑1 (MQ=255)
                                                                      ccACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGa                   <  1:7607050/93‑1 (MQ=255)
                                                                       cACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGAt                  <  1:8147438/93‑1 (MQ=255)
                                                                       cACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGAt                  <  1:2519199/93‑1 (MQ=255)
                                                                       cACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGAt                  <  2:7811596/93‑1 (MQ=255)
                                                                        aCTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATa                 <  1:3618776/93‑1 (MQ=255)
                                                                         cTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAg                <  1:5534537/93‑1 (MQ=255)
                                                                          tGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGc               >  2:6838757/1‑93 (MQ=255)
                                                                          tGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGc               >  1:9035083/1‑93 (MQ=255)
                                                                            gCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCca             >  2:4890119/1‑93 (MQ=255)
                                                                             cAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCcag            <  2:3028073/93‑1 (MQ=255)
                                                                              aaTACCGGAAGGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCcagc           <  1:6630076/93‑1 (MQ=255)
                                                                               aTACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCcagca          <  2:4798923/93‑1 (MQ=255)
                                                                                tACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCcagcag         <  2:2791484/93‑1 (MQ=255)
                                                                                 aCCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAgg        <  1:7173795/93‑1 (MQ=255)
                                                                                  ccGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGa       >  2:4124511/1‑93 (MQ=255)
                                                                                   cGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGAt      <  2:3725518/93‑1 (MQ=255)
                                                                                   cGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGAt      >  1:8649446/1‑93 (MQ=255)
                                                                                    ggAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATc     <  1:3780902/93‑1 (MQ=255)
                                                                                     gAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCg    >  1:3044181/1‑93 (MQ=255)
                                                                                       aTGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCt  >  1:7276429/1‑93 (MQ=255)
                                                                                          |                                                                                         
ATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCT  >  NC_000913/446265‑446444

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