Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,149,812 G→T 16.1% V29F (GTC→TTC)  fabD → malonyl‑CoA‑[acyl‑carrier‑protein] transacylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,149,8120GT16.1% 124.6 / 2.6 59V29F (GTC→TTC) fabDmalonyl‑CoA‑[acyl‑carrier‑protein] transacylase
Reads supporting (aligned to +/- strand):  ref base G (11/37);  new base T (8/3);  total (19/40)
Fisher's exact test for biased strand distribution p-value = 2.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTCTGCGGCCTCTGTCCCGTGCGCGCTGGATGAAGCTGTACGCGACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGCAGCAGCAGGGCGGTAAAGCACCGGCAATGATGGCCGGTC‑ACAGCCTGGGGGAATACTCCGCGCTGGTTTGCGCTGGTGTGATTGATTTC  >  NC_000913/1149581‑1150048
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
aCCCCTGCGGCCTCTGTCCCGTGCGCGCTGGATGAAGCTGTACGCGACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCgaaga                                                                                                                                                                                                                                         <  2:262029/239‑1 (MQ=255)
    ctGCGGCCTCTGTCCCGTGCGCGCTGGATGAAGCTGTACGCGACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGGGGGCA‑TTCCAATTTTCCAAGaaaag                                                                                                                                                                                                                                     >  2:1224443/1‑237 (MQ=255)
    ctGCGGCCTCTGTCCCGTGCGCGCTGGATGAAGCTGTACGCGACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCCCGCTGGTTCGTTTCTAGGATAAGGCTTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCGGATATGGCGGGGCGCT‑ATCCAATTGTTGaaaaaaag                                                                                                                                                                                                                                     >  1:1291094/1‑237 (MQ=255)
     tGCGGCCTCTGTCCCGTGCGCGCTGGATGAAGCTGTACGCGACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGGGCT‑ATCCAATTGTCGAAGaaaaga                                                                                                                                                                                                                                    >  2:439097/1‑236 (MQ=255)
      gCGGCCTCTGTCCCGTGCGCGCTGGATGAAGCTGTACGCGACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGGGCT‑TTCCAATTTTCGACGaaaactt                                                                                                                                                                                                                                   >  2:1074299/1‑235 (MQ=255)
      gCGGCCTCTGTCCCGTGCGCGCTGGATGAAGCTGTACGCGACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGGT‑ATTCAATTGTGGAAGAAAAGtt                                                                                                                                                                                                                                   >  2:243744/1‑239 (MQ=255)
               gTCCCGTGCGCGCTGGATGAAGCTGTACGCGACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGc                                                                                                                                                                                                                          <  2:1288535/239‑1 (MQ=255)
                    gTGCGCGCTGGATGAAGCTGTACGCGACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGGTATGGCGGCGAGCG‑ATCCCATTGTGGAAGAAACGGTTTCTGAAGCTTTTg                                                                                                                                                                                                                     >  2:490746/1‑239 (MQ=255)
                     tGCGCGCTGGATGAAGCTGTACGCGACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGGGCT‑ATCCACTTTTCGCAGAAACGGTTTCTTGAGCTTTTgc                                                                                                                                                                                                                    >  2:1041568/1‑239 (MQ=255)
                     tGCGCGCTGGATGAAGCTGTACGCGACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGGC‑CTCCAAATTTTgaagaacacgtgttggaaaatctttc                                                                                                                                                                                                                    >  2:1279291/1‑219 (MQ=255)
                         gggtGGATGAAGCGGTGCCCAGCGGGCGTCTAAACCCGGGTCAGTTGGTGTGGTTTGAAGCCTTTGGCGGTGGATTCCCCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTg                                                                                                                                                                                                                     <  1:968715/231‑1 (MQ=255)
                            tGGATGAAGCTGTACGGGGAGGGGGGATTAAGCGTGGGTGTGTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTTGTTCGGTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCt                                                                                                                                                                                                              <  1:172406/238‑1 (MQ=255)
                            gtggtgaaggtggtccgggcGGGCGGATTAAGCCGGGCGAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTACCCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCCGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCt                                                                                                                                                                                                              <  1:145220/219‑1 (MQ=255)
                                 tgtcgtgtgccggtgtggCGCATTCAGTCTGTGCACTTGGCTCTGGCTGAACCCTGTGGCGGTGTATTCCCCTGGGGCCTCGGGTCTTTGCGTTTCTAGGATAAGGATTCAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAATCCGTTGGAATTCTGGCTGATATGGCGGCGAGCT‑ATCCCATTGTGGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGGTa                                                                                                                                                                                                        <  1:131791/223‑1 (MQ=255)
                                 gaagttgtagtggaggggCGAATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGTATTCACCTGGGGCTCCGCGCTGGTTCTTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTa                                                                                                                                                                                                        <  1:840964/224‑1 (MQ=255)
                                   aGCTGTACGCGACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGGGCT‑ATCCAATTTTTGCAGAAAAGGTTGCTGAAGCTTCTGAGGCCGTggggccgg                                                                                                                                                                                                      >  2:968715/1‑234 (MQ=255)
                                          cgcgcggggcctactAAGGCGGGGCAGTTTGTTCTGCTTGAAGCCTTTGGCGGGTGATTCACCTGGGGCTGCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAAGTTGCATTTGTGTTCCCTGGACAGGGTTCGCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCtgtg                                                                                                                                                                                               <  2:1124458/225‑1 (MQ=255)
                                             gACGGGGTCGTTAAGCCGGGGTTGTTGGTTTGGTTGGAGGCCTTTGGGGTTGGCTTCACCTGGGGCTCCGCGCTGGTTGGGTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGgc                                                                                                                                                                                            <  1:712123/239‑1 (MQ=255)
                                             gACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAAACGTTGGAATGCTGGCTGATATGGCGGCGAGGT‑ATCCAATTTTTGAAGAAACGTTTGTTTGAGCTTCTGAGGCGGTGGGGTATGACCTGTgggg                                                                                                                                                                                            >  2:552682/1‑238 (MQ=255)
                                              aCGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGGGCT‑ATCCAATTTTCGCAGAAACGTTTGCGTAAGCTTCTGCGGCGCCGGGCTACGACCTTTggggg                                                                                                                                                                                           >  2:281088/1‑237 (MQ=255)
                                                   cgtaTTAGGCCGGGGCTGTTGTTTCTGCTGGGAGCCTTTGGGGGTGGATTCACCTGGGGCTCGGCGCTGGTTCGTTCCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGa                                                                                                                                                                                       <  1:281088/235‑1 (MQ=255)
                                                        tgatgcgtggtctgtGGTTCTTCTTGAAGTCTGTGGCTGTGGATTCACGTGGGGACGCGGGTGGGTTCGTTTCAAGGATAAGGATTAAAACATGACGCAATTTGCATGTGTGTTCACTGGACTGGGTTCTCAAACCGGTGGAATGCTGGCTGATATGGCGGTGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCcagc                                                                                                                                                                                 <  1:1022125/225‑1 (MQ=255)
                                                             cGGGGGAGTTGGTTCTGGTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGGTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAggg                                                                                                                                                                             <  2:92904/238‑1 (MQ=255)
                                                                ggCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGcc                                                                                                                                                                          <  2:1226556/238‑1 (MQ=255)
                                                                         tttttttttgagtgttgggCGGGGGATTCACCTGTGGCTCCGCTCTGGTTCGTTTCTAGGATAAGGATTAAAGCAGGACGCAATTTGCATGTGTGTTCCGTGGGCAGGGTTCTCATACCGGTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGACGAAACGTTTGCTGAAGCTTCTGCGGCGCTTGGCTACGACCTTTGGGCTCTGACCCAGAAGGGGCCAGCTgaagaa                                                                                                                                                                <  1:971912/222‑1 (MQ=255)
                                                                          tccgCTTGAAGCCTTTGGCGGTTGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGGTAAGGATTAAATCATGACGCAATTTGCATTTGGGCTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAACCGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGTCCTGTGGGCGCTGACCCAGCAGGGGCCAGCTgaagaa                                                                                                                                                                <  2:278340/235‑1 (MQ=255)
                                                                              cTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCCATTGTTGAAGAAACGTTTGCTGAAGCTTCTGCGGGGCTGGGCTACGACCCGTGGGGGCTGACCCAGCCGGGGGCAGCTGAAGAACTTaa                                                                                                                                                           >  2:1185149/1‑239 (MQ=255)
                                                                                               tGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACt                                                                                                                                           <  2:1100103/238‑1 (MQ=255)
                                                                                                     cACCTGGGGCTCCGCGGTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTg                                                                                                                                    <  2:398399/239‑1 (MQ=255)
                                                                                                           gggCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAAAGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACACAGCAGGGGCCAGCGGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCGGt                                                                                                                              >  2:893027/1‑239 (MQ=255)
                                                                                                                  ggggtGGTTCGTTTCTAAGATAAGGATTAAAACATGGCGCAATTTGCATTTGTGGTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGc                                                                                                                       <  2:129448/235‑1 (MQ=255)
                                                                                                                  ctgttGGTTCGGTTATTGGATAAGGATTAAAACATGACGCAATTTGAATTTGTGTTCCCTGGACAGGGGTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGc                                                                                                                       <  1:1224443/235‑1 (MQ=255)
                                                                                                                    gttGGTTCGTTTCTTAGATAAGGATTAAAACATGACGCAATTTGGATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCAt                                                                                                                     <  2:1140513/237‑1 (MQ=255)
                                                                                                                          tCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGGGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGAGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGtt                                                                                                                <  2:1276847/238‑1 (MQ=255)
                                                                                                                            ttttcgtgtaTAAGTATTAAAACTTGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAATCCGTTGGATTGCGGGCTGTTATGGCGGCGAGCT‑ATCCAATCTTCGAAGAACCGTTTGCTGAAGCTTCTGCGGCGCGGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCTTGCGCTGTTGACTGCATCTGTTgcg                                                                                                             <  1:1279291/230‑1 (MQ=255)
                                                                                                                            gTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTgcg                                                                                                             <  2:24801/239‑1 (MQ=255)
                                                                                                                            agacagaGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCCTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTgcg                                                                                                             <  2:751366/233‑1 (MQ=255)
                                                                                                                              ttCTAGGTGAAGGATTAAAACTTGCCGCAATTTGATTTTGTGTCCCGTGGACAGGGTTCCAAAACCGTTGGAATGCTTGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGGTGGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCTTCTGTTGCGCt                                                                                                           <  1:1074299/239‑1 (MQ=255)
                                                                                                                                    gATAAGGATTAAAACATGCCGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTgcg                                                                                                             <  2:411052/231‑1 (MQ=255)
                                                                                                                                            ttAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGgc                                                                                            <  2:1237163/239‑1 (MQ=255)
                                                                                                                                            ttAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGgc                                                                                            <  2:1186456/239‑1 (MQ=255)
                                                                                                                                                    tgccGCACTTTTCTTTTTTGTTCCCGGGACTGGGTTCTCAAACCGTTTGAAGGCTGGCTGATATGGCGGCGAGCT‑ATCCAAGTGTCGAAGAAACGGTTGCTGAAGCTTCTGAGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAATCCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGcagcagcag                                                                                    <  1:1185149/236‑1 (MQ=255)
                                                                                                                                                    gttcgcgATTTGCTTTTTGTTTTCCTGGACCGGGTTCTCAACGCGTTGGGATGTGGGCTGCGTGTGCGGCGAGCT‑AGCCAATTGTAGAAAGAACGTTGGCCGAAGGTTCTGCGGCGGTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGcagcagcag                                                                                    <  1:1223958/236‑1 (MQ=255)
                                                                                                                                                      ccgacaTTTTCCTTTTTTTTCCCTGGGCCGGGTTCTCAAACCGTTGGAATGGTGGTTGATATGGCGGCGAGCT‑ACCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGTTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGCAGCAGCAggg                                                                                  <  1:552682/234‑1 (MQ=255)
                                                                                                                                                            tttGCATTTGTGTTCCCTGGACAGGGTTCTCACACCGTTGGAATGCTGGCTGATATGGCGGCCAGCT‑ATCAAATTTTCGAAGAAAAGTGTGCTGAAGCTTCTCCGGCCCCCGCCTACCGCCTTGCGGTCCCGTCCCAACCAGGGCCATCTGAAGTACTGAGTTAAACCTCGCACACTCTACCAGCGCGCTGGAGTGCACCTATTGCGCTTT‑GTCCCGCATGgcacgcttcggggggcac                                                                            >  1:210564/1‑224 (MQ=255)
                                                                                                                                                                                    ggTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGATTTCTGAAGCTTCTGAGGCGCTGGGCTGCGACATGTGGGGGCTGACCCCGCAGGGGCCCCCTGAAGAACAGAATAAAACCTCGCCAACTCCGCCCGCGCTGTTGGCTGCATCTGTTTGTCTGTAATCGCG‑AAGGCAGCCACAGGGCGGTAAAGAACCAGCCAAGATGGCGGGTc‑a                                                   >  1:1038844/1‑239 (MQ=255)
                                                                                                                                                                                     tttCTCTAACCGTTGGAATGCTGGTTGGT‑TGGGGGGGGGATGTTCCAATTTTCGAAGAAACCGTTCTAGAAGCTGCTGCAGGCCTGGGCTAAGGACTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGCACTGAATAAAACCTGGCAAACTCGGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGCAGCAGCAGGGCGGTAAAGCACCGGCAATGATGGCCGGTc‑a                                                   <  1:603583/237‑1 (MQ=255)
                                                                                                                                                                                     gTTCTCAAACCGGTGGAATGCTGGCTGATATGGCGGCGAGCT‑TTCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGCAGCAGCAGGGCGGTAAAGCACCGGCAATGATGGCCGGTc‑a                                                   <  2:512885/238‑1 (MQ=255)
                                                                                                                                                                                     gTTCCCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGCAGCAGCAGGGCGGTAAAGCACCGGCAATGATGGCCGGTc‑a                                                   <  2:1291094/238‑1 (MQ=255)
                                                                                                                                                                                         tcAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGGCAGCTGAAGAACTGAATAAAACCTGGCCAACTCAGCCTGCGCTGTTTACTGCATCTGTTGGGCTGG‑ATCGCGGATGGCAGCAGCAGGGCGGGAAAGAACCGGCAATAGTGGTCgggccaagcc                                               >  2:98356/1‑231 (MQ=255)
                                                                                                                                                                                         ccctgccGTTGGAATTGTGGGTGAGTGGGCGGCGATCT‑TTCCAATTGTCGAAGAACCGTTTGCTGAAGGTTCTGGGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGGAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGCAGCAGCAGGGCGGTAAAGCACCGGCAATGATGGCCGGTC‑ACAGcc                                              <  1:439097/234‑1 (MQ=255)
                                                                                                                                                                                                  tGGAATGCTGGCTGATATGGCGGCGAGGT‑ATCCAATTGTCGAAGAAACGGTTGCTGAAGCTTCTGCGGCGCTGGGCTCCGACCCGTGGGCGCTGACCCAGCAGGGGCCCGCTGAAGCACTGCATAAAACCCGGCAACATCAGCCCGCGCCGTTGACTGCATCTGTTGCGCCGT‑ATTGCGGATTGCAGCGGCAGGGCGGTGAAGCAACGGGGACGAGGGCCGGGC‑ACCGCCCgggggggt                                     >  1:743439/1‑236 (MQ=255)
                                                                                                                                                                                                               tgtagggcggtcggCT‑TTCCAATTGTCGAAGAACCGTTTGCTGGAGCTTCGGCGGCGCTGGGCTACGACCTGGGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGCAGCAGCAGGGCGGTAAAGCACCGGCAATGATGGCCGGTC‑ACAGCCTGGGGGAATACTCCGCGCTGGt                        <  1:243744/226‑1 (MQ=255)
                                                                                                                                                                                                               gATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGCAGCAGCAGGGCGGTAAAGCACCGGCAATGATGGCCGGTC‑ACAGCCTGGGGGAATACTCCGCGCTGGt                        <  2:33727/239‑1 (MQ=255)
                                                                                                                                                                                                                atatGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGCAGCAGCAGGGCGGTAAAGCACCGGCAATGATGGCCGGTC‑ACAGCCTGGGGGAATACTCCGCGCTGGt                        <  2:520146/238‑1 (MQ=255)
                                                                                                                                                                                                                     gggggcgggttttCCAATTTTGGGAGAGTCGTGTGCGTAAGCTTCGGGGGGGCTGGGCTTCGACCTGTGGGCGCTGGCCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGGGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGCAGCAGCAGGGCGGTAAAGCACCGGCAATGATGGCCGGTC‑ACAGCCTGGGGGAATACTCCGCGCTGGTTTgcg                   <  1:1041568/227‑1 (MQ=255)
                                                                                                                                                                                                                      cggcgAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGCAGCAGCAGGGCGGTAAAGCACCGGCAATGATGGCCGGTC‑ACAGCCTGGGGGAATACTCCCCGCTGGgtttg                    >  2:983984/1‑234 (MQ=255)
                                                                                                                                                                                                                          gAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGCAGCAGCAGGGCGGTAAAGCACCGGCAATGATGGCCGGTC‑ACAGCCTGGGGGAATACTCCGCGCTGGTTTGCGCTGgt              <  2:792940/238‑1 (MQ=255)
                                                                                                                                                                                                                                      ttGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGCGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCCGCCTGCGCTGTTGACTGCATCTGTTGCGCTGG‑TTCGCGAATGGCGGCAGCAGGGCGGTAAACCACCGGCAATGTTGG‑CGCTCGACAGCCTGGGGTAATACTCCGCGCTGGTTTACGCTGGTGCGATTGGTTTc  >  1:972248/1‑239 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
ACCTCTGCGGCCTCTGTCCCGTGCGCGCTGGATGAAGCTGTACGCGACGGGCGCATTAAGCCGGGGCAGTTGGTTCTGCTTGAAGCCTTTGGCGGTGGATTCACCTGGGGCTCCGCGCTGGTTCGTTTCTAGGATAAGGATTAAAACATGACGCAATTTGCATTTGTGTTCCCTGGACAGGGTTCTCAAACCGTTGGAATGCTGGCTGATATGGCGGCGAGCT‑ATCCAATTGTCGAAGAAACGTTTGCTGAAGCTTCTGCGGCGCTGGGCTACGACCTGTGGGCGCTGACCCAGCAGGGGCCAGCTGAAGAACTGAATAAAACCTGGCAAACTCAGCCTGCGCTGTTGACTGCATCTGTTGCGCTGT‑ATCGCGTATGGCAGCAGCAGGGCGGTAAAGCACCGGCAATGATGGCCGGTC‑ACAGCCTGGGGGAATACTCCGCGCTGGTTTGCGCTGGTGTGATTGATTTC  >  NC_000913/1149581‑1150048

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