Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,881,623 C→T E63K (GAA→AAA)  dmlR ← DNA‑binding transcriptional regulator DmlR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,881,6230CT93.8% 136.6 / 1.3 51E63K (GAA→AAA) dmlRDNA‑binding transcriptional regulator DmlR
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (19/26);  total (20/31)
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.37e-01

TCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAGTTACGGGGATTAATTCCTTTTCAGTCAATTA  >  NC_000913/1881391‑1881864
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
tCATTGACCACATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACgcggc                                                                                                                                                                                                                                   <  1:527600/250‑1 (MQ=255)
     gcccaaagcgactcatCGGTCTTCCACTCCAAAACGTAATTCCAGTTGCGTATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCTAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCTCCTGACGCACGTCATAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGCTGCCGCGGCGggtg                                                                                                                                                                                                                              <  2:224387/235‑1 (MQ=255)
      aCCAAATCGACCAGTCGGTCTTCCCCATCATAACGGAATTCCAGTTGCGGATACGCCTTCGCCAGCGATAATAACGCCGGAGCCCCCCCCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGggtggt                                                                                                                                                                                                                            <  2:523132/251‑1 (MQ=255)
       ccAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCCTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATACATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGgtggtg                                                                                                                                                                                                                           >  1:3026/1‑251 (MQ=255)
             cGACCAGTCGGTCTTCCACATCAAAACGTACTTCCAGTTGCGGATACGCCTTCGCAAGCGCTAATCACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCCCCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCTGCGGTGGGTGCGgtg                                                                                                                                                                                                                      >  1:397052/1‑250 (MQ=255)
              gaccaggaggTCTTACACATAAAAACGTAATGCCAGTTGCGGATACGCCGTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGAGAATACGCAGCATCCCCTGCGGCACCTGACGCACGGCAGAAAGTTCGTCCAGCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCAGGCATAAATCCGCTCTCCTTCTTGGGTAAGGGTGGCGCGGCGGGTGGTGCGgtgga                                                                                                                                                                                                                    <  2:273386/243‑3 (MQ=255)
              gACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGgtgt                                                                                                                                                                                                                     >  2:108497/1‑250 (MQ=255)
               aCCATTCGGTCTTCCCCAGCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCAAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTaa                                                                                                                                                                                                                  <  2:75406/251‑1 (MQ=255)
                   gTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTaacaac                                                                                                                                                                                                              >  2:176169/1‑251 (MQ=255)
                        tCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAACATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGGGACGCGGCGGGTGGTGCGGTGTAACAAcaccac                                                                                                                                                                                                         >  1:354785/1‑251 (MQ=255)
                          ttCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGCCCGACATCCTGCAAACTCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGCCGCGGCGGGTGGGGCGGTGGAACAACACCACGt                                                                                                                                                                                                       >  2:255260/1‑251 (MQ=255)
                                   aaaaCGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACACCACCACGTTTAGCGttt                                                                                                                                                                                              >  1:49404/1‑251 (MQ=255)
                                          aaTTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCa                                                                                                                                                                                        <  1:255260/250‑1 (MQ=255)
                                             tacAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCAGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGAAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTATTGGTGACGCGGCGGGTGTGGCGGTGTTAAGACCACCAGTTTAGTGTTTGCTCCCGCac                                                                                                                                                                                    >  2:237389/3‑250 (MQ=255)
                                              ccAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAg                                                                                                                                                                                   <  1:463294/251‑1 (MQ=255)
                                                 gTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCg                                                                                                                                                                                >  1:372133/1‑251 (MQ=255)
                                                 gTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGGGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCg                                                                                                                                                                                >  1:456130/1‑251 (MQ=255)
                                                    gccgATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCGGCTGATAATACGCCGCATCCCCTGCGGCCCCAGATGCACGTCAGAAAGGTGTTCCATCATCTGACCTACGTCCTGCAAAATCCCCTGCGCCCATGCAGAAATCCGCTCTCCTTCTTTTGTAATGGTGGCTCTGCGTGTGGTGGGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCTATg                                                                                                                                                                             <  1:329885/248‑1 (MQ=255)
                                                       gATACGCATTCGCCAGAGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATAAGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGAAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGGGGTGCGGGGTAACAACACCCAGTTTAGCGTTTGCCCCACCAAGGcctgtagc                                                                                                                                                                          >  2:327248/1‑244 (MQ=255)
                                                            ggcTTCGCCAGCGATAAGATCGCCTGCGCCACCACCTGACGTCCAAAATCAAAGATGCGGAAAAGACGCAGCATCCCCGGCGGCTTCTAACGGACTGCATAAAGGTCCTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGACCAGTCAGAAAGCCGCTCTCCTTCTTTGGTTACGGTGACGCGGCGGGTGGTGCGGTGTAACAACAACACGTTTAGCTTTGGCTCCAGCAAGGCGATGCGCTTGCt                                                                                                                                                                     <  2:509792/249‑1 (MQ=255)
                                                            gCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACgcgg                                                                                                                                                                                                                                    >  1:104077/1‑189 (MQ=255)
                                                            gCCTTCGCCAGAGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGCCATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTACTGGTAACGCGGCGGGTGGTGCGGGGTAACCCCCCCACGTTTAGCGTTTGATCCCGCAAGGCGCTGCGCTTGCt                                                                                                                                                                     >  2:473055/1‑251 (MQ=255)
                                                            gCCTCCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACgcgg                                                                                                                                                                                                                                    <  2:104077/189‑1 (MQ=255)
                                                             ccTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTCTCGTCCACAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGCCATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCTCTCTCCTTCTTCGGTCATGCTGACGCGGCGGGCGGTGCGGTGGAACAACTCCACCTTTATCGTTTGCATCAGCAAGGCGATGCGCTTGCTg                                                                                                                                                                    >  2:222386/1‑251 (MQ=255)
                                                                   ccAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACACAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGGGTAACAACACCCCGTTTAGCGGTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAAc                                                                                                                                                              >  1:504107/1‑251 (MQ=255)
                                                                           aataaCGCCGGAGCCACCACCCGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCAGCCCCTGCGGCCCCTGACTCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGGTTGCTCCAGCAAGGCGATGCGCGTGCTGACGATCGCCGGTGa                                                                                                                                                      <  2:364607/251‑1 (MQ=255)
                                                                                           ccaccTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTcc                                                                                                                                      <  1:176169/251‑1 (MQ=255)
                                                                                                                       gATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAg                                                                                                          <  1:10240/251‑1 (MQ=255)
                                                                                                                          aataCGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGcc                                                                                                        <  1:464730/250‑1 (MQ=255)
                                                                                                                              cGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCGGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAg                                                                                                   <  1:222386/251‑1 (MQ=255)
                                                                                                                                  gcaTCCCCTGCGGCACCTGACGCACGTCAGAAAGTGCGTCCAGCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTTGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACTCCCAGTTCGTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATa                                                                                               <  1:355921/251‑1 (MQ=255)
                                                                                                                                   caTCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATa                                                                                               >  1:397176/1‑250 (MQ=255)
                                                                                                                                       cccTGCGGCCCCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAAAAACACCCCCTTTAGCGTTTGCTCCCGCACGGCGATACGCTTGCTGACGAACGCCGGTTAAACGCCCAGTTCTTCCACCACGGAGGCAAAACCAGCCAGACGAGccacccacctaacgac                                                                                          >  2:567949/1‑240 (MQ=255)
                                                                                                                                           gCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCAGCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTGTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGAcgcg                                                                                       <  1:388915/250‑1 (MQ=255)
                                                                                                                                            cGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCa                                                                                     <  1:473055/251‑1 (MQ=255)
                                                                                                                                               cACCTGACGCACGTCAGAAAGTTCGTCCACCACCTGACCGACATCCGGCAAAATCCGCTGCGCCCATTCATAAAGCCTCTCTCCTCCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAAt                                                                                  <  2:196358/251‑1 (MQ=255)
                                                                                                                                                   tGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGGGCTTGCTGACGAACGCCGGTGAAAAACCCAGTTCTTCCGCCACGGCGGCAAACACAGCCCGGCGAGCCACCCGCATAAAGAcgaccaactaatt                                                                              >  2:331870/1‑240 (MQ=255)
                                                                                                                                                      cGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCGTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTcag                                                                           <  2:49404/251‑1 (MQ=255)
                                                                                                                                                      cGCACGTCAGAAAGTTCGTCCAACATCTGACCGACATCCTGCAAAATCCGCTGCGTCCATTCATAAATCCGCTCTCCTTTTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTcag                                                                           <  1:280234/251‑1 (MQ=255)
                                                                                                                                                          cGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGGCGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCCGCATAAAGACGCGCAAATCATTCagcagc                                                                       <  2:372133/251‑1 (MQ=255)
                                                                                                                                                               gAAAGTTCGTCCATAATCTTACCGACATCCTGAAAAATCCCCTGCGCCCATTCATAAACCCGCCCCCCTTTCTTGGGAACGGTGACTTGGCTGGTGGTGCGGTGTAAAAACACCTCGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCCTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCCCCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTaa                                                                  <  2:460558/251‑1 (MQ=255)
                                                                                                                                                                   gTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCAGAAAGACGCGCAAATCATTCAGCAGCGGTAAatta                                                              <  1:360764/251‑1 (MQ=255)
                                                                                                                                                                    tacgtccctcctatgcccccaaTCCTGCAATACCCGCTGCCCCCATCCATAAATCCGCCCCCCTCCTTTTGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCCGCAAGGCGATGCGCTTGCGTACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAatta                                                              <  2:450320/229‑1 (MQ=255)
                                                                                                                                                                              atcTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAATTCCGCTCGCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCt                                                                                                                                                                     <  2:166800/137‑1 (MQ=255)
                                                                                                                                                                              atcTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCt                                                                                                                                                                     >  1:166800/1‑137 (MQ=255)
                                                                                                                                                                                  gACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCgt                                                  <  1:254959/248‑1 (MQ=255)
                                                                                                                                                                                   cccGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGttt                                              <  1:331870/250‑1 (MQ=255)
                                                                                                                                                                                    ccgtaaTCCTGCACAATCCGCTGCGCCCATGTATAAATCCGCTTTACTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAGCGCCGGGGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCTCCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGtttt                                             <  2:485548/246‑1 (MQ=255)
                                                                                                                                                                                                   agccGCTGCGCCCATTCATAATTCCGCCCCCCTTCTTGGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGACCGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAGTTa                              <  2:299640/249‑1 (MQ=255)
                                                                                                                                                                                                     ccGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAGTTAc                             <  1:108497/250‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTCCGCTCTCCTTCTTTGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAGTTACGGGGATTAATTCCTTTTc           >  1:321469/1‑251 (MQ=255)
                                                                                                                                                                                                                               cgccctctttGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACCACACCCCGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCATCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAGTTACGGGGATTAATTCCTTTTCAGTCAATta  <  2:190301/246‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
TCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCAGCATCCCCTGCGGCACCTGACGCACGTCAGAAAGTTCGTCCATCATCTGACCGACATCCTGCAAAATCCGCTGCGCCCATTCATAAATCCGCTCTCCTTCTTCGGTAATGGTGACGCGGCGGGTGGTGCGGTGTAACAACACCACGTTTAGCGTTTGCTCCAGCAAGGCGATGCGCTTGCTGACGAACGCCGGTGAAACGCCCAGTTCTTCCGCCACGGCGGCAAAACCAGCCCGGCGAGCCACCAGCATAAAGACGCGCAAATCATTCAGCAGCGGTAAATTATTCATGATTCGTGTTTTATGTTTCACCAGTTACGGGGATTAATTCCTTTTCAGTCAATTA  >  NC_000913/1881391‑1881864

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: