Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,031,044 T→C 8.9% intergenic (‑74/‑116) fadB ← / → pepQ fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase/3‑hydroxyacyl‑CoA dehydrogenase/proline dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,031,0440TC 8.9% 137.2 / 3.9 45intergenic (‑74/‑116)fadB/pepQfused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase/3‑hydroxyacyl‑CoA dehydrogenase/proline dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (17/24);  new base C (2/2);  total (19/26)
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.09e-01

ATCAGCACCGACGATAAAGGCTGCTTTGTTCGAACGCAGCAGCAGCCCTTTTAGATCTGATTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAAC  >  NC_000913/4030764‑4031318
                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                  
atCAGCACCGACGATAAAGGCTGCTTTGTTCGAACGCAGCAGCAGCCCTTTTAGATCTGATTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTGACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGcaaacaa                                                                                                                                                                                                                                                                            <  2:261073/289‑1 (MQ=255)
atCAGCACCGACGATAAAGGCTGCTTTGTTCGAACGCAGCAGCAGCCCTTTTAGATCTGATTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGcaaacaa                                                                                                                                                                                                                                                                            <  1:233807/289‑1 (MQ=255)
atCAGCACCGACGATAAAGGCTGCTTTGTTCGAACGCAGCAGCAGCCCTTTTAGATCTGATTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGcaaacaa                                                                                                                                                                                                                                                                            <  1:256553/289‑1 (MQ=255)
      aCCGACGATAAAGGCTGCTTTGGTCGAACGCAGCAGCAGGCCTTTTAGATCTGATTGCTGTTCCAGCACGTCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGAt                                                                                                                                                                                                                                                                       <  2:485505/288‑1 (MQ=255)
               aaaGGTTGCTTTGTTCAACCGCAGCAGCAGCACTTTTAGATCTGATTGCTGTGCCAGCACGCCGATGGCCCCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATg                                                                                                                                                                                                                                                             <  1:197609/289‑1 (MQ=255)
                     tgctGTGTTCGAACGCAGCAGCAGCCCTTTTAGATCTGATTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGAt                                                                                                                                                                                                                                                        <  1:271967/288‑1 (MQ=255)
                     ggctTTGTTCGGACGCATGCGGAGCCCTTTTAGATCGGGTAGCTATTCCGGCACGCCGATGGCTTCGCGGAGGCTGGCGACGGTCGCGGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAAGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGTTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGAAAACAAAAGATTAAAGAAATGCCGAt                                                                                                                                                                                                                                                        <  1:251910/287‑1 (MQ=255)
                      gctTTGTTCGAACGCAGCAGCAGCCCTTTTAGATCTGATTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGa                                                                                                                                                                                                                                                         >  1:152469/1‑286 (MQ=255)
                                      gcagcagcCCTTTTAGATCTGATTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTACAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAAGATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGa                                                                                                                                                                                                                                      >  1:411043/1‑289 (MQ=255)
                                             cccTTTTAGATCTGATTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCAc                                                                                                                                                                                                                                >  1:139755/1‑288 (MQ=255)
                                               cctttgggtCTGATTGCTGTTCCAGCTCGCCGCGTGCCCGGCCGAGCCTGCCGACGGTCGCAGTGTCGAGTTTAGTAACTGAACCTGGGGAATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGGAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGGTCGGATCCCGCTCGGCAGAAATCACGct                                                                                                                                                                                                                             <  1:119581/281‑1 (MQ=255)
                                                ttttAGATCTGATTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGctc                                                                                                                                                                                                                            <  2:139755/289‑1 (MQ=255)
                                                 tttAGATCTGATTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGc                                                                                                                                                                                                                              >  2:310705/1‑286 (MQ=255)
                                                   tatgtCTGATTGCTGTTACAGCACGCAGATGGCTGCGCCGGGGCTGTCGACGGTCGCAGTGTCGGGTTTATTAACTGAACCTGGGGCATCAAATACCATTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTgg                                                                                                                                                                                                                         <  1:29040/285‑1 (MQ=255)
                                                    aGATCTGATTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGa                                                                                                                                                                                                                        >  1:355639/1‑289 (MQ=255)
                                                          gATTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACg                                                                                                                                                                                                                  <  1:60365/289‑1 (MQ=255)
                                                           aTTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACg                                                                                                                                                                                                                  <  1:490126/288‑1 (MQ=255)
                                                                  ttCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATACGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTa                                                                                                                                                                                                          >  2:248899/1‑289 (MQ=255)
                                                                          cGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCgg                                                                                                                                                                                                  >  2:85813/1‑289 (MQ=255)
                                                                                                  cGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATACGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGa                                                                                                                                                                           <  1:248899/288‑1 (MQ=255)
                                                                                                                             cTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGGAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAAGATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCg                                                                                                                                               <  2:151425/289‑1 (MQ=255)
                                                                                                                             cTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCg                                                                                                                                               <  1:289147/289‑1 (MQ=255)
                                                                                                                                aaCCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATACGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGcc                                                                                                                                                  <  2:81226/283‑1 (MQ=255)
                                                                                                                                aaCCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATACGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGcc                                                                                                                                                  >  1:81226/1‑283 (MQ=255)
                                                                                                                                     ggggCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATaa                                                                                                                                       <  2:152469/289‑1 (MQ=255)
                                                                                                                                      gggCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGa                                                                                                                                                            <  1:341486/267‑1 (MQ=255)
                                                                                                                                      gggCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGa                                                                                                                                                            >  2:341486/1‑267 (MQ=255)
                                                                                                                                       ggCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATaaa                                                                                                                                      <  1:109367/288‑1 (MQ=255)
                                                                                                                                                     gttaggcacacggaacagccagccagTCAAGGTACAGGGTGTCGCCTGTGTAAAGCATGTTCGTCTCCTGAATCCGCAAAGTGATCTAGTCGCAACAGATGAGTCAAAGTGTGTATTTTGTGCTAAAAATATGCAAACAAAAGATTAATGAAATGCCGATCAGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCGCTATAAAAATCATATAGCTAc                                                                                                                        <  1:467615/270‑1 (MQ=255)
                                                                                                                                                       tCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGc                                                                                                                           >  1:296599/1‑283 (MQ=255)
                                                                                                                                                       tCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGc                                                                                                                           <  2:296599/283‑1 (MQ=255)
                                                                                                                                                                                         gggTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAg                                                                                                                                                                                         <  1:405413/187‑1 (MQ=255)
                                                                                                                                                                                         gggTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAg                                                                                                                                                                                         >  2:405413/1‑187 (MQ=255)
                                                                                                                                                                                               cGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTg                                                                             >  2:165099/1‑289 (MQ=255)
                                                                                                                                                                                               cGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGAATCGCGATGCGCTGGAGAGCTTCAAGCTg                                                                             >  1:369883/1‑289 (MQ=255)
                                                                                                                                                                                                               gTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTc                                                             >  2:121305/1‑289 (MQ=255)
                                                                                                                                                                                                               gTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGAGATGAGCTGGCGCGCTTCAAGCTGTCTCCGTTACTTATTc                                                             >  1:42223/1‑289 (MQ=255)
                                                                                                                                                                                                                     agcagggaTCCGCAAAGTGATCTGGTCGTACCAGATTAGGTGAAGTGTTTTTTTTTTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCAAGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACCTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCg                                                       <  1:156340/282‑1 (MQ=255)
                                                                                                                                                                                                                      tcCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCg                                                       >  2:503736/1‑288 (MQ=255)
                                                                                                                                                                                                                                gCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTc                                            <  1:121305/289‑1 (MQ=255)
                                                                                                                                                                                                                                     gTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCCCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCCCTCCGGCGAGCTGtt                                             >  2:156340/1‑283 (MQ=255)
                                                                                                                                                                                                                                                gTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAg                                                  <  1:52578/267‑1 (MQ=255)
                                                                                                                                                                                                                                                gTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAg                                                  >  2:52578/1‑267 (MQ=255)
                                                                                                                                                                                                                                                  aCCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGat                          <  2:301842/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                           tgtgTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAAc  >  2:222906/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                  
ATCAGCACCGACGATAAAGGCTGCTTTGTTCGAACGCAGCAGCAGCCCTTTTAGATCTGATTGCTGTTCCAGCACGCCGATGGCCTCGCCGAGGCTGGCGACGGTCGCAGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAAC  >  NC_000913/4030764‑4031318

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