Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 311,363 G→A intergenic (‑27/‑749) ecpR ← / → ykgL putative transcriptional regulator for the ecp operon/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913311,3630GA99.2% 384.5 / ‑5.9 119intergenic (‑27/‑749)ecpR/ykgLputative transcriptional regulator for the ecp operon/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (55/63);  total (56/63)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 119 total aligned reads displayed.

ATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGT  >  NC_000913/311229‑311490
                                                                                                                                      |                                                                                                                               
ataGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGc                                                                                                                             <  2:739508/139‑1 (MQ=255)
 taGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCa                                                                                                                            <  1:50927/139‑1 (MQ=255)
  aGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCAt                                                                                                                           <  2:1056616/139‑1 (MQ=255)
      gTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAg                                                                                                                       >  2:197366/1‑139 (MQ=255)
         ggCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGttt                                                                                                                    >  1:577153/1‑139 (MQ=255)
         ggCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGttt                                                                                                                    >  2:920757/1‑139 (MQ=255)
         ggCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGttt                                                                                                                    >  1:887907/1‑139 (MQ=255)
         ggCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGttt                                                                                                                    >  1:226803/1‑139 (MQ=255)
         ggCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGttt                                                                                                                    >  1:1234170/1‑139 (MQ=255)
            gTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCt                                                                                                                 <  1:197366/139‑1 (MQ=255)
              cATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTa                                                                                                               >  1:784098/1‑139 (MQ=255)
              cATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTa                                                                                                               >  2:218187/1‑139 (MQ=255)
                    gagaCCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAAc                                                                                                         <  1:1164710/139‑1 (MQ=255)
                       aCCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCAt                                                                                                                           >  1:797678/1‑118 (MQ=255)
                       aCCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCAt                                                                                                                           <  2:797678/118‑1 (MQ=255)
                       aCCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCgg                                                                                                      >  1:1145831/1‑139 (MQ=255)
                         cAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGaa                                                                                                    <  1:1211380/139‑1 (MQ=255)
                               tATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAAt                                                                                              >  2:1315826/1‑139 (MQ=255)
                                  ttATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           <  1:218187/139‑1 (MQ=255)
                                  ttATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           <  2:1066520/139‑1 (MQ=255)
                                  ttATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           <  1:1185185/139‑1 (MQ=255)
                                   tATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGt                                                                                                                      >  2:331491/1‑111 (MQ=255)
                                   tATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGt                                                                                                                      <  1:331491/111‑1 (MQ=255)
                                    aTCAGAACGGTTTTGACATTCCATATTGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGAt                                                                                         <  2:607715/139‑1 (MQ=255)
                                    aTCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGAt                                                                                         <  2:239934/139‑1 (MQ=255)
                                     tCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATa                                                                                        <  2:724616/139‑1 (MQ=255)
                                      cAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGttt                                                                                                                    >  2:979432/1‑110 (MQ=255)
                                      cAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGttt                                                                                                                    <  1:979432/110‑1 (MQ=255)
                                      cAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATAc                                                                                       >  1:1230436/1‑139 (MQ=255)
                                      cAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATAc                                                                                       >  1:999104/1‑139 (MQ=255)
                                       aGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      <  2:577153/139‑1 (MQ=255)
                                       aGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      <  2:1230436/139‑1 (MQ=255)
                                       aGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      <  2:1234170/139‑1 (MQ=255)
                                         aaCGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      <  1:572514/137‑1 (MQ=255)
                                         aaCGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      >  2:572514/1‑137 (MQ=255)
                                         aaCGGTTTTAACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACTGGAAGTAAATAAGATACgta                                                                                    <  2:999104/139‑1 (MQ=255)
                                          aCGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACgtag                                                                                   <  2:1145831/139‑1 (MQ=255)
                                           cGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACgtagt                                                                                  >  2:1396455/1‑139 (MQ=255)
                                            ggTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACgtagta                                                                                 <  2:705431/139‑1 (MQ=255)
                                            ggTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACgtagta                                                                                 <  1:1396455/139‑1 (MQ=255)
                                                 tGACATTCCATATGGTTTTTAACCTCATAGTCCCTTCTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACatat                                                                            <  2:887907/139‑1 (MQ=255)
                                                   aCATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      <  2:4135/127‑1 (MQ=255)
                                                   aCATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      >  1:4135/1‑127 (MQ=255)
                                                      ttCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATa                                                                       <  2:1204398/139‑1 (MQ=255)
                                                       tCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCgg                                                                                                      >  1:599066/1‑107 (MQ=255)
                                                       tCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCgg                                                                                                      <  2:599066/107‑1 (MQ=255)
                                                        ccATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATaaa                                                                     >  1:354659/1‑139 (MQ=255)
                                                             tGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           <  1:124225/112‑1 (MQ=255)
                                                                  tttAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           >  1:301386/1‑107 (MQ=255)
                                                                  tttAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           <  2:301386/107‑1 (MQ=255)
                                                                   ttAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           <  1:31111/106‑1 (MQ=255)
                                                                   ttAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           >  2:31111/1‑106 (MQ=255)
                                                                   ttAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACgtagta                                                                                 <  2:119460/116‑1 (MQ=255)
                                                                   ttAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTAGTAAACCGGAAGTAAATAAGATACGTAGTACata                                                                             >  2:51986/1‑119 (MQ=255)
                                                                    tAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACgaga                                                         <  2:784098/139‑1 (MQ=255)
                                                                          cATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTc                                                                                                                   <  1:1459240/75‑1 (MQ=255)
                                                                          cATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGa                                                                                                     <  1:1050417/89‑1 (MQ=255)
                                                                          cATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGa                                                                                                     >  2:1050417/1‑89 (MQ=255)
                                                                              gTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTaa                                               <  2:354659/139‑1 (MQ=255)
                                                                                   tACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      <  2:478115/95‑1 (MQ=255)
                                                                                   tACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      >  1:478115/1‑95 (MQ=255)
                                                                                   tACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAattatt                                          >  1:1351615/1‑139 (MQ=255)
                                                                                    aCTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGt                                                                                                             >  2:60755/1‑71 (MQ=255)
                                                                                    aCTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGt                                                                                                             <  1:60755/71‑1 (MQ=255)
                                                                                     cTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           >  1:710809/1‑88 (MQ=255)
                                                                                     cTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           <  2:710809/88‑1 (MQ=255)
                                                                                        tAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGCTTAAGTTTAATTATTTTCGa                                     >  1:1275587/1‑139 (MQ=255)
                                                                                         aaTCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      >  1:848075/1‑89 (MQ=255)
                                                                                         aaTCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGAt                                                                                         >  2:723058/1‑86 (MQ=255)
                                                                                         aaTCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGAt                                                                                         <  1:723058/86‑1 (MQ=255)
                                                                                         aaTCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      <  2:848075/89‑1 (MQ=255)
                                                                                           tCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGc                                                                                                                             <  1:1460755/48‑1 (MQ=255)
                                                                                           tCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGc                                                                                                                             >  2:1460755/1‑48 (MQ=255)
                                                                                           tCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      >  1:1434932/1‑87 (MQ=255)
                                                                                           tCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      <  2:1434932/87‑1 (MQ=255)
                                                                                             aCTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           >  1:1070838/1‑80 (MQ=255)
                                                                                             aCTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           <  2:1070838/80‑1 (MQ=255)
                                                                                               tttGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      >  1:191545/1‑83 (MQ=255)
                                                                                               tttGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                      <  2:191545/83‑1 (MQ=255)
                                                                                               tttGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTaaa                                                              >  1:562370/1‑107 (MQ=255)
                                                                                               tttGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTaaa                                                              <  2:562370/107‑1 (MQ=255)
                                                                                               tttGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCt                              >  2:1416350/1‑139 (MQ=255)
                                                                                                  gCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGc                                                                                                                             >  1:322962/1‑41 (MQ=255)
                                                                                                  gCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGc                                                                                                                             <  2:322962/41‑1 (MQ=255)
                                                                                                  gCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATa                                                                       >  2:1439674/1‑95 (MQ=255)
                                                                                                  gCCATGTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTaa                           >  1:1421158/1‑139 (MQ=255)
                                                                                                       gTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           >  2:495827/1‑70 (MQ=255)
                                                                                                       gTCACTACTTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           <  1:495827/70‑1 (MQ=255)
                                                                                                             actTTCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGtttt                <  1:1416350/139‑1 (MQ=255)
                                                                                                                 tCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           >  1:976414/1‑60 (MQ=255)
                                                                                                                 tCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                           <  2:976414/60‑1 (MQ=255)
                                                                                                                 tCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGtttttttt            >  2:120737/1‑139 (MQ=255)
                                                                                                                 tCCAAACCTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGtttttttt            >  1:995600/1‑139 (MQ=255)
                                                                                                                       ccTGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGAt      <  2:995600/139‑1 (MQ=255)
                                                                                                                         tGTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGa    <  2:1351615/139‑1 (MQ=255)
                                                                                                                          gTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACgtagta                                                                                 >  1:509230/1‑61 (MQ=255)
                                                                                                                          gTAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACgtagta                                                                                 <  2:509230/61‑1 (MQ=255)
                                                                                                                           tAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGaa                                                                                                    >  1:976917/1‑41 (MQ=255)
                                                                                                                           tAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGaa                                                                                                    <  2:976917/41‑1 (MQ=255)
                                                                                                                           tAATTTTTACAAGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAgt  <  1:1060850/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                               
ATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGT  >  NC_000913/311229‑311490

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