Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 19,014 C→T G100G (GGC→GGT)  nhaR → DNA‑binding transcriptional activator NhaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091319,0140CT96.4% 153.8 / ‑2.1 56G100G (GGC→GGT) nhaRDNA‑binding transcriptional activator NhaR
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (30/24);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-02

TAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGA  >  NC_000913/18768‑19259
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
taAAAAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGAAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAAATGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGACTCACATTTATTGTTTGACATTGGTGTgg                                                                                                                                                                                                                                                   >  2:376078/1‑251 (MQ=255)
           ccGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACtt                                                                                                                                                                                                                                        <  2:89923/251‑1 (MQ=255)
             gTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTc                                                                                                                                                                                                                                      >  1:531824/1‑251 (MQ=255)
               ggTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTcc                                                                                                                                                                                                                                     <  2:319793/250‑1 (MQ=255)
                          cGGAGGCTCTTTAATTAACACCACAAACCATTACAGGACAGATTCGAGCGCTGGAAGAGAGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGGCTCGAACCCAGCGAGCTGGGTGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCTGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCATTTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTc                                                                                                                                                                                                                         <  1:31222/251‑1 (MQ=255)
                           ggAGGCGATTTATTTAACTAGACACACAATTACCGGACAGATTCGAGCGCTGGAAGAGAGCCTGAAAGGCAAATTATTTAAACGCAAGGGCAGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCACAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTc                                                                                                                                                                                                                         >  2:57043/1‑250 (MQ=255)
                                 gcTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCg                                                                                                                                                                                                                  >  1:406224/1‑251 (MQ=255)
                                          aaCTCCACAAACAATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGCCGTGGTCTCGAACCCAGCGAGCTGGGCGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGATTGACGTTGGTGTGGCTGATGAACTTTCCCAACGCCTGGTCAGTAGCGTACTTaccg                                                                                                                                                                                                         >  2:305666/1‑248 (MQ=255)
                                               cacaAACCATTACCGGCAAGATTCGAGCGCTGGAAGACCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGCCTCGAACCAAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCCATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAAcgccgc                                                                                                                                                                                                     >  2:372909/1‑250 (MQ=255)
                                                     ccATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGt                                                                                                                                                                                               >  2:88442/1‑250 (MQ=255)
                                                     ccATTAACGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCATTAGCGTACTTAACGCTGCAGTGGTa                                                                                                                                                                                              <  2:628939/251‑1 (MQ=255)
                                                                    ttCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCAtt                                                                                                                                                                               >  2:172481/1‑251 (MQ=255)
                                                                       gAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCCGTAGCGTACTTAACGCCGCCGTGGTAGAAGGCGAGCCCATTCAt                                                                                                                                                                            >  2:419715/1‑251 (MQ=255)
                                                                        aGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATc                                                                                                                                                                           >  1:465419/1‑251 (MQ=255)
                                                                                        cTGCAAGGCAAATTATTTAAACGCAAGGGACGGGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGATTCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGGAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATc                                                                                                                                                           >  2:535210/1‑251 (MQ=255)
                                                                                          gCAAGGCAAAGTATTTAAACGCAAGGGACGTGGTCGAGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTGCACCTTAAGCCATGAAATGCTGGATCTTGTGAACTGTCGCAAAGAATGCAATTTATTGGTTGACGTTGGTGTTGCTGATGCACTTTCCAAGCGCCTAGGCAGTATCGTACTTGACGCCTCAGTTGTAGAAGGCGAGCCCATTCAACTTCTCTGCTTCGAATc                                                                                                                                                           <  2:413907/249‑1 (MQ=255)
                                                                                            aaGGCAAATTATGTAACCGCAGGGGACGGGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCTCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATGGTGAACTATGGCAAAGAATCCAATTTATTGTTTGACGGTGGTGTGGCTGATGCACTTTCCAAACGCCTTGTCATTAGGGTACTGAACGCCGCAGTGGTACAAGGCGATGCCTTTCATCTTCGCTTCTTAGAAGccacc                                                                                                                                                       <  2:589288/251‑1 (MQ=255)
                                                                                              ggCAAATTATTTAAACGCAATGGACGTGGTCTCGAACACAGCGAGCTGTGAGAATTGATCTATCGCTATGCCGATAGAATGTTCATCGTAAGCTAGGAATTGCTTTATATTGTGAACGTTCGCAAAGTATCCAAATTATTGTTTGACGTTGGTGTGGATGACGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACTCCGCAGTTGTCGAAGGCGATCTCATTCAGCTTTGCTGTTTCGAACcctccca                                                                                                                                                     <  2:110838/251‑1 (MQ=255)
                                                                                                       tttAACCGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAAtgct                                                                                                                                            >  1:133538/1‑251 (MQ=255)
                                                                                                           aaCGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAAATGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGctgctg                                                                                                                                        >  2:625001/1‑251 (MQ=255)
                                                                                                               cAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGc                                                                                                                                    <  1:338256/251‑1 (MQ=255)
                                                                                                                aaGGGACGTGGTCTCGCACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATCAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGACCTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACttt                                                                                                                                                                                                                                       >  2:289084/1‑151 (MQ=255)
                                                                                                                aaGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACttt                                                                                                                                                                                                                                       <  1:289084/151‑1 (MQ=255)
                                                                                                                   ggACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAAtt                                                                                                                                >  1:435427/1‑251 (MQ=255)
                                                                                                                         ggTCCCGAACACAGCGAGCTGTGAGAACTGGACTAGCGCTATGCCGATAAATTGTTCAACTTAAGCCAGGAAATGCTGGATATTGGGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCATTAGCGTCCTTAGCGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCTCTGCTTCGAACCCACCCACGAAATGCTGCTGGAGCAATTAAGTCa                                                                                                                          <  2:346642/251‑1 (MQ=255)
                                                                                                                              cGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGGTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCGTCGCTGCTTCGAATCCACCCACGAAATGCTGCGGGAGCAATTAAGTCAGCa                                                                                                                       <  1:236043/249‑1 (MQ=255)
                                                                                                                                   ccAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATGCATCTTCGCTGCTTCGAATCCACCCACGAAATGCGGCTGGAGCAATTAAGTCAGCATAAACTg                                                                                                                <  2:517334/251‑1 (MQ=255)
                                                                                                                                       cGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTAGCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACGGGata                                                                                                            <  1:575498/251‑1 (MQ=255)
                                                                                                                                       cGAGCTGGGAGAACTGGGCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGata                                                                                                            <  1:380246/251‑1 (MQ=255)
                                                                                                                                                           aTCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGt                                                                                           <  2:133538/248‑1 (MQ=255)
                                                                                                                                                             cGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGAt                                                                                      <  1:187971/251‑1 (MQ=255)
                                                                                                                                                             cGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGAt                                                                                      >  2:562701/1‑251 (MQ=255)
                                                                                                                                                                 aTGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGa                                                                                   <  2:531824/250‑1 (MQ=255)
                                                                                                                                                                    ccGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGActct                                                                               >  2:362830/1‑251 (MQ=255)
                                                                                                                                                                        tAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACgc                                                                           >  1:560924/1‑251 (MQ=255)
                                                                                                                                                                         aaaaTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCAACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACgc                                                                           >  1:82012/1‑250 (MQ=255)
                                                                                                                                                                                 aaCCTTAAGACAGGAAATGCTGGATATTGTGAACTATAGCAAAGCATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGATGATTCGAATCCACCCCCGCAATGCTGCAGGAGCCATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGa                                                                    >  2:380564/2‑249 (MQ=255)
                                                                                                                                                                                  aCCTTACGCCAGGAAATGCTGGATCTTGTGACCTATCGCAAAGAATACAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGCAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGACAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCATGAAgc                                                                 >  2:59559/1‑250 (MQ=255)
                                                                                                                                                                                  aCCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAgg                                                                 >  1:383161/1‑251 (MQ=255)
                                                                                                                                                                                  aCCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAgg                                                                 <  1:24862/251‑1 (MQ=255)
                                                                                                                                                                                        aGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCCGCAGGAAGGCCTGtt                                                           >  2:452877/1‑251 (MQ=255)
                                                                                                                                                                                         gCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTtc                                                          >  1:283459/1‑251 (MQ=255)
                                                                                                                                                                                           cAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGCTGCACTTTCCAAAAGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCAAGCCCATTCCTCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGCGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGCCTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTtctc                                                        >  2:323184/1‑251 (MQ=255)
                                                                                                                                                                                                 aTGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGCGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTg                                                    >  1:416586/1‑249 (MQ=255)
                                                                                                                                                                                                   gCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTTGAGCAAGTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGaa                                                <  1:376078/251‑1 (MQ=255)
                                                                                                                                                                                                     tGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTACGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAAtt                                              >  1:418360/1‑251 (MQ=255)
                                                                                                                                                                                                          atTGTGAACAATCGCAAAGAATCCAATTTATTGTTTGACGATGGTGTGGCTGATGCACTTTCCAAACGCCTGGTAAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCGTAGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCg                                          >  2:225592/1‑250 (MQ=255)
                                                                                                                                                                                                              tgAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAAtgt                                     >  2:97999/1‑251 (MQ=255)
                                                                                                                                                                                                                       gCAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTTGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGctggagt                            <  2:417625/251‑7 (MQ=255)
                                                                                                                                                                                                                        cAAAGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGtt                           <  1:601217/251‑1 (MQ=255)
                                                                                                                                                                                                                           aGAATCCAATTTATTTTTTGCCGTGGGTGTGTCGGATGCTCTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCGTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAATTCAGCATAAACTGGATATGATCATTGCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTTGCGAATGTGGCGTGAGTTTCt                        <  2:568070/251‑1 (MQ=255)
                                                                                                                                                                                                                           aGAATCCAATTTATTGTTTGACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCt                        <  1:260948/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ttGTTTGACGTTGGTGTTGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGCAATGCTGCTGGAGCAATTATGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCCGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATcc           >  2:229886/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gACGTTGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATGTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATccaccacc     <  1:95667/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 gTTGGTGTGGCTGATGCACCTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGTATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGa  <  1:534034/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttGGTGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGa  >  1:439894/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
TAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTACCGGACAGATTCGAGCGCTGGAAGAGCGCCTGCAAGGCAAATTATTTAAACGCAAGGGACGTGGTCTCGAACCCAGCGAGCTGGGAGAACTGGTCTATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCGTACTTAACGCCGCAGTGGTAGAAGGCGAGCCCATTCATCTTCGCTGCTTCGAATCCACCCACGAAATGCTGCTGGAGCAATTAAGTCAGCATAAACTGGATATGATCATTTCTGACTGTCCGATAGACTCTACGCAGCAGGAAGGCCTGTTCTCCGTGAGAATTGGCGAATGTGGCGTGAGTTTCTGGTGTACAAATCCACCACCAGA  >  NC_000913/18768‑19259

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑