Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,643,110 T→A 5.1% K44N (AAA→AAT) ‡ ydfT ← Qin prophage; putative antitermination protein Q

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,643,1100TA 5.1% 548.7 / 11.5 167K44N (AAA→AAT) ‡ydfTQin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base T (84/74);  new base A (6/3);  total (90/77)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 172 total aligned reads displayed.

TAAGGAATAATCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTA  >  NC_000913/1642961‑1643257
                                                                                                                                                     |                                                                                                                                                   
tAAGGAATAATCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCACCCCAAGTGGTGCt                                                                                                                                                     >  1:701632/1‑150 (MQ=255)
tAAGGAATAATCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCt                                                                                                                                                     >  1:893077/1‑150 (MQ=255)
    gAATAATCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCg                                                                                                                                                 <  1:1618375/150‑1 (MQ=255)
            tGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCCACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGcccc                                                                                                                                         >  2:822573/1‑150 (MQ=255)
               aaCCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCAGTCGAAAGCACCact                                                                                                                                      <  1:2121169/150‑3 (MQ=255)
                aCCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATcc                                                                                                                                     >  1:798774/1‑150 (MQ=255)
                 ccGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCAGTCGAAAGCACCacttg                                                                                                                                    <  1:224347/150‑5 (MQ=255)
                  cGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATcccc                                                                                                                                   <  1:2266856/150‑1 (MQ=255)
                  cGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGGCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATcccc                                                                                                                                   <  1:1038712/150‑1 (MQ=255)
                    gATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTg                                                                                                                                 >  1:1947874/1‑150 (MQ=255)
                    gATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCAGTCGAAAGCACCacttgggt                                                                                                                                 <  2:1074598/150‑8 (MQ=255)
                     aTTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTgc                                                                                                                                >  1:912238/1‑150 (MQ=255)
                      tttCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTgcg                                                                                                                               >  2:383318/1‑150 (MQ=255)
                         cttcACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATGGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTg                                                                                                                            >  2:1623779/1‑150 (MQ=255)
                          ttcACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGc                                                                                                                           >  2:1218308/1‑150 (MQ=255)
                          ttcACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTACTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGc                                                                                                                           <  2:2465250/150‑1 (MQ=255)
                           tcACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGcc                                                                                                                          <  1:329235/150‑1 (MQ=255)
                            cACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCa                                                                                                                         >  2:76138/1‑150 (MQ=255)
                              cTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACCCCAACTTTTGCCAGAAATAAGGCCACCCCAAGTGGTGCTTTAGACTGCACCATCCCCTGCGCTGCCAtt                                                                                                                       >  2:693797/1‑150 (MQ=255)
                              cTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCAtt                                                                                                                       >  2:62207/1‑150 (MQ=255)
                              cTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCAtt                                                                                                                       >  2:1526257/1‑150 (MQ=255)
                              cTCAACTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCAtt                                                                                                                       >  1:496497/1‑150 (MQ=255)
                               tCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTa                                                                                                                      >  1:1788418/1‑150 (MQ=255)
                               tCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACATTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTa                                                                                                                      <  1:2329430/150‑1 (MQ=255)
                                cAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTa                                                                                                                      >  2:320109/1‑149 (MQ=255)
                                cAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTa                                                                                                                      <  1:320109/149‑1 (MQ=255)
                                cAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGGCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTAc                                                                                                                     <  1:2147076/150‑1 (MQ=255)
                                  aTTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACAt                                                                                                                   <  2:422894/150‑1 (MQ=255)
                                         cAATTTCGGGTTATCCATTGCCATCGCGTAATTTAGGAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCATGTGGTGCTTTCGACTGCACCATCCCCTGCGATGCCATTACATCCGTAAt                                                                                                            <  2:98434/150‑1 (MQ=255)
                                         cAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAAt                                                                                                            >  2:360385/1‑150 (MQ=255)
                                          aaTGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAAtt                                                                                                           <  2:383289/150‑1 (MQ=255)
                                           aTGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTg                                                                                                          <  2:1349403/150‑1 (MQ=255)
                                            tGTCTGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGt                                                                                                         <  1:147945/150‑1 (MQ=255)
                                              tCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTa                                                                                                       <  2:190900/150‑1 (MQ=255)
                                                         aGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTg                                                                                            >  2:1049447/1‑150 (MQ=255)
                                                          gTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGc                                                                                           <  2:893077/150‑1 (MQ=255)
                                                           tGCCATCGCGTAATTTAGCAGGCCTTCAATCGCACACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACTTCCGTAATTGTTAAACAACCGGTGcc                                                                                          <  1:440604/150‑1 (MQ=255)
                                                           tGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGcc                                                                                          <  1:1399994/150‑1 (MQ=255)
                                                               aTCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTc                                                                                      >  1:1108634/1‑150 (MQ=255)
                                                                 cgcgTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGc                                                                                    >  1:956840/1‑150 (MQ=255)
                                                                  gcgTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCCCCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCt                                                                                   >  2:831962/1‑150 (MQ=255)
                                                                  gcgTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCt                                                                                   >  1:1952327/1‑150 (MQ=255)
                                                                  gcgTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCt                                                                                   >  1:80538/1‑150 (MQ=255)
                                                                  gcgTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAAAACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCt                                                                                   <  2:2276381/150‑1 (MQ=255)
                                                                   cgTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTg                                                                                  <  1:591219/150‑1 (MQ=255)
                                                                    gTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTgg                                                                                 <  2:2149145/150‑1 (MQ=255)
                                                                      aaTTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAg                                                                               >  2:778551/1‑150 (MQ=255)
                                                                        tttAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCg                                                                             >  1:2097538/1‑150 (MQ=255)
                                                                             cAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTTTACTATTTCCTGATCTGGCATTACATGAGTAATTGTTATACAAACTGTGCTTGTCGCTGGTTCGTCATc                                                                        <  2:654854/150‑1 (MQ=255)
                                                                                 ccTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTc                                                                    >  1:214964/1‑150 (MQ=255)
                                                                                  cTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCa                                                                   >  2:1714127/1‑150 (MQ=255)
                                                                                   ttCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCAGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCaa                                                                  >  2:1882309/1‑150 (MQ=255)
                                                                                       aTCGCAAACTGAGGGTCCAGAACACCAAATTTTGCCATAAATAAGCCCAACACAAGTGCTGATTTCGACTACACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAACCAACCGGTGCCTGTCGCTGGAGCGTCCTCGCAAAAtttt                                                              >  2:1282693/1‑150 (MQ=255)
                                                                                        tCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGACATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTg                                                             >  2:1699411/1‑150 (MQ=255)
                                                                                          gCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGga                                                           >  2:2327025/1‑150 (MQ=255)
                                                                                          gCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGga                                                           >  1:293100/1‑150 (MQ=255)
                                                                                                    ggTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTgggg                                                 <  2:253554/150‑1 (MQ=255)
                                                                                                    ggTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTgggg                                                 <  2:1188511/150‑1 (MQ=255)
                                                                                                    ggTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTgggg                                                 <  1:1709282/150‑1 (MQ=255)
                                                                                                     gTCCTGAACCCCTACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTggggt                                                >  2:640084/1‑149 (MQ=255)
                                                                                                     gTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGa                                                >  1:1558386/1‑150 (MQ=255)
                                                                                                       ccTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGt                                              <  2:1947874/150‑1 (MQ=255)
                                                                                                        cTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTa                                             >  1:1974968/1‑150 (MQ=255)
                                                                                                         tGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTgggg                                                 >  2:1176541/1‑145 (MQ=255)
                                                                                                         tGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTgggg                                                 <  1:1176541/145‑1 (MQ=255)
                                                                                                             cacCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATAtt                                        >  1:1362613/1‑150 (MQ=255)
                                                                                                               ccAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATAtttt                                      >  1:1847936/1‑150 (MQ=255)
                                                                                                                  aCTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGt                                   <  2:2310766/150‑1 (MQ=255)
                                                                                                                       tGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGc                              <  1:62207/150‑1 (MQ=255)
                                                                                                                        gCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCt                             <  1:76138/150‑1 (MQ=255)
                                                                                                                        gCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCt                             >  1:879005/1‑150 (MQ=255)
                                                                                                                        gCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCt                             >  1:2458048/1‑150 (MQ=255)
                                                                                                                         ccAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATACCCTGCGATGCCATAACATCAGAAATTGTTAAACAACCGGTGCCTGTCGCTGCAGCGTCATCGCTCAATTTTGGAGATTTAGGGGAGTTATATTTTGGTAAGGCtt                            >  1:1730903/1‑150 (MQ=255)
                                                                                                                              aaTAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAgg                       <  1:1913017/150‑1 (MQ=255)
                                                                                                                              aaTAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAgg                       <  1:1817816/150‑1 (MQ=255)
                                                                                                                                  aGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCa                   <  1:1255146/150‑1 (MQ=255)
                                                                                                                                  aGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCa                   <  2:1108634/150‑1 (MQ=255)
                                                                                                                                   ggCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCAt                  >  2:1943430/1‑150 (MQ=255)
                                                                                                                                    gCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATg                 >  1:2003131/1‑150 (MQ=255)
                                                                                                                                    gCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATg                 >  2:706301/1‑150 (MQ=255)
                                                                                                                                    gCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATg                 >  2:2212521/1‑150 (MQ=255)
                                                                                                                                      cAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCt               <  1:536061/150‑1 (MQ=255)
                                                                                                                                         cccAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGt            <  2:515599/150‑1 (MQ=255)
                                                                                                                                         cccAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGt            >  1:302367/1‑150 (MQ=255)
                                                                                                                                          ggggatggtgCAGTCGAAAGCACCATCCCCTGCGATGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTAATCGCTCAATTTTCGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGtt           >  2:1395368/6‑150 (MQ=255)
                                                                                                                                           cAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTtct         <  2:2097538/150‑1 (MQ=255)
                                                                                                                                            ggatggtgCAGTCGAAAGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTACTATTTTGGTAAGGCTTCAAGGTTCATGCTCGTtct         >  2:1590702/4‑150 (MQ=255)
                                                                                                                                            ggatggtgCAGTCGAAAGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTtct         >  1:485343/4‑150 (MQ=255)
                                                                                                                                            ggatggtgCAGTCGAAAGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTtct         >  1:521800/4‑150 (MQ=255)
                                                                                                                                            ggatggtgCAGTCGAAAGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTtct         >  1:2438390/4‑150 (MQ=255)
                                                                                                                                            aaGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTtct         <  2:2017660/150‑1 (MQ=255)
                                                                                                                                              gtggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTcc       <  2:296319/150‑1 (MQ=255)
                                                                                                                                              gtggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTcc       >  2:1823491/1‑150 (MQ=255)
                                                                                                                                              gtggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTcc       <  1:822573/150‑1 (MQ=255)
                                                                                                                                              gtggtgCTTTCGACTGCACCATCCCCGGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGGGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTcc       <  2:2257708/150‑1 (MQ=255)
                                                                                                                                               tggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCa      >  1:1272589/1‑150 (MQ=255)
                                                                                                                                                ggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCAc     <  1:1489169/150‑1 (MQ=255)
                                                                                                                                                ggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCAc     <  2:30365/150‑1 (MQ=255)
                                                                                                                                                  tgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACtt   <  2:2131904/150‑1 (MQ=255)
                                                                                                                                                   gCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTa  >  2:277383/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TAAGGAATAATCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTA  >  NC_000913/1642961‑1643257

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