Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,630,551 A→G 11.3% S121S (TCT→TCC ydfI ← putative NAD‑dependent D‑mannonate oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,630,5510AG11.3% 133.4 / 2.3 56S121S (TCT→TCCydfIputative NAD‑dependent D‑mannonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (21/28);  new base G (4/3);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAGGAAGTGACAACGTCACGCATAACATGACCGTTTTCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCATAATATCCCCAGTCACTGAAATGTTCCGTAGCAAGAATA  >  NC_000913/1630316‑1630788
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  
cGTGGGAAGTGACAACCTACCGCATAACTGGACCGTTTTCGGGTATGTTGTCACATGCAATGCAGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCCGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCagaga                                                                                                                                                                                                                                                  <  2:489102/238‑1 (MQ=255)
        gTGACAACGTCACGCATAACATGACCGTTTTCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCTGCAGCTACCATCGGGAGATCGAGCATTAACTGTCCGGTCGCCGGGGAGTGGTAATACCCTTTTTCGGAGATTGTCAGGGAGACAAAcgc                                                                                                                                                                                                                                         >  1:268149/1‑239 (MQ=255)
           aCAACGTCACGCATAACATGACCGTTTTCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGAt                                                                                                                                                                                                                                     >  1:146477/1‑239 (MQ=255)
                   aCGCATAACATGACCGTTTTCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGGGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGt                                                                                                                                                                                                                              >  2:668011/1‑238 (MQ=255)
                          aCATGACCGTTTTCTGGCATGTTGTCACATAACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGCGAGGGGAAATACCCTTTTTCGGGAATTGTCAGAGAGAAAATCGCCATTTGCGGTTccacac                                                                                                                                                                                                                       >  1:254332/1‑233 (MQ=255)
                           cattacttttttCTGGCGTTTTGGCACGTGGCAGGACGGTAAAGTCGGGAAGTCTGGGCGCGTTCCGGGGAGCCAGCGCATCAACAATCACCCCTGTTGGGGTTTTCGGCTGGTGGGGATTTTGCAAGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATc                                                                                                                                                                                                                    <  2:19268/231‑1 (MQ=255)
                               aCCGTTTTCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGGTTGTAAGGAAGACAATCACGATTTGCGGCTcacaccacact                                                                                                                                                                                                                 >  1:376382/1‑232 (MQ=255)
                               aCCGTTTTCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGGGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGGGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTTATTGCAAGAGAGGCAACCGCGATTTGGGGTTcccaactccct                                                                                                                                                                                                                 >  2:47120/1‑231 (MQ=255)
                                     ttCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCCCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACACTCGCGATTTGCGGTTCACACATCGCTGCCAac                                                                                                                                                                                                          >  1:25794/1‑238 (MQ=255)
                                      tCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAacac                                                                                                                                                                                                        <  1:353133/239‑1 (MQ=255)
                                               gtggcacaTTACATGGCGGTAAATGCGGGGAGGCCTGCCGCTTTACGGCGAGCCCGCGCCTCAACAATCCCCCCTGGTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTa                                                                                                                                                                                                <  1:10973/234‑1 (MQ=255)
                                                tttCACATGGCATGACCGTAAATGCGGCAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGGCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAg                                                                                                                                                                                              <  2:481602/237‑1 (MQ=255)
                                                            tgacGGTAAATGGGGGAAGTCCGGTCGCTTTACGGCGAGCCAGCGCCTCAACAATCGCCCCTGTTGCCGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGt                                                                                                                                                                                 <  2:146816/239‑1 (MQ=255)
                                                            tgacGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTTTCAAGGAGACAATCCCGATTTGCGGGTCCCACCCCGCTCCCCAACCGGTTTCTAAGCAATCTATCTGt                                                                                                                                                                                 >  2:687716/1‑239 (MQ=255)
                                                                ggTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTg                                                                                                                                                                                                                                  >  1:10752/1‑189 (MQ=255)
                                                                ggTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTg                                                                                                                                                                                                                                  <  2:10752/189‑1 (MQ=255)
                                                                  gataTGCGGGAAGTTTTGCCGTGTTGCGGCGAGCCAGCCCCTCAAAAATCACCCCTGTTGCTGGTTTCGGCTGGTGGGGATTTTGGACGTCGGCAGATACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGc                                                                                                                                                                           <  2:556731/236‑1 (MQ=255)
                                                                                ccTGCCCGTTTACGGGAGAGCAGCGCCTCAACAATTCCCCCTGGTGCTGGTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAac                                                                                                                                                             <  1:855507/239‑1 (MQ=255)
                                                                                         ttACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGCAATACCCTTTTTCGGTGATTGTCAGAGAGGCAATCGCGATTTGCGGTTCACCACTTCCTGCCAACACGGTTTCTCAGCCCTCTTTTTTTACGGCCAAGGCCTTTTTAaagaccaccaa                                                                                                                                                    >  2:39069/1‑233 (MQ=255)
                                                                                                  aGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGa                                                                                                                                             <  1:775942/238‑1 (MQ=255)
                                                                                                   gCCATCGCCTCAACAATCCGCTCGTTTGCTGTTTTCGGCGTGTGGGGATTTGGCAGGTCGGCAGCTTCTATCGGGTTATCGAGTATTATCTGTCCGGTCTCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTATAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAg                                                                                                                                            <  2:268149/238‑1 (MQ=255)
                                                                                                        cgcCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAAAGACAATCGCGATTTGCGGTTCCTACATCGCTGCCAACACGGTTTCCAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACCCGACCCGTc                                                                                                                                      >  1:630813/1‑238 (MQ=255)
                                                                                                             taaaaaaccacccTGGTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCCTCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCcacaca                                                                                                                                <  1:238093/229‑1 (MQ=255)
                                                                                                                   tCACCCCTGTTGTTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGGTCGAGCATTAACTGTCCGGTCGACGGAGAGTTGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGcc                                                                                                                          <  1:422901/239‑1 (MQ=255)
                                                                                                                   tCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTCCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACACCGCCAACGACCCGAGCCGTCCCCAAATAGGc                                                                                                                           >  1:345991/1‑238 (MQ=255)
                                                                                                                     aCCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGGGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTTTAAGCCCTCTATCTTTACGTGCAAGGCTTTTTTAACGACCCCAACACCGCCAACCGTCCcacacacgcccgc                                                                                                                         >  2:883535/1‑232 (MQ=255)
                                                                                                                       cccTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTTTCAGAGAAACAATGCGGATTTGCGGTTCACCATCGCTTGCCACAACGGTTCTTAGCCCTCTATTCTGTCGGTCCAAGCCTTTTTTAACAAGGCCACGGACCCGAGCGGCCCCCACTCGGcccgcca                                                                                                                     >  1:140902/1‑234 (MQ=255)
                                                                                                                          tGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCCCGATTTGCGGTTCACGCATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGGGCAAGGCTTTTTTAAAGACCCCAACGAACCAGACCGTCCACAAAACGGCCGACaat                                                                                                                    >  2:358727/1‑236 (MQ=255)
                                                                                                                           gTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAAAATGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGTTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCACCGACGCGAGCCGTCACCACATCGGACGCCAttt                                                                                                                   >  1:34481/1‑238 (MQ=255)
                                                                                                                             tGCTGTTTTCGGCTGGGGGGGATTTTGCACGTCGGAAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCg                                                                                                                <  1:890002/239‑1 (MQ=255)
                                                                                                                               gttttttCGGCTGGGTGGGTTTTTGCCCGTCGGCAGCTGCCATCGGGTGGTCGGGCATTAACTGTCCGGTCCCCGGGGGATGGAAATTCCCTTTTTCGGTGAGTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGGCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCa                                                                                                              <  2:443795/236‑1 (MQ=255)
                                                                                                                               cTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACCCGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACACCCGGAGCCGTCCACACCTCGGCCGACCTTTCCGCa                                                                                                              >  1:462686/1‑239 (MQ=255)
                                                                                                                                ttttttCGGCTGGGGTGGGTTTTGCCCGGCGGGAGCCTCCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATGCCCTTTTTCGGGTATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCaa                                                                                                             <  2:621496/237‑1 (MQ=255)
                                                                                                                                tGTTTTCGGCTGGTTGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGGCCGGTCGCCGGGGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCa                                                                                                              <  1:269919/238‑1 (MQ=255)
                                                                                                                                    ttCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCCCACATAGCTGCCAACACGGTTTCTAAGCCATCTATCTTGACGTGCAAGGCTTTTTTACCGGCGCCAACAACACGAGCCGTCCACCCATCCGCCCAACTTTCCGCacgggg                                                                                                         >  2:709287/1‑234 (MQ=255)
                                                                                                                                                  aTTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATTCCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCt                                                                                            <  1:39069/238‑1 (MQ=255)
                                                                                                                                                                     aCCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCCCACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAAt                                                                        >  1:348861/1‑239 (MQ=255)
                                                                                                                                                                                  ctggCATTAACTGTCCGGGCGCCGGAGAGGGGAAATCCCCTTTTCCGGTGATTGTCAGAGAGACAATCGCTATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCGGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTcgcc                                                            <  1:443718/235‑1 (MQ=255)
                                                                                                                                                                                       aTTAACTGTCCGGTCGCCGGAGAGTGGAAATCCCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGa                                                       <  1:47120/238‑1 (MQ=255)
                                                                                                                                                                                           aCTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGCCATTTCCGCAACGGTATAAAGATTATCTTGTTGTTGTAAATCGGCAATTTGCTTTTCCCCCCCGATTaaa                                                  >  2:541012/1‑238 (MQ=255)
                                                                                                                                                                                             tGTCCGGTCGGCGGAGGGTGGAAATCCCCTTTTGGGGGGATTGGCAGAGTGACAATCGCGATTTGGGGTTCACACATCGCTGCCAACCCGGTTTATAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACTTCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGtt                                                <  2:345991/239‑1 (MQ=255)
                                                                                                                                                                                             ctttcgttCGCCGGGGATTTAAAATACCCTTTTTCGGTGTTTGCCAGAAAGACTAGCCGATTTTGCGGTTCACAACTCGCTTCCATCACGGTTTCAAAGCCAACTATCTGTACGTGCAAGGCTTTTTTAACGGCGCCAACGACGCGAGCCGTCCCCACTTCAGCCGACATATCCGCAACGGTATAAAGATTATCTTGCTTTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGtt                                                <  2:376382/232‑1 (MQ=255)
                                                                                                                                                                                              gTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGGGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGAAGCGAGCCGTCCACACATCGGCCGACCTTTCCGCCACGGTCTAAAGATTATCTTTATGTTGTAAATCGGCAATTTTCTGTTGGCCCACGGTTaattg                                                >  2:294034/1‑235 (MQ=255)
                                                                                                                                                                                                 cgggccccGGGGAGTGGAAATCCCCTTTTTCGGTTATTGTGAGGGATACAAACGGGTTTTGCGGTTCCCAAATGCGTGCCAACACGGGTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGTCGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGAcc                                            <  2:25794/233‑1 (MQ=255)
                                                                                                                                                                                                 cGGTCGCCGGGGAGTGGAAATTCCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACAGCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGAcc                                            <  2:420480/239‑1 (MQ=255)
                                                                                                                                                                                                 cGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACCCGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGAcc                                            <  1:687716/239‑1 (MQ=255)
                                                                                                                                                                                                  ggTCGCCGGAGAGTGGAAATTCCCTTTTTCGGTGATTGTCAGAGAGACAATCGAGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCt                                           <  1:169191/239‑1 (MQ=255)
                                                                                                                                                                                                  ggTCCCCGGAGAGGGGAAATACCCTTTTTCGGTGATTGTCAGGGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACAAGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGAAGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCt                                           <  1:294034/239‑1 (MQ=255)
                                                                                                                                                                                                     cGCCGGAGAGTTGAAATACCCTTTTGCGGTTATTGGCAGGAAGCCAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCat                                        <  1:901985/239‑1 (MQ=255)
                                                                                                                                                                                                              gTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGACGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTTTAAATCGGCAATTTGCTGTTCACAGCCGATTAAGTTGACCTAATAAAATcccc                               >  2:693436/1‑239 (MQ=255)
                                                                                                                                                                                                                ggAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCATAATATcccc                               <  1:668011/237‑1 (MQ=255)
                                                                                                                                                                                                                   aaTCCCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCATAATATCCCCAGTCa                          <  1:709287/239‑1 (MQ=255)
                                                                                                                                                                                                                       cccTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAAc                                                                                                            <  1:319699/153‑1 (MQ=255)
                                                                                                                                                                                                                       cccTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAAc                                                                                                            >  2:319699/1‑153 (MQ=255)
                                                                                                                                                                                                                       cccTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCATAATATCCCCAGTCACTGa                      >  1:109078/1‑239 (MQ=255)
                                                                                                                                                                                                                            tttCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCATAATATCCCCAGTCACTGCAATGt                 >  1:273867/1‑239 (MQ=255)
                                                                                                                                                                                                                                     ttttCAGAGAGTCAATCGCGATTTGCGGTTCCCACACCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAACCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCATAATATCCCCAGTCACTGAAATGTTCCGTAGc         <  2:348861/235‑1 (MQ=255)
                                                                                                                                                                                                                                           gagagaCAATCGCGATTTGCGGTTCACGCATCGCTGCCAACCCGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCATAATATCCCCAGTCACTGAAATGTTCCGTAGCAAGAata  <  1:358727/239‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  
CGTAGGAAGTGACAACGTCACGCATAACATGACCGTTTTCTGGCATGTTGTCACATGACATGACGGTAAATGCGGGAAGTCCTGCCGCTTTACGGCGAGCCAGCGCCTCAACAATCACCCCTGTTGCTGTTTTCGGCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGGTATAAAGATTATCTTGCTGTTGTAAATCGGCAATTTGCTGTTCGCCGCCGATTAAGTTGACCTCATAATATCCCCAGTCACTGAAATGTTCCGTAGCAAGAATA  >  NC_000913/1630316‑1630788

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