Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,636,423 T→C 48.4% intergenic (‑56/‑333) nohQ ← / → ynfO pseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin/uncharacterized protein, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,636,4230TC48.4% 49.1 / 75.5 79intergenic (‑56/‑333)nohQ/ynfOpseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin/uncharacterized protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base T (20/19);  new base C (10/25);  total (30/49)
Fisher's exact test for biased strand distribution p-value = 5.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGAT  >  NC_000913/1636275‑1636571
                                                                                                                                                     |                                                                                                                                                    
tCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGt                                                                                                                                                      >  8:251817/1‑149 (MQ=40)
    aGAACGGGCATTCCCTGTTCCTGACAGTTCTGAATGGTACGTATACTCGCACCTAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTc                                                                                                                                                  <  4:44707/149‑1 (MQ=255)
              ttCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTc                                                                                                                                        >  2:376768/1‑149 (MQ=255)
                     ttacTGCCAGTTCTGAATGTTACGGATACTCGCACCTAAAATGTCAGCCAGCTGCTTTTTGTTGACATCCATGGTTCATTACACGGCCAAAAACAGAGAAAGGAAACGACTGAGGACAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTc                                                                                                                                   <  2:84423/144‑1 (MQ=11)
                        cTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCctttc                                                                                                                              <  5:147119/149‑1 (MQ=255)
                        cTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATCCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCctttc                                                                                                                              >  6:44836/1‑149 (MQ=255)
                         tGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCctttct                                                                                                                             >  1:89721/1‑149 (MQ=255)
                           ccAGTTCTGAATGGTCCCGATACTCGCACCGAAAAAGTCAGCCAGCTGCTTTTTGTTGACTTCCAATGTTCATTCCACGGCCAAATACAGAGAAATGAAACGAGAGAGGCCAAT‑AAGCTCGCTTTCAGCACCTGTCGTTTGctttcttt                                                                                                                           <  6:237630/149‑1 (MQ=11)
                           ccAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCTCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCctttcttt                                                                                                                           >  5:383508/1‑149 (MQ=255)
                           ccAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCctttcttt                                                                                                                           >  3:52298/1‑149 (MQ=255)
                            cAGTTCTGAATGGTACGGATACTCGCACCGATAATGTCTGCAAGCTGATTTTTGTTGACTTCCATTGTTCATTCCACGTCCAAAAACAGAGAACGGCAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCctttcttt                                                                                                                           >  2:275562/1‑148 (MQ=11)
                               ttCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCTCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAg                                                                                                                       <  6:383508/149‑1 (MQ=255)
                                 cttaaTGGTACTGATACTCGCACCTAAAATGTCAGCCTGCTGGTTTTTGTTGACTTCCATTGTTCATTCCACTTCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAgg                                                                                                                      <  7:153846/145‑1 (MQ=14)
                                  tGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTc                                                                                                                                   <  8:389924/134‑1 (MQ=255)
                                  tGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTc                                                                                                                                   >  7:389924/1‑134 (MQ=255)
                                     aTGGTACGGATACTCGCTCCGAAAATGTCAGCCAGCTTCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTa                                                                                                                 <  6:40814/149‑1 (MQ=11)
                                        gtccGGATACTCGCACCGAAACTGTCAGCCAGCGGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACCGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAttt                                                                                                              <  6:78401/146‑1 (MQ=11)
                                        ggaCGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAACACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCACGCTTTCAGCACCTGTCGGTTCCTTTCTATTCAGAGGGTAatt                                                                                                              >  6:25510/3‑149 (MQ=11)
                                        gTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGTGGGTAttt                                                                                                              >  6:204474/1‑149 (MQ=255)
                                        gTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAttt                                                                                                              >  6:216963/1‑149 (MQ=12)
                                            ggATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCGAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTtaaa                                                                                                          >  6:251632/1‑149 (MQ=255)
                                             gATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTtaaat                                                                                                         <  1:241067/149‑1 (MQ=12)
                                                 cTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGG‑CCAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATaaaaa                                                                                                    <  1:8371/149‑1 (MQ=255)
                                                    gCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAACACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATaaaaa                                                                                                    <  5:263841/147‑1 (MQ=255)
                                                    gCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAAc                                                                                                   <  7:251817/148‑1 (MQ=255)
                                                      aCCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGGCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACAt                                                                                                 >  8:369665/1‑148 (MQ=255)
                                                      aCCGAAAATGTCACCCAGCTGCTTTTTGGTGACTTCCATTGTTCATTCCACGGCCCAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACAtt                                                                                                <  2:112802/149‑1 (MQ=11)
                                                       ccGAAAATGTCAGCCATCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGTACAAAAACAGAGTATGGAATCGAAAGAGGCCAAA‑AAGCTCGATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGTTATTTTAAATAAAAACATTa                                                                                               <  1:380758/149‑1 (MQ=11)
                                                        cGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGTCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATaaa                                                                                              >  1:341879/1‑149 (MQ=11)
                                                        cGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATCCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTaa                                                                                              <  5:44836/149‑1 (MQ=255)
                                                         gaaacggTCAGCCAGCTGCTATTTGTGGACTTCCTTTGTTCATTCCACGGACAAACACTGAGACATGCAACGACAGAGGCCAAC‑ACGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAg                                                                                             <  4:60150/143‑1 (MQ=1)
                                                            actggcagccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTa                                                                                          <  8:169777/144‑1 (MQ=1)
                                                                tCAGCCAGCTGCTTTTTGTTGACTTCGATTGTTCATTCCACGGTCTAAATCTGAGAAAGGAAACTACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTTTTTTAAATAAAATCATTAAGTTatgac                                                                                      <  1:258641/149‑1 (MQ=11)
                                                                tCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTacgac                                                                                      >  1:364116/1‑149 (MQ=255)
                                                                   gccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCGAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaa                                                                                   <  2:130298/149‑1 (MQ=11)
                                                                   gccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGG‑CCAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTatgacg                                                                                     >  7:63144/1‑146 (MQ=11)
                                                                     cagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaaga                                                                                 <  8:313983/149‑1 (MQ=11)
                                                                      agcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACgaag                                                                                  >  2:384765/1‑147 (MQ=255)
                                                                      agcTGCTTTTTGTTGAACTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTAGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaagaa                                                                                <  8:297019/149‑1 (MQ=11)
                                                                           gtttttGGTGTCTTACATTGTTCATTCCTCGGAGAAATAGAGAGAATGGAAAATACAGAGGCCAAA‑AACCTCGCTTTGTGCTCCTGTGGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGATGAAGAAc                                                                            <  7:146030/147‑1 (MQ=255)
                                                                           ctatttGTTCACTTCCATTGTTCATTCCACCGACAAAACCAGAGAAACGAAACGACAGAGGCCAAC‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaagaagaa                                                                             <  6:363682/144‑1 (MQ=1)
                                                                           ctTTTTGTTGACTTCCATTGTTCCTTCCACGGCCAAAAACACAGAAAGGAAACGACAGAGGCCAAA‑ACGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGCTATTTTAAATAAAAACATTACGTTATGACgaagaagaa                                                                             <  4:215376/147‑1 (MQ=11)
                                                                           ctTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACg                                                                           <  1:384765/149‑1 (MQ=255)
                                                                           ctTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACg                                                                           <  1:78433/149‑1 (MQ=11)
                                                                             ttttGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACgg                                                                          >  2:203603/1‑148 (MQ=255)
                                                                              tttGGTGACTTCCATTCTTCATTCCTCGGAGAAAAACAGAGAAAGGAAACGTCTGTTGCCAAA‑AAGCTCTATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAGTTTAAATAAAATCATTAAGTTATGACGAAGAAGAACg                                                                           <  2:292256/146‑1 (MQ=11)
                                                                                tgttgACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGaaa                                                                       <  1:280104/148‑1 (MQ=255)
                                                                                tgttgACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGaaa                                                                       <  2:89721/148‑1 (MQ=255)
                                                                                     tctcccATTGTTCATTCCACTGACAAAAACAGAGAACCCCAACGACAGTGGCCAAA‑AAGCTCGCTATCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCtt                                                                 <  4:302118/145‑1 (MQ=1)
                                                                                     aCTTCCATTGTTCATTCCACGGACAAAAACAGTCAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAAAGCCtt                                                                 <  6:42106/149‑1 (MQ=11)
                                                                                     aCTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAATGCCt                                                                  <  3:94868/148‑1 (MQ=11)
                                                                                     aCTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGG‑CCAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCtt                                                                 <  2:62553/148‑1 (MQ=11)
                                                                                        tCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCtt                                                                 >  1:335296/1‑146 (MQ=255)
                                                                                        tCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCtt                                                                 <  2:335296/146‑1 (MQ=255)
                                                                                          cATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAAcc                                                            >  5:204851/1‑149 (MQ=255)
                                                                                          cATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAAcc                                                            >  7:91303/1‑149 (MQ=11)
                                                                                           aTTGTTCATTCCACGGTCAAAAACAGAGAAAGGAAACGACATCTGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCg                                                           <  6:31553/149‑1 (MQ=11)
                                                                                           aTTGTTCATTCCACGGCCAAAAACAGAGACAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCg                                                           <  6:251868/149‑1 (MQ=11)
                                                                                             tGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCgg                                                          >  4:66417/1‑148 (MQ=255)
                                                                                               ttCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGaa                                                        >  8:59345/1‑148 (MQ=255)
                                                                                               ttCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGGTTTCAGCACCTGGCGGTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGaaa                                                       <  6:5456/149‑1 (MQ=255)
                                                                                                     ccACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAAtttt                                                  >  1:21041/1‑148 (MQ=255)
                                                                                                         ggacAAAATCAGATAAATGAAACGACAGAGTCCAAC‑AAGCTCGCTTTCAGAACCTGTCGTTTACTTTCTTTTCATAGGTTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCTTAAAATTTTCATaa                                             <  3:261272/146‑1 (MQ=255)
                                                                                                         ggCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATaa                                             <  1:376768/149‑1 (MQ=255)
                                                                                                         ggCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATa                                              >  1:192156/1‑148 (MQ=255)
                                                                                                           ccAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAAt                                           <  2:21041/149‑1 (MQ=255)
                                                                                                            cAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATa                                          <  4:263139/149‑1 (MQ=255)
                                                                                                                                 aCAGAGGCCAAA‑ACGCTCGCTTTCCGCGCCTGTCGGGTCCTTTCTTTTCAGCGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATATTCAATAATAGCGAAAACCCGCGACGTcgcc                     >  8:74781/1‑149 (MQ=255)
                                                                                                                                    cttgCCCAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCgccgcc                  <  3:43886/146‑1 (MQ=11)
                                                                                                                                     gagCCCAACAAGCTCGCTATCAGCACCTGTCGTTTACGTTCTTTTGAGAGTGTATTTTAAATAAAAACAATTACTAATGACGAAGAAGAACGGAAAAGCCTTAACCGTTAAAATTTTCATAAATAGCGAAAACCCGCGAGTTcgccggc                  <  2:241642/147‑3 (MQ=2)
                                                                                                                                     agGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGccc                 <  5:53733/149‑1 (MQ=255)
                                                                                                                                          cAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa            >  8:166212/1‑149 (MQ=40)
                                                                                                                                          cAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa            <  2:16884/149‑1 (MQ=40)
                                                                                                                                          cAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa            >  2:128215/1‑149 (MQ=40)
                                                                                                                                           aa‑aaGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAAc           <  7:59345/149‑1 (MQ=39)
                                                                                                                                           aa‑aaGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAAc           <  3:66417/149‑1 (MQ=39)
                                                                                                                                           aa‑aaGCCCGCTTTCAGCACATGTCGTTTCCTTTCTTTTCATAGTGTATTTAAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAAAGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAAc           <  1:241001/149‑1 (MQ=11)
                                                                                                                                                gCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGt       >  3:201465/1‑149 (MQ=39)
                                                                                                                                                  tCGCTTTCAGCACCTGTCGTTTCCATTCTTTTAAGAGGGGATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAAAGCCGTAAACCGGAAAATTTAAATAAATAGCGAAACCCCGCGAGGTCGCCGCCCCGTAACCTGTCg     >  2:183564/1‑149 (MQ=11)
                                                                                                                                                    gTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGa   >  6:376556/1‑149 (MQ=38)
                                                                                                                                                     ttttCAGCACCTGTCGTTTCCTTTCTTTTCAGTGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGACTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGATGTCGCCGCCCCGTAACCTGTCGGAt  <  5:204474/149‑1 (MQ=255)
                                                                                                                                                     ttttCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGAt  >  2:72504/1‑149 (MQ=38)
                                                                                                                                                     ttttCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGAt  <  1:72504/149‑1 (MQ=38)
                                                                                                                                                     |                                                                                                                                                    
TCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGAT  >  NC_000913/1636275‑1636571

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