Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,643,111 T→G 5.4% K44T (AAA→ACA) ‡ 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,1110TG 5.4% 508.6 / 15.1 167K44T (AAA→ACA) ‡ydfTQin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base T (84/73);  new base G (6/3);  total (90/77)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 171 total aligned reads displayed.

AGGAATAATCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGC  >  NC_000913/1642963‑1643260
                                                                                                                                                    |                                                                                                                                                     
atgtttaatCTGTAAGCGGATTTCTTCACTCAATTTTTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGAAAACTGAGGGTCCTTAACACCAACTTTTGGCAGAAATAAGGCCAAGCCAAGTTGTGCttt                                                                                                                                                      <  2:1278538/145‑1 (MQ=255)
  gAATAATCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCg                                                                                                                                                    >  1:515599/1‑150 (MQ=255)
  gAATAATCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCg                                                                                                                                                    <  2:404771/150‑1 (MQ=255)
  gAATAATCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCg                                                                                                                                                    <  1:1618375/150‑1 (MQ=255)
    ataatCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGAc                                                                                                                                                  >  2:1489169/1‑150 (MQ=255)
         cTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCAc                                                                                                                                             <  1:674078/150‑1 (MQ=255)
          tGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCCACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGcccc                                                                                                                                            >  2:822573/1‑150 (MQ=255)
             aaCCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCAGTCGAAAGCACCact                                                                                                                                         <  1:2121169/150‑3 (MQ=255)
                cGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATcccc                                                                                                                                      <  1:2266856/150‑1 (MQ=255)
                cGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGGCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATcccc                                                                                                                                      <  1:1038712/150‑1 (MQ=255)
                 ggATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCt                                                                                                                                     >  1:30365/1‑150 (MQ=255)
                  gATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTg                                                                                                                                    >  1:1947874/1‑150 (MQ=255)
                  gATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCAGTCGAAAGCACCacttgggt                                                                                                                                    <  2:1074598/150‑8 (MQ=255)
                   aTTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTgc                                                                                                                                   >  1:912238/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)
                          cACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCa                                                                                                                            <  2:1560202/150‑1 (MQ=255)
                            cTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACCCCAACTTTTGCCAGAAATAAGGCCACCCCAAGTGGTGCTTTAGACTGCACCATCCCCTGCGCTGCCAtt                                                                                                                          >  2:693797/1‑150 (MQ=255)
                            cTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCAtt                                                                                                                          >  2:62207/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)
                              cAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGGCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTAc                                                                                                                        <  1:2147076/150‑1 (MQ=255)
                                aTTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACAt                                                                                                                      <  2:422894/150‑1 (MQ=255)
                                 tttGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATc                                                                                                                     >  2:13491/1‑150 (MQ=255)
                                 tttGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATc                                                                                                                     <  2:1904886/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                                                                                                             <  1:1803251/150‑1 (MQ=255)
                                          tGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGt                                                                                                            >  2:260037/1‑150 (MQ=255)
                                            tCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTa                                                                                                          <  2:190900/150‑1 (MQ=255)
                                                        gTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGc                                                                                              <  2:893077/150‑1 (MQ=255)
                                                         tGCCATCGCGTAATTTAGCAGGCCTTCAATCGCACACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACTTCCGTAATTGTTAAACAACCGGTGcc                                                                                             <  1:440604/150‑1 (MQ=255)
                                                         tGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGcc                                                                                             <  1:1399994/150‑1 (MQ=255)
                                                           ccATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTg                                                                                           <  2:1095416/150‑1 (MQ=255)
                                                             aTCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTc                                                                                         >  1:1108634/1‑150 (MQ=255)
                                                                gcgTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCCCCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCt                                                                                      >  2:831962/1‑150 (MQ=255)
                                                                gcgTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCt                                                                                      >  1:80538/1‑150 (MQ=255)
                                                                gcgTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCt                                                                                      >  1:1952327/1‑150 (MQ=255)
                                                                gcgTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAAAACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCt                                                                                      <  2:2276381/150‑1 (MQ=255)
                                                                   tAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGa                                                                                   >  1:2017660/1‑150 (MQ=255)
                                                                    aaTTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAg                                                                                  >  2:778551/1‑150 (MQ=255)
                                                                        tAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTc                                                                              >  1:2374423/1‑150 (MQ=255)
                                                                         aGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCa                                                                             >  1:1458171/1‑150 (MQ=255)
                                                                           cAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTTTACTATTTCCTGATCTGGCATTACATGAGTAATTGTTATACAAACTGTGCTTGTCGCTGGTTCGTCATc                                                                           <  2:654854/150‑1 (MQ=255)
                                                                             ggCCTTCAATCGCAAACTGAGGGTCATGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGc                                                                         >  1:2150619/1‑150 (MQ=255)
                                                                               ccTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTc                                                                       >  1:214964/1‑150 (MQ=255)
                                                                                cTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCa                                                                      >  2:1714127/1‑150 (MQ=255)
                                                                                    aaTCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAAttt                                                                  >  2:1142035/1‑150 (MQ=255)
                                                                                      tCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGACATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTg                                                                >  2:1699411/1‑150 (MQ=255)
                                                                                        gCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGga                                                              >  1:293100/1‑150 (MQ=255)
                                                                                        gCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGga                                                              >  2:2327025/1‑150 (MQ=255)
                                                                                              tGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGAtttt                                                        >  1:802682/1‑150 (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)
                                                                                                      cTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTa                                                >  1:1974968/1‑150 (MQ=255)
                                                                                                       tGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTgggg                                                    <  1:1176541/145‑1 (MQ=255)
                                                                                                           cacCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATAtt                                           >  1:574310/1‑150 (MQ=255)
                                                                                                           cacCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATAtt                                           >  1:1362613/1‑150 (MQ=255)
                                                                                                             ccAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATAtttt                                         >  1:1847936/1‑150 (MQ=255)
                                                                                                                aCTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGt                                      <  2:2310766/150‑1 (MQ=255)
                                                                                                                aCTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCCCTCAATTTTGGAGATTTTGGGGATTAATATTTTGGt                                      >  2:1041029/1‑150 (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)
                                                                                                                                aGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCa                      <  1:1255146/150‑1 (MQ=255)
                                                                                                                                aGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCa                      >  2:754511/1‑150 (MQ=255)
                                                                                                                                 ggCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCAt                     >  2:1943430/1‑150 (MQ=255)
                                                                                                                                  gCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATg                    >  2:2212521/1‑150 (MQ=255)
                                                                                                                                  gCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATg                    >  1:2003131/1‑150 (MQ=255)
                                                                                                                                   ccAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGc                   <  1:952974/150‑1 (MQ=255)
                                                                                                                                       cccAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGt               >  1:302367/1‑150 (MQ=255)
                                                                                                                                       cccAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGt               <  2:515599/150‑1 (MQ=255)
                                                                                                                                        ggggatggtgCAGTCGAAAGCACCATCCCCTGCGATGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTAATCGCTCAATTTTCGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGtt              >  2:1395368/6‑150 (MQ=255)
                                                                                                                                         cAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTtct            <  2:2097538/150‑1 (MQ=255)
                                                                                                                                          aaGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTtct            <  2:2017660/150‑1 (MQ=255)
                                                                                                                                            gtggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTcc          >  2:1823491/1‑150 (MQ=255)
                                                                                                                                            gtggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTcc          <  2:296319/150‑1 (MQ=255)
                                                                                                                                            gtggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTcc          <  1:822573/150‑1 (MQ=255)
                                                                                                                                            gtggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTcc          >  1:531376/1‑150 (MQ=255)
                                                                                                                                            gtggtgCTTTCGACTGCACCATCCCCGGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGGGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTcc          <  2:2257708/150‑1 (MQ=255)
                                                                                                                                             tggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCa         >  1:1272589/1‑150 (MQ=255)
                                                                                                                                              ggtgCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGTAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCAc        <  2:214964/150‑1 (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:1976164/1‑150 (MQ=255)
                                                                                                                                                 gCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTa     >  2:277383/1‑150 (MQ=255)
                                                                                                                                                    ttCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGc  >  1:46889/1‑150 (MQ=255)
                                                                                                                                                    ttCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGc  <  1:344456/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGGAATAATCTGTAACCGGATTTCTTCACTCAATTTGTTCAATGTCGGGTTATCCAGTGCCATCGCGTAATTTAGCAGGCCTTCAATCGCAAACTGAGGGTCCTGAACACCAACTTTTGCCAGAAATAAGGCCAACCCAAGTGGTGCTTTCGACTGCACCATCCCCTGCGCTGCCATTACATCCGTAATTGTTAAACAACCGGTGCCTGTCGCTGGAGCGTCATCGCTCAATTTTGGAGATTTTGGGGAGTAATATTTTGGTAAGGCTTCAAGGTTCATGCTCGTTCTCCACTTACGC  >  NC_000913/1642963‑1643260

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: