Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 298,885 G→A P684S (CCG→TCG)  paoC ← aldehyde dehydrogenase: molybdenum cofactor‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913298,8850GA92.5% 157.3 / ‑3.4 53P684S (CCG→TCG) paoCaldehyde dehydrogenase: molybdenum cofactor‑binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base A (32/17);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 2.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGTCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACC  >  NC_000913/298655‑299126
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
gtggtCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTat                                                                                                                                                                                                                               >  2:401654/1‑251 (MQ=255)
 tggtCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTatc                                                                                                                                                                                                                              >  1:20074/1‑251 (MQ=255)
   gtCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGATGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTatcat                                                                                                                                                                                                                            <  1:25538/251‑1 (MQ=255)
     cATCGGTGAGAAACGGCTGCTCGGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGGATGACATTGGTCCTCCTTAAACCTCATCCGGCAGCTTATCGAGCAGCTTATCCAGATTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGGTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACGCCTTTGG‑CCTTCATCGGGGAGGATATCTGGTCGGTACCATcc                                                                                                                                                                                                                          <  2:76120/251‑1 (MQ=255)
     cATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGACGGTCTCATcc                                                                                                                                                                                                                          >  2:89106/1‑251 (MQ=255)
      aTCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATcc                                                                                                                                                                                                                          >  2:73048/1‑250 (MQ=255)
      aTCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCa                                                                                                                                                                                                                         <  1:3363/251‑1 (MQ=255)
         ggTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAgg                                                                                                                                                                                                                       >  2:309954/1‑250 (MQ=255)
          gttataaaCGGCTGCTCTGGCTTTTGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATCCACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGaa                                                                                                                                                                                                                     <  2:423717/246‑1 (MQ=255)
                   gCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAcctc                                                                                                                                                                                                            >  1:137053/1‑251 (MQ=255)
                       gcgctGGCTTATGCCAGTGTTCGTCTTGGTCAAAAAGCGGGTATGACAGTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCGCGCCGCACAGGCCCAGCTCACCGACACCTTGGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGt                                                                                                                                                                                                        <  2:181344/248‑1 (MQ=255)
                          ctGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGCATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGACAATCACCTCCTGtttt                                                                                                                                                                                                     >  1:70968/1‑251 (MQ=255)
                            ggCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGCATATGGGGTc                                                                                                                                                                                                                                    >  1:482738/1‑218 (MQ=255)
                            ggCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGCATATGGGGTc                                                                                                                                                                                                                                    <  2:482738/218‑1 (MQ=255)
                                   gCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATcc                                                                                                                                                                                            <  2:20074/251‑1 (MQ=255)
                                     cAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGc                                                                                                                                                                                          <  1:216864/251‑1 (MQ=255)
                                      aGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAAGCACCTCCTGTTTTGGGATATCCGCa                                                                                                                                                                                         <  1:178246/251‑1 (MQ=255)
                                                        aaatcGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGTGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATTCACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACATCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAc                                                                                                                                                                       <  1:467070/247‑1 (MQ=255)
                                                         aaGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAcc                                                                                                                                                                      <  1:418793/251‑1 (MQ=255)
                                                          aGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAccc                                                                                                                                                                     <  1:83181/251‑1 (MQ=255)
                                                           gCGGGTATGACATTGTTCCTCCTTAAACCACATACGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAcccc                                                                                                                                                                    >  1:157847/1‑251 (MQ=255)
                                                                 aTGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGCAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCCCCTCATACCCCGccct                                                                                                                                                               >  2:244821/1‑248 (MQ=255)
                                                                  tGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GAGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCatat                                                                                                                                                             >  1:333500/1‑251 (MQ=255)
                                                                      aTTGTTCCTCCTTAAACCACATCCGGCAGCTTCTCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATAGTg                                                                                                                                                         >  1:388137/1‑251 (MQ=255)
                                                                                  tAAACCACATCCGGCAGCTTATCGAGAAGCTTATCCAGAGTGATGGGATCATCCCGTACCCGAATACCGGAGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCCTCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGCTTACCGAAGt                                                                                                                                              >  2:164422/1‑250 (MQ=255)
                                                                                   aaaCCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGt                                                                                                                                              >  1:390642/1‑249 (MQ=255)
                                                                                    aaCCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGc                                                                                                                                           >  1:137364/1‑251 (MQ=255)
                                                                                             cggcggCTTATCGAGAAGTTTACCGAGAGTGGGGGGATACCCCCGTACCCGAATCCCGGGGGCGTTATACACC‑GCGGTGGCGATACCCGCGCTCGCGCCGCACAGGCCCAGCTCACCGCCACCTTTGG‑CCTTCATCGAGGAGGATAGGGGGGCGTTACCACCCAGGAAAACCCCCCCCTGTTTTGGGATAGCCGAATGACCCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGt                                                                                                                                  <  2:157847/251‑1 (MQ=255)
                                                                                               gcagtttctcgggcatcttcgccAGGGTGATGGGATATTCCCTTCCCCGATTACGGGGGGCTTTATACCCC‑GGGTTGGCGAGAGCCGCGCTCCCGCCGCACAGGCCCAGCTCACCGACCCCTTTGG‑CCTTCAGCGAGGAGGATAGGGGGTCGGGATAATCCGGGAAAACCCCCTCCTGTTTTGGGATATCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCa                                                                                                                                <  1:457383/230‑1 (MQ=255)
                                                                                               gCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTCGCCCAAACGGTc                                                                                                                                 >  2:380180/1‑250 (MQ=255)
                                                                                                           gcagcTTATCCAGAGGGATGGGATAATCCCGTCCCCGAATACCGGCGGCATTATACCCC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGACGAGGATATGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCTCCATATCGTGATTAACGACGTCGCCAAAACGGTCATACACCGCCaac                                                                                                                    >  2:509962/1‑249 (MQ=255)
                                                                                                              gtttATCCAGAGTGATGGTATAACCCCGTCCCAGAATCCCGGTGGCGTAATACCCC‑GCGTGGGCGATGGCGGCCCTCCCGCCCCAAAGGCCCAGCCCACCGACCCTTTTGG‑CCTTCACCGAGGAGAATAGGGGGTGGGTATCACCCAGAAAAACCCCCCCCTTTTTGGGGATACCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc                                                                                                                 <  1:169504/249‑1 (MQ=255)
                                                                                                                  actcAGAGGGATGGGAAAACCCCGCCCCCAAATGCCGGTGGGGTTAT‑CACCCGCGTGGGCGATAGCCGCGCTCCCGCCGCACGGCCCCAGCTCCCCGAC‑CCTTTGGACTTTCATGGGGGAGGAGATGGGGTCGTAATCACCAAGGAAAACCCCCCCCTGTTTTGGAGTACCCTAAGGACCCGGCCCCCGAGCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCtcctcc                                                                                                              <  2:34268/247‑1 (MQ=255)
                                                                                                                   tCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCAt                                                                                                            >  1:129860/1‑251 (MQ=255)
                                                                                                                                         gTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGACATCCACCGCCAGCTCCTCCATCAGTGCCGCgcacatccccat                                                                                       >  2:391211/1‑250 (MQ=255)
                                                                                                                                            cccGAATACCGGTGGCGTTATACACC‑GCGTTGGAGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                                                   >  1:453794/1‑251 (MQ=255)
                                                                                                                                                          gCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATCTCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc                                                                     >  2:182168/1‑251 (MQ=255)
                                                                                                                                                            gTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCAACGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCAATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                                                   >  1:206065/1‑251 (MQ=255)
                                                                                                                                                              tataCACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCCATAACCTgg                                                                 >  1:366800/1‑251 (MQ=255)
                                                                                                                                                                 acacC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg                                                              >  1:246721/1‑251 (MQ=255)
                                                                                                                                                                     c‑gcgTTGGCGATACCCGCGCCCACGCCCCACAGCCCCGCCCCACCGACCCCTTTGG‑CCTTCATGGGGGAGGATAGGGGGGCGGTATCACCCAGAAAAACCACCCCCTGTTTGGGGATACCCGAATGACCCGGCACCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcg                                                          <  2:137364/251‑1 (MQ=255)
                                                                                                                                                                        cgTTGGCGATAGCCGCACTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCAAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGCAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCCGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                                                         >  2:181832/1‑250 (MQ=255)
                                                                                                                                                                             gCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAACTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGAGCCGCGCCCATGCCCCTAGTCATTGCGCCAATGCCCTGTCTGCGCGCAGtttt                                                   >  1:504681/1‑251 (MQ=255)
                                                                                                                                                                              cGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc                                                  >  1:509423/1‑251 (MQ=255)
                                                                                                                                                                                   gCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCCATGACCTGGCTGCGCCCAGTTTTCGGAtt                                             >  2:429989/1‑251 (MQ=255)
                                                                                                                                                                                     cgcgCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGACCATAGTCATTGCGCCAATGACCTGGCTGCGAGCAGTTTTCGGATTc                                            >  2:383830/1‑250 (MQ=255)
                                                                                                                                                                                          caccgccgcCCGGGCCCGGCTCACCGAC‑CCTTTGGGCCTTCATCGAGGAGGATAGGGGGTGGGTATCAGCCAGAAAAACCACCCCCTGTTTGGGGATACCCGCATGAACCGGCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg                                      <  1:161719/248‑1 (MQ=255)
                                                                                                                                                                                            ccgccgcACAGGCCCAGCCCCCCGCCCCCTTGGG‑CCTTCATCGAGGAGGATAGGGGGCCGGTATCATCCAGGAAAATCCCCCCCTTTTTTGGGATATCCGCAGAAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCg                                    <  1:401654/250‑1 (MQ=255)
                                                                                                                                                                                                   caggccccAGCCCCCCGACCCCTTGGG‑CCTTCATCGGGGGGGTTATGGGGTCGGTATCCTCCAGAAAAACCCCCTCCTGTTTTGGGATATCCGCATGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgca                             <  1:341256/246‑1 (MQ=255)
                                                                                                                                                                                                    caGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc                                     >  2:89868/1‑242 (MQ=255)
                                                                                                                                                                                                       gCCCAGCTCACCGACACCTTTGG‑CCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac                         >  2:494664/1‑251 (MQ=255)
                                                                                                                                                                                                                 cccgacccTTTGGGCCTTCATCGAGGAGGATATGGGGTCGGTACCACCCAGAAAAATCCCCCCCTGTTTTGGGTTATCCGCATAACCCGCCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGa                <  1:167421/245‑1 (MQ=255)
                                                                                                                                                                                                                               ccTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACAGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGAcc  >  1:366851/1‑251 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
GTGGTCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACC‑GCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACC  >  NC_000913/298655‑299126

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