Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 502,653 C→T intergenic (‑191/+47) ybaL ← / ← fsr putative transporter with NAD(P)‑binding Rossmann‑fold domain/putative fosmidomycin efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913502,6530CT99.1% 2324.5 / NA 645intergenic (‑191/+47)ybaL/fsrputative transporter with NAD(P)‑binding Rossmann‑fold domain/putative fosmidomycin efflux system
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (395/241);  minor base G (6/0);  total (404/241)
Fisher's exact test for biased strand distribution p-value = 8.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 660 total aligned reads displayed.

CATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAG  >  NC_000913/502555‑502750
                                                                                                  |                                                                                                 
  tGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACt                                                                                               <  1:3367197/101‑1 (MQ=255)
    aataatTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGt                                                                                             <  1:185576/101‑1 (MQ=255)
     nnnattattGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTa                                                                                            <  2:7711809/98‑1 (MQ=255)
      taatTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaa                                                                                           >  1:12555429/1‑101 (MQ=255)
       aatTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaaa                                                                                          <  1:5840181/101‑1 (MQ=255)
        attaatGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaaaa                                                                                         >  1:14552825/1‑101 (MQ=255)
         ttattGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAAc                                                                                        >  1:322471/1‑101 (MQ=255)
         ttattGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAAc                                                                                        <  2:2896758/101‑1 (MQ=255)
            ttGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGATTATGGTGACTGTAAAAcnnn                                                                                     >  2:13170099/1‑98 (MQ=255)
            tagAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAAcnnn                                                                                     >  2:9461269/3‑98 (MQ=255)
                aTTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAg                                                                                 <  2:332770/101‑1 (MQ=255)
                    tGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaa                                                                             >  1:6006443/1‑101 (MQ=255)
                     nnngAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaaa                                                                            <  2:13774613/98‑1 (MQ=255)
                       aGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGt                                                                          >  1:2125553/1‑101 (MQ=255)
                        nnnaaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGtt                                                                         <  2:10532897/98‑1 (MQ=255)
                         aaaaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAgnnn                                                                        >  2:14575016/1‑98 (MQ=255)
                           aaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGc                                                                      <  1:5264781/101‑1 (MQ=255)
                            aaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGtttnnn                                                                     >  2:11544155/1‑98 (MQ=255)
                            aaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGcc                                                                     >  1:1487662/1‑101 (MQ=255)
                             aTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCt                                                                    >  1:6708126/1‑101 (MQ=255)
                               ggAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGccnnn                                                                  >  2:8903912/1‑98 (MQ=255)
                                gAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTg                                                                 >  2:1525302/1‑101 (MQ=255)
                                 aaTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGGTTGCCTTTgg                                                                >  2:4386350/1‑101 (MQ=255)
                                  aTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTATAACGATGCAGCCAAAGTTTGCCtttnnn                                                               >  2:10371140/1‑98 (MQ=255)
                                  aTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGc                                                               >  2:1434675/1‑101 (MQ=255)
                                   ttAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTgnnn                                                              >  2:10179158/1‑98 (MQ=255)
                                    tAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTg                                                             >  1:9826256/1‑101 (MQ=255)
                                     aGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGc                                                            >  2:5777854/1‑101 (MQ=255)
                                      nnngCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCa                                                           <  2:14041927/98‑1 (MQ=255)
                                      gCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCa                                                           >  2:5540295/1‑101 (MQ=255)
                                      gCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCa                                                           >  1:4980332/1‑101 (MQ=255)
                                       nnncAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCaa                                                          <  2:12055590/98‑1 (MQ=255)
                                       cTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCaa                                                          >  1:10771711/1‑101 (MQ=255)
                                        tGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCTAAGTTTGCCTTTGGCTGCAAt                                                         >  1:11766616/1‑101 (MQ=255)
                                         gCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATg                                                        >  1:3899527/1‑101 (MQ=255)
                                         gCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATg                                                        >  2:3786348/1‑101 (MQ=255)
                                          cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGa                                                       >  1:4543702/1‑101 (MQ=255)
                                           nnnaaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGaa                                                      <  2:10065253/98‑1 (MQ=255)
                                           aaaaaaaGCACGGATCGGAAAACAGAAGGCGGTTTCACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGaa                                                      <  1:11972824/101‑1 (MQ=255)
                                           aaaaaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGaa                                                      >  1:3668923/1‑101 (MQ=255)
                                                aggCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAAtnnn                                                 >  2:8037285/3‑98 (MQ=255)
                                                aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAg                                                 >  1:8514684/1‑101 (MQ=255)
                                                  gCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTc                                               >  1:4416057/1‑101 (MQ=255)
                                                  gCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCannn                                               >  2:13948457/1‑98 (MQ=255)
                                                   cACTGATCGGAAAACAGAAGGCGTTATAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCANTCt                                              <  2:254632/101‑1 (MQ=255)
                                                    aCGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCtt                                             >  1:2912904/1‑101 (MQ=255)
                                                    aCGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCtt                                             <  1:10026417/101‑1 (MQ=255)
                                                      ggATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTa                                           >  1:14238457/1‑101 (MQ=255)
                                                        aTCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCtttnnn                                         >  2:13432516/1‑98 (MQ=255)
                                                         tCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGc                                        >  1:10996365/1‑101 (MQ=255)
                                                          cGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGct                                       >  1:12221415/1‑100 (MQ=255)
                                                           ggAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATgnnn                                      >  2:9266117/1‑98 (MQ=255)
                                                           ggAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCg                                      >  2:2275858/1‑101 (MQ=255)
                                                           ggAAAACAGAAGGCGGTTTAACAGAGGAAAGGATTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTCGGCTGCAATGAAATCAGTCTTTATGCCg                                      <  1:8002146/101‑1 (MQ=255)
                                                            gAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTGTGcnnn                                     >  2:10702475/1‑98 (MQ=255)
                                                            gAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGcnnn                                     >  2:8328550/1‑98 (MQ=255)
                                                            gAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCgg                                     <  2:6372397/101‑1 (MQ=255)
                                                            gAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCgg                                     >  1:478762/1‑101 (MQ=255)
                                                              aaaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGtt                                   >  2:7109768/1‑101 (MQ=255)
                                                              aaaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGtt                                   <  2:6071365/101‑1 (MQ=255)
                                                               aaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCCATGAAATCAGTCTTTATGCCggnnn                                  >  2:12817904/1‑98 (MQ=255)
                                                               aaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTa                                  <  1:10239665/101‑1 (MQ=255)
                                                                nnngAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTAt                                 <  2:7472366/98‑1 (MQ=255)
                                                                 nnnaaGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTATATGCCGGTTAtt                                <  2:11747277/98‑2 (MQ=255)
                                                                 cAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATc                                <  2:4637330/101‑1 (MQ=255)
                                                                 cAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATc                                >  2:4081031/1‑101 (MQ=255)
                                                                  aGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATca                               <  2:5362331/101‑1 (MQ=255)
                                                                  aGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGCCTGCAATGAAATCAGTCTTTATGCCGGTTATca                               <  1:6903765/101‑1 (MQ=255)
                                                                   gAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcag                              <  1:13692304/101‑1 (MQ=255)
                                                                    aaGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCACAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcagg                             <  2:4795242/101‑1 (MQ=255)
                                                                      ggCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggca                           <  1:13429506/101‑1 (MQ=255)
                                                                       gCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGTCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggcag                          <  1:5630890/101‑1 (MQ=255)
                                                                       gCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggcag                          >  1:10398075/1‑101 (MQ=255)
                                                                        cGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggcagg                         <  1:5777854/101‑1 (MQ=255)
                                                                        cGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggcagg                         >  1:10512656/1‑101 (MQ=255)
                                                                        cGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggcagg                         >  2:669679/1‑101 (MQ=255)
                                                                         nnnttAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGa                        <  2:14461127/98‑1 (MQ=255)
                                                                          gTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGaa                       >  1:11553849/1‑101 (MQ=255)
                                                                           tttAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAat                      <  2:1229810/101‑1 (MQ=255)
                                                                           tttAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAat                      <  1:13053463/101‑1 (MQ=255)
                                                                           tttAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAat                      <  1:14024109/101‑1 (MQ=255)
                                                                           tttAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAat                      <  1:8248016/101‑1 (MQ=255)
                                                                              aaCAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAatnnn                   >  2:7537800/1‑98 (MQ=255)
                                                                               aCAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATgg                  <  1:14790544/101‑1 (MQ=255)
                                                                               aCAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATgg                  <  2:1848566/101‑1 (MQ=255)
                                                                                 agagGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTc                >  2:3107370/1‑101 (MQ=255)
                                                                                 agagGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGGTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTc                >  1:1445339/1‑101 (MQ=255)
                                                                                  gtggAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATggnnn               >  2:10968134/3‑98 (MQ=255)
                                                                                  gagGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCa               >  1:9342757/1‑101 (MQ=255)
                                                                                   agGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCaa              >  1:7280789/1‑101 (MQ=255)
                                                                                      aaaGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACAt           <  1:2705030/101‑1 (MQ=255)
                                                                                        aGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATcc         <  1:4167670/101‑1 (MQ=255)
                                                                                        aGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATcc         <  1:6353651/101‑1 (MQ=255)
                                                                                         ggCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATccc        >  1:12845169/1‑101 (MQ=255)
                                                                                           nnnaTGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGTCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCa      <  2:11161704/98‑1 (MQ=255)
                                                                                              aTGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATTAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATa   >  1:7647298/1‑101 (MQ=255)
                                                                                              aTGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATa   <  1:879430/101‑1 (MQ=255)
                                                                                              aTGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATa   <  2:4980332/101‑1 (MQ=255)
                                                                                              aTGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATa   <  1:2297296/101‑1 (MQ=255)
                                                                                              aTGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATa   >  1:11293342/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
CATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAG  >  NC_000913/502555‑502750

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