Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,636,484 C→T 93.0% intergenic (‑117/‑272) 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,4840CT93.0% 1499.4 / 62.2 675intergenic (‑117/‑272)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 C (23/25);  new base T (300/325);  total (325/350)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 679 total aligned reads displayed.

TCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGG  >  NC_000913/1636339‑1636630
                                                                                                                                                 |                                                                                                                                                  
     cagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaaga                                                                                                                                            <  8:313983/149‑1 (MQ=11)
           ctatttGTTCACTTCCATTGTTCATTCCACCGACAAAACCAGAGAAACGAAACGACAGAGGCCAACAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaagaagaa                                                                                                                                        <  6:363682/144‑1 (MQ=1)
                     aCTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCtt                                                                                                                            <  2:62553/148‑1 (MQ=11)
                        tCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCtt                                                                                                                            <  2:335296/146‑1 (MQ=255)
                               ttCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGaa                                                                                                                   >  8:59345/1‑148 (MQ=255)
                                         ggCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATa                                                                                                         >  1:192156/1‑148 (MQ=255)
                                                                               gCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGt                                                                  >  3:201465/1‑149 (MQ=39)
                                                                                   gctttCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCtt                                                                                                                            >  3:194353/3‑87 (MQ=11)
                                                                                   gctttCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCtt                                                                                                                            >  3:119251/3‑87 (MQ=11)
                                                                                   gctttCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCg                                                                <  3:348760/145‑1 (MQ=2)
                                                                                   gctttCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGa                                                              >  3:279718/3‑149 (MQ=2)
                                                                                   gctttCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAACAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGACGTCGCAGCCCCGTAACCTGTCgg                                                               >  2:141268/3‑148 (MQ=11)
                                                                                     tttCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAATTTATGACGAAGAAGAACGGAAACGCCTAAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGAt                                                             >  4:249545/1‑148 (MQ=11)
                                                                                     tttCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGAt                                                             >  5:282719/1‑148 (MQ=12)
                                                                                     tttCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCAGCCCCGTAACCTGTCGGATc                                                            >  6:85214/1‑149 (MQ=255)
                                                                                     tttAAGCAAATGTCGTTTCCTTTCTTTTAACAGATTATTTTAAATAAAAACATTAAATTAAGACCAACAACAACTGAAACGCCTTAAAACCAAAAAATTTCATAAATAAATAAAACCCGCGAGCACGCCGCACCGAACCCTGACGGAt                                                             >  4:172490/1‑148 (MQ=255)
                                                                                       tCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATaaa                                                                                                       >  7:64137/1‑104 (MQ=11)
                                                                                       tCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCgccgcc                                                                             <  2:387647/130‑1 (MQ=12)
                                                                                       tCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCgccgcc                                                                             >  1:387647/1‑130 (MQ=12)
                                                                                       tCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCgccgcc                                                                             >  8:375825/1‑130 (MQ=12)
                                                                                         aGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCTGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGGAACCGGTCGGATCACCg                                                        >  5:363787/1‑149 (MQ=11)
                                                                                         aGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTcgccg                                                                               >  2:168120/1‑126 (MQ=12)
                                                                                           cACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCgg                                                       >  7:390000/1‑148 (MQ=255)
                                                                                             ccTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCtt                                                                                                                            >  5:118161/1‑77 (MQ=11)
                                                                                             ccTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATATCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGaaa                                                    >  6:178958/1‑149 (MQ=255)
                                                                                             ccTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGaaa                                                    <  2:210561/149‑1 (MQ=255)
                                                                                             ccTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGaaa                                                    <  1:224143/149‑1 (MQ=255)
                                                                                             ccTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGaaa                                                    <  7:192956/149‑1 (MQ=255)
                                                                                              cTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGcc                                                                                                                              <  5:319640/74‑1 (MQ=11)
                                                                                              cTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCg                                                                                                                      <  7:96576/82‑1 (MQ=11)
                                                                                              cTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAAAATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGTAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCgccgcc                                                                             <  8:84117/123‑1 (MQ=11)
                                                                                               tGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAg                                                   >  7:133214/1‑148 (MQ=255)
                                                                                                gTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTcg                                                                                  <  4:12344/116‑1 (MQ=255)
                                                                                                gTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGCAAGGa                                                 >  7:263197/1‑149 (MQ=255)
                                                                                                 tCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTc                                                                                                            <  4:19619/89‑1 (MQ=11)
                                                                                                 tCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGa                                                 <  1:58255/148‑1 (MQ=255)
                                                                                                 tCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGAc                                                <  4:44628/149‑1 (MQ=255)
                                                                                                 cgctctCCTTTCTTTTCAGGGGGTATTTTAACTAAAAACATTAAGTTACGACGAAGACGACCGGAAATGCCTTAAACCGGAAACTTTTCATAAATCGCGAAACCCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGa                                                 <  1:291714/143‑1 (MQ=14)
                                                                                                   gTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAAtttt                                                                                                             >  5:221306/1‑86 (MQ=255)
                                                                                                   gTTTCCTTTATTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAAtttt                                                                                                             <  6:221306/86‑1 (MQ=11)
                                                                                                    tttCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAACCGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCg                                             >  2:77831/1‑149 (MQ=255)
                                                                                                       ccTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAAcc                                                                                                                       <  3:215852/72‑1 (MQ=255)
                                                                                                        ctttctttTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTaaa                                         <  5:201386/149‑1 (MQ=255)
                                                                                                        ctttctttTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTaaa                                         <  3:309673/149‑1 (MQ=255)
                                                                                                         tttctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCTCCGCCCCGTAACCTGTCGGATCACCGGAAAGGa                                                 >  1:178432/1‑140 (MQ=255)
                                                                                                         tttctttTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATa                                                                                                         >  6:272181/1‑84 (MQ=11)
                                                                                                         tttctttTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCg                                                        <  6:380572/133‑1 (MQ=255)
                                                                                                         tttctttTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAg                                                   >  8:185992/1‑138 (MQ=255)
                                                                                                         tttctttTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAg                                        >  7:302312/1‑149 (MQ=255)
                                                                                                          ttctttTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGt                                       >  1:79073/1‑149 (MQ=255)
                                                                                                          ttctttTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTAATAAATAGCGAAAACCCGCGAGGTCGCAGCCCCGTAAACTGTCTGATCACCGGAAAGGACCCGTAAAGt                                       >  5:5334/1‑149 (MQ=14)
                                                                                                           tctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGt                                       <  6:204851/148‑1 (MQ=255)
                                                                                                           tctttTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAACCCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGt                                       <  3:240519/148‑1 (MQ=255)
                                                                                                            ctttTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCgccgcc                                                                             >  2:357403/1‑109 (MQ=255)
                                                                                                               ttCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGa                                     <  1:347806/146‑1 (MQ=255)
                                                                                                               ttCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATa                                   <  1:157379/148‑1 (MQ=255)
                                                                                                                tCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATaa                                  <  1:125708/148‑1 (MQ=255)
                                                                                                                 cAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAAt                                 <  7:221807/148‑1 (MQ=255)
                                                                                                                  agagggTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGAccc                                              <  8:336781/131‑1 (MQ=255)
                                                                                                                  agagggTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGAccc                                              >  7:336781/4‑134 (MQ=255)
                                                                                                                  agagggTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGa                               >  2:258890/4‑149 (MQ=255)
                                                                                                                  agagggTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGa                               <  8:270914/146‑1 (MQ=255)
                                                                                                                  agagggTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCAGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGa                               >  8:353972/4‑149 (MQ=255)
                                                                                                                    agtgTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGAtt                             <  4:109813/146‑1 (MQ=255)
                                                                                                                     gggTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGCCCCGTAAAGTGATAATGATTa                            >  4:271102/1‑149 (MQ=255)
                                                                                                                     gggTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTa                            <  6:342324/149‑1 (MQ=255)
                                                                                                                       ggaTTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAACCCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGTAAAGGACCCGTAAAGTGATAATGATTatc                          <  2:113803/147‑1 (MQ=255)
                                                                                                                       gTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTatc                          <  6:129304/149‑1 (MQ=255)
                                                                                                                          ttttAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGAGAATGATTatcctc                       >  6:46265/1‑149 (MQ=255)
                                                                                                                            ttAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCt                      <  2:186697/148‑1 (MQ=255)
                                                                                                                              aaataaaAACATTAAGTTATGACGAAGAAGAACGGAAAGGCCTTAATCCGGAAAATTTTCATTATTAGCTAACACCCGAGAGGTCGCCGCCCCGTAACCGGTCGGATCACCGGAAAGGACGCGTAAAGTGATAATGATTATCATCTAc                    >  6:300442/1‑148 (MQ=14)
                                                                                                                              aaataaaAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCg                                             >  2:89394/1‑123 (MQ=255)
                                                                                                                              aaataaaAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACa                   >  1:98463/1‑149 (MQ=255)
                                                                                                                              aaataaaAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACATGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTAc                    <  3:70332/148‑1 (MQ=255)
                                                                                                                              aaataaaAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACa                   >  7:366207/1‑149 (MQ=255)
                                                                                                                               aagaaaaaCATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACa                   <  4:217836/145‑1 (MQ=255)
                                                                                                                                ataaaAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGaaa                                                    <  8:210617/114‑1 (MQ=255)
                                                                                                                                ataaaAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTaaa                                         >  5:186018/1‑125 (MQ=255)
                                                                                                                                ataaaAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACata                 >  8:125979/1‑149 (MQ=255)
                                                                                                                                 taaaAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACatat                >  5:344189/1‑149 (MQ=255)
                                                                                                                                 taaaAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACata                 >  5:257360/1‑148 (MQ=255)
                                                                                                                                  aaaaaCATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGAt                                                             <  4:335058/103‑1 (MQ=12)
                                                                                                                                   aaaaCATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATc               <  1:258890/148‑1 (MQ=255)
                                                                                                                                   aaaaCATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCAAAAATAGCGAAAACCCCCGAGGTCGCCGCCCCGTAACCTGGCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCCTCTACGAATca              >  6:32062/1‑149 (MQ=14)
                                                                                                                                    aaaCATTAAGTTATGACGAAGAATAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATcac             <  5:8963/149‑1 (MQ=255)
                                                                                                                                     aaCATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCcg                                                                                          >  7:106394/1‑71 (MQ=9)
                                                                                                                                     aaCATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATcaca            <  1:370018/149‑1 (MQ=255)
                                                                                                                                      aCATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCGCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACaa           <  6:152296/149‑1 (MQ=255)
                                                                                                                                       cATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAc          >  5:158817/1‑149 (MQ=255)
                                                                                                                                       cATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAc          >  7:51002/1‑149 (MQ=255)
                                                                                                                                         ttAAGTTATGACGAAGAAGAACGGAAACGCCTTATACCTGAAAATTTTCATAAATAGCTAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAcgt        <  6:248898/149‑1 (MQ=14)
                                                                                                                                         ttAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAcgt        >  3:16228/1‑149 (MQ=255)
                                                                                                                                            aGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAcgtgc      <  1:3218/148‑1 (MQ=255)
                                                                                                                                             gTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAcgtgcgt    <  8:234740/149‑1 (MQ=255)
                                                                                                                                             gTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAcgtgcgt    <  5:246699/149‑1 (MQ=255)
                                                                                                                                             gTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAcgtgcgt    <  4:210133/149‑1 (MQ=255)
                                                                                                                                             gTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAcgtgcgt    <  1:100139/149‑1 (MQ=255)
                                                                                                                                              ttATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAcgtgcgtg   >  5:381027/1‑149 (MQ=255)
                                                                                                                                              ttATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGGCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAcgtgcgtg   <  2:304515/149‑1 (MQ=255)
                                                                                                                                                acgacgAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTgg  <  4:308352/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGG  >  NC_000913/1636339‑1636630

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