Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,814,595 C→T V47M (GTG→ATG)  luxS ← S‑ribosylhomocysteine lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,814,5950CT95.0% 163.2 / 0.3 62V47M (GTG→ATG) luxSS‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (22/35);  total (23/39)
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.02e-01

ACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTC  >  NC_000913/2814350‑2814837
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
acacTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCAtttctt                                                                                                                                                                                                                                               >  1:520779/1‑251 (MQ=255)
acacTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCAtttctt                                                                                                                                                                                                                                               >  2:509575/1‑251 (MQ=255)
  acTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGCGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCCAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTg                                                                                                                                                                                                                                             >  2:193117/1‑251 (MQ=255)
     ggTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGtt                                                                                                                                                                                                                                           <  1:640672/250‑1 (MQ=255)
       tAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCgg                                                                                                                                                                                                                                        >  1:665716/1‑251 (MQ=255)
        aGACGTTCAGTTCCGGGTTCTGATTCTGATCCTGCACTTTCAGCACTTCTTCCATTGCCGCTGTCCAGGCAGCAGCAACACGCTGCCCATCTGGCGTACCAATCAGACTCATATAAAAACCGTTGCCTCAGCCCATTGGCTAGACATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCACCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGc                                                                                                                                                                                                                                       <  1:417296/251‑1 (MQ=255)
         gACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCa                                                                                                                                                                                                                                      <  1:598540/251‑1 (MQ=255)
                 gTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATACAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGa                                                                                                                                                                                                                                >  2:407003/1‑249 (MQ=255)
                    ccGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAgcg                                                                                                                                                                                                                            <  1:318798/250‑1 (MQ=255)
                                 ccgatCCTGCACTTTCATCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCGCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCAACATGTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACg                                                                                                                                                                                                               <  2:662999/248‑1 (MQ=255)
                                   gatCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACg                                                                                                                                                                                                               <  2:172561/248‑1 (MQ=255)
                                      ccTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGAt                                                                                                                                                                                                          <  1:647765/250‑1 (MQ=255)
                                      ccTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGAtt                                                                                                                                                                                                         <  1:186547/251‑1 (MQ=255)
                                              ttCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCg                                                                                                                                                                                                   <  2:554482/249‑1 (MQ=255)
                                                 aGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAAACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGCTGGTTCCGCATAAAACCCGCAAACAGGTGCTCCAGGGTCTGGATCCCTCTTTCTGGCATCATTTCTTTTTTCGGCACGCAGAAGCGCATATAGAACACGGGTCTTGCGTCGCCATg                                                                                                                                                                                              >  2:172147/1‑251 (MQ=255)
                                                   cACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCg                                                                                                                                                                                            <  2:549493/251‑1 (MQ=255)
                                                   cACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCAAGGGTATGGATCTCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCg                                                                                                                                                                                            >  2:501926/1‑251 (MQ=255)
                                                   cACGTCTGTCATTGCCGCTCTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGGGCGGCAGCCCAGTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCGTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCg                                                                                                                                                                                            <  2:59568/251‑1 (MQ=255)
                                                       tCTTCCATTGCCGCTTTCAAGGCATCAGCAACACGCTGCTCATCTGGCGTAACAATCAGCCTCATATAAAAACCGGTGCGGCAGCCCATTGCGAGATATCGATAATCTCTACACCATTACCGATAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCACTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGACCACGGTGATTGCATCGCCCTGCgggg                                                                                                                                                                                         >  2:97771/1‑249 (MQ=255)
                                                                     tttCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTc                                                                                                                                                                           <  2:252104/250‑1 (MQ=255)
                                                                     tttCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTc                                                                                                                                                                           >  1:252104/1‑250 (MQ=255)
                                                                                                    ggCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCg                                                                                                                                              <  1:407003/248‑1 (MQ=255)
                                                                                                    gccggtccAATCAGTCACATATAAAAATCGGTTCGGCAGTCCATGGGCGAGATATCGATAATCGCTACACCATTACCGTTAAGAGTGTTACGCATAAAACCAGCAAACAGTTGCTCCAGGGTATGGATCCCTCTTTCTGGCAGTATTTCTTTGTTCGGCACGCAGAAGCGGAGATCGAACGCGGTGATTTCGTCGCCATGCGGGGGGTTCATTGTTTTCGCCTCCCGATTTGCAGGCTCTTCCATACgg                                                                                                                                             <  2:236708/243‑1 (MQ=255)
                                                                                                              ccatatgcttattaaaaCCGGTGCGGCAGCCCATTGGCGAGATATCGAGAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGAGTGCGGCGCCATGCGGGGTGTTCATTGCTTTCGCCACCCGAACTGCATTCGCTCCCATCCGGTTATGATCGACt                                                                                                                                 <  2:552962/238‑1 (MQ=255)
                                                                                                              acAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACt                                                                                                                                 <  1:314445/250‑1 (MQ=255)
                                                                                                                     cATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGc                                                                                                                          >  2:386491/1‑251 (MQ=255)
                                                                                                                         aaaaaaCCGGTGCGGCAGCTCATTGGCGATATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCAGTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCGTCCATCCGGGTATGATCGACTGTGAAGCTaca                                                                                                                      <  2:116075/250‑3 (MQ=255)
                                                                                                                          aaaaaCCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCg                                                                                                                                                                                                                     >  1:90399/1‑155 (MQ=255)
                                                                                                                          aaaaaCCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCg                                                                                                                                                                                                                     <  2:90399/155‑1 (MQ=255)
                                                                                                                             aaCCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGCTGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATAGACTGTGAAGCTATCTaac                                                                                                                  >  2:85243/1‑251 (MQ=255)
                                                                                                                               ccGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTaacaa                                                                                                                <  1:509575/251‑1 (MQ=255)
                                                                                                                                       gCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCAttt                                                                                                        <  1:68939/251‑1 (MQ=255)
                                                                                                                                       gCAGCCCATTGGCGAGATAGCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATGTCTTTGTTCGGCACGCAGAAGCGCAGATCGATCACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCAttt                                                                                                        <  1:399494/251‑1 (MQ=255)
                                                                                                                                        cAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTa                                                                                                       <  1:386491/251‑1 (MQ=255)
                                                                                                                                          gCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGc                                                                                                     >  1:335330/1‑251 (MQ=255)
                                                                                                                                                    gagaTATCGATAATCTCTACACCATTACCGTTAAGAGGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTGCTGGCATCATTTCTTTGTGCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGt                                                                                           <  1:65406/251‑1 (MQ=255)
                                                                                                                                                     agaTATCGATAAGCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATGGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTa                                                                                          <  1:250879/251‑1 (MQ=255)
                                                                                                                                                        aatCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAAtt                                                                                       <  1:97771/250‑1 (MQ=255)
                                                                                                                                                         atcgtgaacctctACCCCATTACCGTTAAGATGGTGACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCTCCCTTCCTGGCATCATTTCTTTGTGCGTCACGCAGAAGCGCAGTGCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTGCCATCCGGGTATGAGCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCCCCGGTAAtt                                                                                       <  2:679034/241‑1 (MQ=255)
                                                                                                                                                         atCGATAAACGCTACACCATTACCGTTAAGATGGTTTCGCATAAAACCAGCAAGCAGGTTCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGCACACGGTGATTGCGTTGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCTGGTATGATCGACTGTGAAGCTATCTAACAACGGTATTTAGCCACCTCCGGTaagt                                                                                       <  2:469908/250‑3 (MQ=255)
                                                                                                                                                           cGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTTCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACGGTGAAGCTATCTAACAACGTCATTTa                                                                                                       <  2:95139/232‑1 (MQ=255)
                                                                                                                                                           cGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGCCTGTGAAGCTATCTAACAACGTCATTTa                                                                                                       >  1:95139/1‑232 (MQ=255)
                                                                                                                                                                    ctACAACATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGACTCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTACCAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAAtt                                                                           >  2:311101/1‑251 (MQ=255)
                                                                                                                                                                             tACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAAct                                                                  >  1:112075/1‑251 (MQ=255)
                                                                                                                                                                             tACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAAc                                                                   <  1:602361/250‑1 (MQ=255)
                                                                                                                                                                                 tttaatattgTTACTCATAAAACCAGCAAACAGTTGCGCCAGGGTATGGATTCCTCATTTTTGCAGCATTTCTTTGTTCGGCACGCAGAATCGCAGATTGAACACGTTTATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGACACCCGAACTGCAGGTGCTTCCATCCGGGTATGAGCGACTGTGAAGCTATCTAACAACGTCATTTAGCCACCTACGGTAACTTTTTTAATAATTTTTTGAATtctt                                                               <  2:546948/241‑1 (MQ=255)
                                                                                                                                                                                 gTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACtcttt                                                              >  2:677270/1‑251 (MQ=255)
                                                                                                                                                                                         ggTTACGCATAAAACCAGCAAACCGGTGCTCCAGGGTATGGAGCCCTCTTTCTGGCATCATGTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGATATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAg                                                      <  2:208836/251‑1 (MQ=255)
                                                                                                                                                                                         ggTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGTATCCCTCTTTCTGGCATCCGTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGGGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCTCCCGAACTGCAGGCGCTTCCATCCGGTTATGATCGACTGTGAAGCTATCTAACATCGGCATTTAGCCACCTCCGGTAATTTGCTTAAAAATTTTCTGAACTCTTTCTTCCCAg                                                      <  1:472041/251‑1 (MQ=255)
                                                                                                                                                                                          tttACGCATAAAAACCGCCACCAAGTACTCCAGGGTATGGATACCTCGTTCTGGCATCATTTATATGTTCGGCACGCAGAAGCGCAGCTCTAACACGGTGATAGCGTCGCCATGCGGGGTGTTCATTGGTTTCGCCACCCGCACTGCAGGCGCTTCCATCCGGGTATGAGCGAATGTGAAGCTATCTAACAACGGCATTTAGCCACTGACGGTACTTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAgc                                                     >  2:95872/2‑250 (MQ=255)
                                                                                                                                                                                              cGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCg                                                   >  1:11876/1‑249 (MQ=255)
                                                                                                                                                                                                cATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTc                                               >  1:573415/1‑251 (MQ=255)
                                                                                                                                                                                                      aCCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGAGTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCTCCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTTTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTTTGAGta                                         <  2:372041/251‑1 (MQ=255)
                                                                                                                                                                                                            aaaCAGGTGCTCCAGGGGATGGCTCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGGGTCTGAGTATTTGa                                    >  2:68380/1‑250 (MQ=255)
                                                                                                                                                                                                              acatgTGCTCCATGGTATGGATCCCTCTTTCTGTCATCATTTCTTTGGTCGGCACGCAGAAGCGCTGATCGAACACGGTGATTGCTTCGCCATTCGGGTTGTTCATTGTTTGCGCCACCCGAACTGCAGGCGCTTCCATCCGTGTATGATCGACTGTGAAGCTATCTAACATCGGCATTTAGCCCCCCCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGaaac                                 <  2:228182/247‑2 (MQ=255)
                                                                                                                                                                                                                    gCTCCAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTCAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACgcgc                           >  1:586102/1‑251 (MQ=255)
                                                                                                                                                                                                                       ccAGGGTATGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCAt                         <  1:170722/250‑1 (MQ=255)
                                                                                                                                                                                                                              aTGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTatc                 >  1:416358/1‑251 (MQ=255)
                                                                                                                                                                                                                               tGGATCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCg                                                                                                                                    <  1:13376/135‑1 (MQ=255)
                                                                                                                                                                                                                                 gagtccgctTTCTGGCATCATTTCTTTGTTCGACACGCAGAAGCGCAGAGCGAACACGGTGATTGCGGCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCGTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTatcatc              <  2:11876/244‑1 (MQ=255)
                                                                                                                                                                                                                                   tCCCTCTTTCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTatcatcat            >  2:478646/1‑251 (MQ=255)
                                                                                                                                                                                                                                         tttCTGGCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTCATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTg       >  2:632093/1‑250 (MQ=255)
                                                                                                                                                                                                                                              ggCATCATTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAAtt   <  1:632093/249‑1 (MQ=255)
                                                                                                                                                                                                                                              ggCATCATTTCTTTGTGCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAAtt   <  2:573415/249‑1 (MQ=255)
                                                                                                                                                                                                                                                catcatTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTc  >  1:542378/1‑248 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
ACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTC  >  NC_000913/2814350‑2814837

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: