Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,031,045 T→C 9.1% 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,0450TC 9.1% 135.6 / 3.9 44intergenic (‑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/23);  new base C (2/2);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

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

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