Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,201 A→T 7.6% intergenic (+44/+27) smrB → / ← yfcO putative endonuclease SmrB/DUF2544 domain‑containing protein YfcO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,2010AT 7.6% 485.0 / 20.2 158intergenic (+44/+27)smrB/yfcOputative endonuclease SmrB/DUF2544 domain‑containing protein YfcO
Reads supporting (aligned to +/- strand):  ref base A (69/76);  new base T (5/7);  total (75/83)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 158 total aligned reads displayed.

CGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCG  >  NC_000913/2449054‑2449349
                                                                                                                                                   |                                                                                                                                                    
cGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAAt                                                                                                                                                    >  1:1570673/1‑150 (MQ=255)
acatgcggATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAAt                                                                                                                                                    >  1:735595/8‑150 (MQ=255)
 gaaTGTGATGGACATTCATCAGGCACCAAAAGAGTATGGAGGTGATTCGGCGTTGTTGGTGATGATTGAAGTCGAAGAGTGGCTGCCCACGGAGTTGACCTGAGGAGTTGAGCGCGTTGCGCCAGCTACGGGTAATGGCAGGGGAGAATg                                                                                                                                                   >  2:1121255/3‑150 (MQ=255)
     gtgATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGcccc                                                                                                                                               >  1:222180/1‑150 (MQ=255)
     gtgATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGcccc                                                                                                                                               <  2:1388115/150‑1 (MQ=255)
       gATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCcctg                                                                                                                                             <  2:616132/150‑1 (MQ=255)
        aTGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCATTCGCCcctgc                                                                                                                                            <  2:852392/150‑1 (MQ=255)
        aTGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCATTCGCCcctgc                                                                                                                                            <  1:1651465/150‑1 (MQ=255)
         tGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCcctgcc                                                                                                                                           <  2:581077/150‑1 (MQ=255)
         tGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCcctgcc                                                                                                                                           >  2:1201171/1‑150 (MQ=255)
         tGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGCGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCcctgcc                                                                                                                                           <  1:2001649/150‑1 (MQ=255)
           gCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCcctgcctg                                                                                                                                         >  2:1993930/1‑150 (MQ=255)
           gCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCcctgcctg                                                                                                                                         >  2:2400826/1‑150 (MQ=255)
           gCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCATTCGCCcctgccat                                                                                                                                         <  1:1745997/150‑3 (MQ=255)
              tttCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCGTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                      <  2:186618/150‑1 (MQ=255)
              tttCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                      <  1:655465/150‑1 (MQ=255)
              tttCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                      <  1:2106865/150‑1 (MQ=255)
              tttCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGGGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                      <  2:1472987/150‑1 (MQ=255)
                   tcaGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                      >  1:232560/1‑145 (MQ=255)
                   tcaGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGGGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTGGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCt                                                                                                                                 <  2:537780/150‑1 (MQ=255)
                      ggCACCAAAAGAGTTTGGCGGTGTTGCGGCGTTGTTGGTGTTGATTGAAGGCGAAGCGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGTTACGGGTAATGGCAGGGGCGAATGCCCCGGCCTGGTGTATCTAAg                                                                                                                              <  1:634005/150‑1 (MQ=255)
                       gCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAgg                                                                                                                             >  2:1438806/1‑150 (MQ=255)
                        cACCAAAAGAGTATGGCGGTGATGCTGCTTTGTTGGTGTTGATTTAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAggg                                                                                                                            <  2:2250347/150‑1 (MQ=255)
                        cACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAggg                                                                                                                            >  2:2321110/1‑150 (MQ=255)
                           cAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTa                                                                                                                         >  2:344121/1‑150 (MQ=255)
                            aaaaGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTAc                                                                                                                        <  2:86382/150‑1 (MQ=255)
                             aaaGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACa                                                                                                                       <  1:202896/150‑1 (MQ=255)
                              aaGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACaa                                                                                                                      <  1:714470/150‑1 (MQ=255)
                                  gTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCa                                                                                                                  <  2:36061/150‑1 (MQ=255)
                                    aTGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCAtt                                                                                                                <  2:339077/150‑1 (MQ=255)
                                     tGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTc                                                                                                               <  2:1930256/150‑1 (MQ=255)
                                     tGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTc                                                                                                               >  2:2121567/1‑150 (MQ=255)
                                      ggCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCa                                                                                                              <  1:1475042/150‑1 (MQ=255)
                                      ggCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCa                                                                                                              <  2:2196087/150‑1 (MQ=255)
                                       gCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAg                                                                                                             <  1:2199443/150‑1 (MQ=255)
                                        cGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGt                                                                                                            >  2:2483614/1‑150 (MQ=255)
                                        cGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGt                                                                                                            >  2:833999/1‑150 (MQ=255)
                                         ggTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTa                                                                                                           >  1:1013252/1‑150 (MQ=255)
                                           tGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATg                                                                                                         <  1:1727238/150‑1 (MQ=255)
                                            taTGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGt                                                                                                        <  1:1447114/149‑1 (MQ=255)
                                              tGCGGCGTTGTTGGTGTTGATTGAAGTCTAACAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCCTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTAc                                                                                                      >  2:2364585/1‑150 (MQ=255)
                                              tGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTAc                                                                                                      >  2:760487/1‑150 (MQ=255)
                                              tGCGGCGTTGTTGGTGTTAATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTAc                                                                                                      <  1:97274/150‑1 (MQ=255)
                                                cggcgTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTg                                                                                                    <  1:279077/150‑1 (MQ=255)
                                                     ttgttgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAAcc                                                                                               <  1:720658/150‑1 (MQ=255)
                                                          ggTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCaa                                                                                          >  2:2095359/1‑150 (MQ=255)
                                                             gtTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTa                                                                                       >  2:1197601/1‑150 (MQ=255)
                                                                gattgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGaa                                                                                    >  1:1224717/1‑150 (MQ=255)
                                                                 attgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAAt                                                                                   >  1:1799658/1‑150 (MQ=255)
                                                                  ttgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCGTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATc                                                                                  <  2:1072250/150‑1 (MQ=255)
                                                                   tgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCg                                                                                 <  2:1570673/150‑1 (MQ=255)
                                                                    gaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGa                                                                                >  2:478896/1‑150 (MQ=255)
                                                                     aaGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAt                                                                               <  2:1307256/150‑1 (MQ=255)
                                                                      aGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAtt                                                                              >  2:1256293/1‑150 (MQ=255)
                                                                      aGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAtt                                                                              >  1:973538/1‑150 (MQ=255)
                                                                       gTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGTGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTg                                                                             >  1:1503772/1‑150 (MQ=255)
                                                                       gTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGAGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTg                                                                             <  2:1239404/150‑1 (MQ=255)
                                                                         cGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACTGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGc                                                                           >  2:1760454/1‑150 (MQ=255)
                                                                         cGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGc                                                                           >  2:1960202/1‑150 (MQ=255)
                                                                           aaGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTg                                                                         <  1:874337/150‑1 (MQ=255)
                                                                             gagTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTgga                                                                       <  1:1277732/150‑3 (MQ=255)
                                                                             gagTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCa                                                                       <  2:186374/150‑1 (MQ=255)
                                                                             gagTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCa                                                                       <  2:1985526/150‑1 (MQ=255)
                                                                               gTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATa                                                                     >  1:2317489/1‑150 (MQ=255)
                                                                                   cTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGACCCGGCAAGTAGAATCGATTGGCTGCATActgc                                                                 >  1:1903430/1‑150 (MQ=255)
                                                                                   cTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgc                                                                 >  1:1808654/1‑150 (MQ=255)
                                                                                   cTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgc                                                                 >  2:128901/1‑150 (MQ=255)
                                                                                      ccgccgGAGTTGCCCTCAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTgg                                                              >  2:49041/1‑150 (MQ=255)
                                                                                       cgccgGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGt                                                             >  2:1454294/1‑150 (MQ=255)
                                                                                       cgccgGAGTGGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGt                                                             <  2:2068493/150‑1 (MQ=255)
                                                                                        gccgGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTa                                                            >  1:170287/1‑150 (MQ=255)
                                                                                          cgGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCGACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATc                                                          <  2:1681667/150‑1 (MQ=255)
                                                                                                 gCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgc                                                   >  1:1956462/1‑150 (MQ=255)
                                                                                                 gCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgc                                                   <  1:1201171/150‑1 (MQ=255)
                                                                                                  cccTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                  >  2:1545009/1‑150 (MQ=255)
                                                                                                  cccTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCACTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                  >  2:158110/1‑150 (MQ=255)
                                                                                                   ccTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgcac                                                 <  2:2034986/150‑1 (MQ=255)
                                                                                                      gaggagTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGAc                                              <  2:959234/150‑1 (MQ=255)
                                                                                                      gaggagTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGAc                                              >  1:1195255/1‑150 (MQ=255)
                                                                                                       aggagTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGAcc                                             <  1:579096/150‑1 (MQ=255)
                                                                                                       aggagTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGAcc                                             <  2:1508756/150‑1 (MQ=255)
                                                                                                             tGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCg                                       <  2:117731/150‑1 (MQ=255)
                                                                                                              gAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCgt                                      >  2:35072/1‑150 (MQ=255)
                                                                                                                gcgcGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCgtgt                                    >  2:623705/1‑150 (MQ=255)
                                                                                                                 cgcgTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTa                                   >  2:729743/1‑150 (MQ=255)
                                                                                                                     ttGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCACTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGTCTGCTTACTGCTGGTATCAACACGCACGACCATATCGTGTAACTg                               <  1:1121255/150‑1 (MQ=255)
                                                                                                                      tGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGt                              >  1:1176106/1‑150 (MQ=255)
                                                                                                                             ggTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCg                       <  2:1713208/150‑1 (MQ=255)
                                                                                                                               tACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGtt                     <  1:1782503/150‑1 (MQ=255)
                                                                                                                               tACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTta                    >  1:871867/1‑150 (MQ=255)
                                                                                                                                   ggTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAg                 <  1:344121/150‑1 (MQ=255)
                                                                                                                                    caccaggCAGGGGCATTCGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAgg                >  2:1900737/6‑150 (MQ=255)
                                                                                                                                     tAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAggg               >  1:1678376/1‑150 (MQ=255)
                                                                                                                                        tGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACCCGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGgtgt            >  2:469899/1‑150 (MQ=255)
                                                                                                                                        tGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGgtgt            >  2:241640/1‑150 (MQ=255)
                                                                                                                                         ggCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTc           >  2:1666335/1‑150 (MQ=255)
                                                                                                                                         ggCAGGGGCATTCGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGCCCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTc           >  2:146384/1‑150 (MQ=255)
                                                                                                                                             ggggCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGAtt       <  1:1449998/150‑1 (MQ=255)
                                                                                                                                                 cGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGc   <  1:1571415/150‑1 (MQ=255)
                                                                                                                                                 cGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGc   <  2:2386435/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCG  >  NC_000913/2449054‑2449349

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