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

TCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGAT  >  NC_000913/1636275‑1636571
                                                                                                                                                     |                                                                                                                                                    
tCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGt                                                                                                                                                      >  8:251817/1‑149 (MQ=40)
    aGAACGGGCATTCCCTGTTCCTGACAGTTCTGAATGGTACGTATACTCGCACCTAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTc                                                                                                                                                  <  4:44707/149‑1 (MQ=255)
              ttCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTc                                                                                                                                        >  2:376768/1‑149 (MQ=255)
                     ttacTGCCAGTTCTGAATGTTACGGATACTCGCACCTAAAATGTCAGCCAGCTGCTTTTTGTTGACATCCATGGTTCATTACACGGCCAAAAACAGAGAAAGGAAACGACTGAGGACAAAAAGCTCGCTTTCAGCACCTGTCGTTTc                                                                                                                                   <  2:84423/144‑1 (MQ=11)
                        cTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCctttc                                                                                                                              <  5:147119/149‑1 (MQ=255)
                        cTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATCCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCctttc                                                                                                                              >  6:44836/1‑149 (MQ=255)
                         tGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCctttct                                                                                                                             >  1:89721/1‑149 (MQ=255)
                           ccAGTTCTGAATGGTCCCGATACTCGCACCGAAAAAGTCAGCCAGCTGCTTTTTGTTGACTTCCAATGTTCATTCCACGGCCAAATACAGAGAAATGAAACGAGAGAGGCCAATAAGCTCGCTTTCAGCACCTGTCGTTTGctttcttt                                                                                                                           <  6:237630/149‑1 (MQ=11)
                           ccAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCTCAAAAGCTCGTTTTCAGCACCTGTCGTTTCctttcttt                                                                                                                           >  5:383508/1‑149 (MQ=255)
                           ccAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCctttcttt                                                                                                                           >  3:52298/1‑149 (MQ=255)
                            cAGTTCTGAATGGTACGGATACTCGCACCGATAATGTCTGCAAGCTGATTTTTGTTGACTTCCATTGTTCATTCCACGTCCAAAAACAGAGAACGGCAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCctttcttt                                                                                                                           >  2:275562/1‑148 (MQ=11)
                               ttCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCTCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAg                                                                                                                       <  6:383508/149‑1 (MQ=255)
                                 cttaaTGGTACTGATACTCGCACCTAAAATGTCAGCCTGCTGGTTTTTGTTGACTTCCATTGTTCATTCCACTTCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAgg                                                                                                                      <  7:153846/145‑1 (MQ=14)
                                  tGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTc                                                                                                                                   <  8:389924/134‑1 (MQ=255)
                                  tGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTc                                                                                                                                   >  7:389924/1‑134 (MQ=255)
                                     aTGGTACGGATACTCGCTCCGAAAATGTCAGCCAGCTTCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTa                                                                                                                 <  6:40814/149‑1 (MQ=11)
                                        gtccGGATACTCGCACCGAAACTGTCAGCCAGCGGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACCGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAttt                                                                                                              <  6:78401/146‑1 (MQ=11)
                                        ggaCGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAACACAGAGAAAGGAAACGACAGAGGCCAAAAAGCACGCTTTCAGCACCTGTCGGTTCCTTTCTATTCAGAGGGTAatt                                                                                                              >  6:25510/3‑149 (MQ=11)
                                        gTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGTGGGTAttt                                                                                                              >  6:204474/1‑149 (MQ=255)
                                        gTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAttt                                                                                                              >  6:216963/1‑149 (MQ=12)
                                            ggATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCGAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTtaaa                                                                                                          >  6:251632/1‑149 (MQ=255)
                                             gATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTtaaat                                                                                                         <  1:241067/149‑1 (MQ=12)
                                                 cTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATaaaaa                                                                                                    <  1:8371/149‑1 (MQ=255)
                                                    gCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAACACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATaaaaa                                                                                                    <  5:263841/147‑1 (MQ=255)
                                                    gCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAAc                                                                                                   <  7:251817/148‑1 (MQ=255)
                                                      aCCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGGCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACAt                                                                                                 >  8:369665/1‑148 (MQ=255)
                                                      aCCGAAAATGTCACCCAGCTGCTTTTTGGTGACTTCCATTGTTCATTCCACGGCCCAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACAtt                                                                                                <  2:112802/149‑1 (MQ=11)
                                                       ccGAAAATGTCAGCCATCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGTACAAAAACAGAGTATGGAATCGAAAGAGGCCAAAAAGCTCGATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGTTATTTTAAATAAAAACATTa                                                                                               <  1:380758/149‑1 (MQ=11)
                                                        cGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGTCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATaaa                                                                                              >  1:341879/1‑149 (MQ=11)
                                                        cGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATCCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTaa                                                                                              <  5:44836/149‑1 (MQ=255)
                                                         gaaacggTCAGCCAGCTGCTATTTGTGGACTTCCTTTGTTCATTCCACGGACAAACACTGAGACATGCAACGACAGAGGCCAACACGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAg                                                                                             <  4:60150/143‑1 (MQ=1)
                                                            actggcagccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTa                                                                                          <  8:169777/144‑1 (MQ=1)
                                                                tCAGCCAGCTGCTTTTTGTTGACTTCGATTGTTCATTCCACGGTCTAAATCTGAGAAAGGAAACTACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTTTTTTAAATAAAATCATTAAGTTatgac                                                                                      <  1:258641/149‑1 (MQ=11)
                                                                tCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTacgac                                                                                      >  1:364116/1‑149 (MQ=255)
                                                                   gccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCGAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaa                                                                                   <  2:130298/149‑1 (MQ=11)
                                                                   gccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTatgacg                                                                                     >  7:63144/1‑146 (MQ=11)
                                                                     cagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaaga                                                                                 <  8:313983/149‑1 (MQ=11)
                                                                      agcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACgaag                                                                                  >  2:384765/1‑147 (MQ=255)
                                                                      agcTGCTTTTTGTTGAACTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTAGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaagaa                                                                                <  8:297019/149‑1 (MQ=11)
                                                                           gtttttGGTGTCTTACATTGTTCATTCCTCGGAGAAATAGAGAGAATGGAAAATACAGAGGCCAAAAACCTCGCTTTGTGCTCCTGTGGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGATGAAGAAc                                                                            <  7:146030/147‑1 (MQ=255)
                                                                           ctatttGTTCACTTCCATTGTTCATTCCACCGACAAAACCAGAGAAACGAAACGACAGAGGCCAACAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaagaagaa                                                                             <  6:363682/144‑1 (MQ=1)
                                                                           ctTTTTGTTGACTTCCATTGTTCCTTCCACGGCCAAAAACACAGAAAGGAAACGACAGAGGCCAAAACGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGCTATTTTAAATAAAAACATTACGTTATGACgaagaagaa                                                                             <  4:215376/147‑1 (MQ=11)
                                                                           ctTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACg                                                                           <  1:384765/149‑1 (MQ=255)
                                                                           ctTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACg                                                                           <  1:78433/149‑1 (MQ=11)
                                                                             ttttGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACgg                                                                          >  2:203603/1‑148 (MQ=255)
                                                                              tttGGTGACTTCCATTCTTCATTCCTCGGAGAAAAACAGAGAAAGGAAACGTCTGTTGCCAAAAAGCTCTATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAGTTTAAATAAAATCATTAAGTTATGACGAAGAAGAACg                                                                           <  2:292256/146‑1 (MQ=11)
                                                                                tgttgACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGaaa                                                                       <  1:280104/148‑1 (MQ=255)
                                                                                tgttgACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGaaa                                                                       <  2:89721/148‑1 (MQ=255)
                                                                                     tctcccATTGTTCATTCCACTGACAAAAACAGAGAACCCCAACGACAGTGGCCAAAAAGCTCGCTATCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCtt                                                                 <  4:302118/145‑1 (MQ=1)
                                                                                     aCTTCCATTGTTCATTCCACGGACAAAAACAGTCAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAAAGCCtt                                                                 <  6:42106/149‑1 (MQ=11)
                                                                                     aCTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAATGCCt                                                                  <  3:94868/148‑1 (MQ=11)
                                                                                     aCTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCtt                                                                 <  2:62553/148‑1 (MQ=11)
                                                                                        tCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCtt                                                                 >  1:335296/1‑146 (MQ=255)
                                                                                        tCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCtt                                                                 <  2:335296/146‑1 (MQ=255)
                                                                                          cATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAAcc                                                            >  5:204851/1‑149 (MQ=255)
                                                                                          cATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAAcc                                                            >  7:91303/1‑149 (MQ=11)
                                                                                           aTTGTTCATTCCACGGTCAAAAACAGAGAAAGGAAACGACATCTGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCg                                                           <  6:31553/149‑1 (MQ=11)
                                                                                           aTTGTTCATTCCACGGCCAAAAACAGAGACAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCg                                                           <  6:251868/149‑1 (MQ=11)
                                                                                             tGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCgg                                                          >  4:66417/1‑148 (MQ=255)
                                                                                               ttCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGaa                                                        >  8:59345/1‑148 (MQ=255)
                                                                                               ttCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGGTTTCAGCACCTGGCGGTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGaaa                                                       <  6:5456/149‑1 (MQ=255)
                                                                                                     ccACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAAtttt                                                  >  1:21041/1‑148 (MQ=255)
                                                                                                         ggacAAAATCAGATAAATGAAACGACAGAGTCCAACAAGCTCGCTTTCAGAACCTGTCGTTTACTTTCTTTTCATAGGTTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCTTAAAATTTTCATaa                                             <  3:261272/146‑1 (MQ=255)
                                                                                                         ggCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATaa                                             <  1:376768/149‑1 (MQ=255)
                                                                                                         ggCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATa                                              >  1:192156/1‑148 (MQ=255)
                                                                                                           ccAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAAt                                           <  2:21041/149‑1 (MQ=255)
                                                                                                            cAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATa                                          <  4:263139/149‑1 (MQ=255)
                                                                                                                                 aCAGAGGCCAAAACGCTCGCTTTCCGCGCCTGTCGGGTCCTTTCTTTTCAGCGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATATTCAATAATAGCGAAAACCCGCGACGTcgcc                     >  8:74781/1‑149 (MQ=255)
                                                                                                                                    cttgCCCAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCgccgcc                  <  3:43886/146‑1 (MQ=11)
                                                                                                                                     gagCCCAACAAGCTCGCTATCAGCACCTGTCGTTTACGTTCTTTTGAGAGTGTATTTTAAATAAAAACAATTACTAATGACGAAGAAGAACGGAAAAGCCTTAACCGTTAAAATTTTCATAAATAGCGAAAACCCGCGAGTTcgccggc                  <  2:241642/147‑3 (MQ=2)
                                                                                                                                     agGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGccc                 <  5:53733/149‑1 (MQ=255)
                                                                                                                                          cAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa            >  8:166212/1‑149 (MQ=40)
                                                                                                                                          cAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa            <  2:16884/149‑1 (MQ=40)
                                                                                                                                          cAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa            >  2:128215/1‑149 (MQ=40)
                                                                                                                                           aaaaGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAAc           <  7:59345/149‑1 (MQ=39)
                                                                                                                                           aaaaGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAAc           <  3:66417/149‑1 (MQ=39)
                                                                                                                                           aaaaGCCCGCTTTCAGCACATGTCGTTTCCTTTCTTTTCATAGTGTATTTAAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAAAGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAAc           <  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)
                                                                                                                                                     |                                                                                                                                                    
TCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGAT  >  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: