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

GTGGTCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACC  >  NC_000913/298655‑299126
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
gtggtCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTat                                                                                                                                                                                                                               >  2:401654/1‑251 (MQ=255)
 tggtCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTatc                                                                                                                                                                                                                              >  1:20074/1‑251 (MQ=255)
   gtCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGATGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTatcat                                                                                                                                                                                                                            <  1:25538/251‑1 (MQ=255)
     cATCGGTGAGAAACGGCTGCTCGGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGGATGACATTGGTCCTCCTTAAACCTCATCCGGCAGCTTATCGAGCAGCTTATCCAGATTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGGTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACGCCTTTGGCCTTCATCGGGGAGGATATCTGGTCGGTACCATcc                                                                                                                                                                                                                          <  2:76120/251‑1 (MQ=255)
     cATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGACGGTCTCATcc                                                                                                                                                                                                                          >  2:89106/1‑251 (MQ=255)
      aTCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATcc                                                                                                                                                                                                                          >  2:73048/1‑250 (MQ=255)
      aTCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCa                                                                                                                                                                                                                         <  1:3363/251‑1 (MQ=255)
         ggTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAgg                                                                                                                                                                                                                       >  2:309954/1‑250 (MQ=255)
          gttataaaCGGCTGCTCTGGCTTTTGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATCCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGaa                                                                                                                                                                                                                     <  2:423717/246‑1 (MQ=255)
                   gCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAcctc                                                                                                                                                                                                            >  1:137053/1‑251 (MQ=255)
                       gcgctGGCTTATGCCAGTGTTCGTCTTGGTCAAAAAGCGGGTATGACAGTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCGCGCCGCACAGGCCCAGCTCACCGACACCTTGGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGt                                                                                                                                                                                                        <  2:181344/248‑1 (MQ=255)
                          ctGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGCATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGACAATCACCTCCTGtttt                                                                                                                                                                                                     >  1:70968/1‑251 (MQ=255)
                            ggCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGCATATGGGGTc                                                                                                                                                                                                                                    >  1:482738/1‑218 (MQ=255)
                            ggCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGCATATGGGGTc                                                                                                                                                                                                                                    <  2:482738/218‑1 (MQ=255)
                                   gCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATcc                                                                                                                                                                                            <  2:20074/251‑1 (MQ=255)
                                     cAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGc                                                                                                                                                                                          <  1:216864/251‑1 (MQ=255)
                                      aGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAAGCACCTCCTGTTTTGGGATATCCGCa                                                                                                                                                                                         <  1:178246/251‑1 (MQ=255)
                                                        aaatcGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGTGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATTCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACATCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAc                                                                                                                                                                       <  1:467070/247‑1 (MQ=255)
                                                         aaGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAcc                                                                                                                                                                      <  1:418793/251‑1 (MQ=255)
                                                          aGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAccc                                                                                                                                                                     <  1:83181/251‑1 (MQ=255)
                                                           gCGGGTATGACATTGTTCCTCCTTAAACCACATACGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAcccc                                                                                                                                                                    >  1:157847/1‑251 (MQ=255)
                                                                 aTGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGCAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCCCCTCATACCCCGccct                                                                                                                                                               >  2:244821/1‑248 (MQ=255)
                                                                  tGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGAGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCatat                                                                                                                                                             >  1:333500/1‑251 (MQ=255)
                                                                      aTTGTTCCTCCTTAAACCACATCCGGCAGCTTCTCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATAGTg                                                                                                                                                         >  1:388137/1‑251 (MQ=255)
                                                                                  tAAACCACATCCGGCAGCTTATCGAGAAGCTTATCCAGAGTGATGGGATCATCCCGTACCCGAATACCGGAGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCCTCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGCTTACCGAAGt                                                                                                                                              >  2:164422/1‑250 (MQ=255)
                                                                                   aaaCCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGt                                                                                                                                              >  1:390642/1‑249 (MQ=255)
                                                                                    aaCCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGc                                                                                                                                           >  1:137364/1‑251 (MQ=255)
                                                                                             cggcggCTTATCGAGAAGTTTACCGAGAGTGGGGGGATACCCCCGTACCCGAATCCCGGGGGCGTTATACACCGCGGTGGCGATACCCGCGCTCGCGCCGCACAGGCCCAGCTCACCGCCACCTTTGGCCTTCATCGAGGAGGATAGGGGGGCGTTACCACCCAGGAAAACCCCCCCCTGTTTTGGGATAGCCGAATGACCCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGt                                                                                                                                  <  2:157847/251‑1 (MQ=255)
                                                                                               gcagtttctcgggcatcttcgccAGGGTGATGGGATATTCCCTTCCCCGATTACGGGGGGCTTTATACCCCGGGTTGGCGAGAGCCGCGCTCCCGCCGCACAGGCCCAGCTCACCGACCCCTTTGGCCTTCAGCGAGGAGGATAGGGGGTCGGGATAATCCGGGAAAACCCCCTCCTGTTTTGGGATATCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCa                                                                                                                                <  1:457383/230‑1 (MQ=255)
                                                                                               gCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTCGCCCAAACGGTc                                                                                                                                 >  2:380180/1‑250 (MQ=255)
                                                                                                           gcagcTTATCCAGAGGGATGGGATAATCCCGTCCCCGAATACCGGCGGCATTATACCCCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGACGAGGATATGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCTCCATATCGTGATTAACGACGTCGCCAAAACGGTCATACACCGCCaac                                                                                                                    >  2:509962/1‑249 (MQ=255)
                                                                                                              gtttATCCAGAGTGATGGTATAACCCCGTCCCAGAATCCCGGTGGCGTAATACCCCGCGTGGGCGATGGCGGCCCTCCCGCCCCAAAGGCCCAGCCCACCGACCCTTTTGGCCTTCACCGAGGAGAATAGGGGGTGGGTATCACCCAGAAAAACCCCCCCCTTTTTGGGGATACCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc                                                                                                                 <  1:169504/249‑1 (MQ=255)
                                                                                                                  actcAGAGGGATGGGAAAACCCCGCCCCCAAATGCCGGTGGGGTTATCACCCGCGTGGGCGATAGCCGCGCTCCCGCCGCACGGCCCCAGCTCCCCGACCCTTTGGACTTTCATGGGGGAGGAGATGGGGTCGTAATCACCAAGGAAAACCCCCCCCTGTTTTGGAGTACCCTAAGGACCCGGCCCCCGAGCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCtcctcc                                                                                                              <  2:34268/247‑1 (MQ=255)
                                                                                                                   tCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCAt                                                                                                            >  1:129860/1‑251 (MQ=255)
                                                                                                                                         gTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGACATCCACCGCCAGCTCCTCCATCAGTGCCGCgcacatccccat                                                                                       >  2:391211/1‑250 (MQ=255)
                                                                                                                                            cccGAATACCGGTGGCGTTATACACCGCGTTGGAGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                                                   >  1:453794/1‑251 (MQ=255)
                                                                                                                                                          gCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATCTCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc                                                                     >  2:182168/1‑251 (MQ=255)
                                                                                                                                                            gTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCAACGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCAATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                                                   >  1:206065/1‑251 (MQ=255)
                                                                                                                                                              tataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCCATAACCTgg                                                                 >  1:366800/1‑251 (MQ=255)
                                                                                                                                                                 acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg                                                              >  1:246721/1‑251 (MQ=255)
                                                                                                                                                                     cgcgTTGGCGATACCCGCGCCCACGCCCCACAGCCCCGCCCCACCGACCCCTTTGGCCTTCATGGGGGAGGATAGGGGGGCGGTATCACCCAGAAAAACCACCCCCTGTTTGGGGATACCCGAATGACCCGGCACCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcg                                                          <  2:137364/251‑1 (MQ=255)
                                                                                                                                                                        cgTTGGCGATAGCCGCACTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCAAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGCAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCCGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                                                         >  2:181832/1‑250 (MQ=255)
                                                                                                                                                                             gCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAACTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGAGCCGCGCCCATGCCCCTAGTCATTGCGCCAATGCCCTGTCTGCGCGCAGtttt                                                   >  1:504681/1‑251 (MQ=255)
                                                                                                                                                                              cGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc                                                  >  1:509423/1‑251 (MQ=255)
                                                                                                                                                                                   gCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCCATGACCTGGCTGCGCCCAGTTTTCGGAtt                                             >  2:429989/1‑251 (MQ=255)
                                                                                                                                                                                     cgcgCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGACCATAGTCATTGCGCCAATGACCTGGCTGCGAGCAGTTTTCGGATTc                                            >  2:383830/1‑250 (MQ=255)
                                                                                                                                                                                          caccgccgcCCGGGCCCGGCTCACCGACCCTTTGGGCCTTCATCGAGGAGGATAGGGGGTGGGTATCAGCCAGAAAAACCACCCCCTGTTTGGGGATACCCGCATGAACCGGCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg                                      <  1:161719/248‑1 (MQ=255)
                                                                                                                                                                                            ccgccgcACAGGCCCAGCCCCCCGCCCCCTTGGGCCTTCATCGAGGAGGATAGGGGGCCGGTATCATCCAGGAAAATCCCCCCCTTTTTTGGGATATCCGCAGAAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCg                                    <  1:401654/250‑1 (MQ=255)
                                                                                                                                                                                                   caggccccAGCCCCCCGACCCCTTGGGCCTTCATCGGGGGGGTTATGGGGTCGGTATCCTCCAGAAAAACCCCCTCCTGTTTTGGGATATCCGCATGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgca                             <  1:341256/246‑1 (MQ=255)
                                                                                                                                                                                                    caGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc                                     >  2:89868/1‑242 (MQ=255)
                                                                                                                                                                                                       gCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac                         >  2:494664/1‑251 (MQ=255)
                                                                                                                                                                                                                 cccgacccTTTGGGCCTTCATCGAGGAGGATATGGGGTCGGTACCACCCAGAAAAATCCCCCCCTGTTTTGGGTTATCCGCATAACCCGCCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGa                <  1:167421/245‑1 (MQ=255)
                                                                                                                                                                                                                               ccTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACAGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGAcc  >  1:366851/1‑251 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
GTGGTCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACC  >  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: