Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,427 C→A 26.8% Q433K (CAA→AAA)  gtrS → serotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,470,4270CA26.8% 182.8 / 23.8 95Q433K (CAA→AAA) gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base C (26/41);  new base A (26/2);  total (52/43)
Fisher's exact test for biased strand distribution p-value = 7.53e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGC  >  NC_000913/2470284‑2470575
                                                                                                                                               |                                                                                                                                                     
taTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaa                                                                                                                                                <  1:1851498/150‑1 (MQ=255)
  tctcTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                                                                                              <  2:236616/150‑1 (MQ=255)
  tctcTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                                                                                              <  2:234204/150‑1 (MQ=255)
    tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                                                                                            <  1:114633/150‑1 (MQ=255)
    tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                                                                                            <  1:1339187/150‑1 (MQ=255)
    tctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATAAAAAAATTAAAAGaaaaaa                                                                                                                                            >  1:2037786/1‑147 (MQ=255)
       ttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                                                                                                         <  1:2390994/150‑1 (MQ=255)
       ttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                                                                                                         <  1:2388523/150‑1 (MQ=255)
         ttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACa                                                                                                                                       >  2:1531439/1‑150 (MQ=255)
         ttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAATAAATAAAGaaa                                                                                                                                       >  2:1434521/1‑148 (MQ=255)
            gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTAAAAGAAATaaataaaata                                                                                                                                    >  1:2126076/1‑150 (MQ=255)
            gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAtaaataaaaaaaaaa                                                                                                                                    >  1:2126112/1‑143 (MQ=255)
            gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGACaaaa                                                                                                                                    >  2:2334487/1‑148 (MQ=255)
            gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATAAAAGAAATAAAGaaaaaa                                                                                                                                    >  1:266363/1‑148 (MQ=255)
                 tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGAt                                                                                                                               <  2:960346/150‑1 (MQ=255)
                   tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAAAAAAGACAAAAGCGATGt                                                                                                                             >  1:552166/1‑150 (MQ=255)
                   tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATCAAAGAAAAAAAGAAAAAAACAATGt                                                                                                                             >  2:1519369/1‑150 (MQ=255)
                     aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                                                                                                           <  2:959689/150‑1 (MQ=255)
                     aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                                                                                                           <  2:951833/150‑1 (MQ=255)
                     aaTTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                                                                                                           <  2:951853/150‑1 (MQ=255)
                       ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                         <  1:1038871/150‑1 (MQ=255)
                          tATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCAc                                                                                                                      <  1:34179/150‑1 (MQ=255)
                               gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATTAAATAAATAAAAGAAAAAAATAAAAAATCGATGTAAAAAACTaatt                                                                                                                 >  1:1788250/1‑148 (MQ=255)
                               gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGCGATGTAAAGCACTAAGt                                                                                                                 >  2:394934/1‑150 (MQ=255)
                               gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAAAAAGAAAAAAGGGATGTTAAGCACTAAGt                                                                                                                 >  1:2385682/1‑150 (MQ=255)
                                 cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAAATTATAAAATAAATAAAATAAATAAAGACAATAGCGATGTAAAGCACAAAGatt                                                                                                               >  2:1876800/1‑150 (MQ=255)
                                  aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAACAAAACAAGACAAAAGCGATGTAAAGCACTAAGGTTa                                                                                                              >  2:260489/1‑150 (MQ=255)
                                  aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAAAAAATAAAAAAAAAAAAAAAAAAAACGATGTAAAAAACAAAATTTa                                                                                                              >  1:1625642/1‑150 (MQ=255)
                                   tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAAACAAAAGAGATGTAAAGCACTAATTTTaa                                                                                                             >  2:1245647/1‑150 (MQ=255)
                                   tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAAAAAGACAAAAGAGCTGTAAAGCACAAAGTTTaa                                                                                                             >  1:2386225/1‑150 (MQ=255)
                                    ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGACTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAACTAAACAAAAAAAAAAACGAGTAAAGCACAAAGTTTaaa                                                                                                            >  1:1823878/1‑150 (MQ=255)
                                    ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                                                                                            >  2:1018639/1‑150 (MQ=255)
                                    ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAAAAATAAAAACAATATCGATTGAAAAAAAAAAGTTTaaa                                                                                                            >  1:1926183/1‑150 (MQ=255)
                                    ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATAAAATAAATAAAAGAAATAAAAACAAAAAAGAAGTAAATAACTAAGTTTaaa                                                                                                            >  1:704320/1‑150 (MQ=255)
                                    ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAAAAGTCGTTATTTTAAATTATAAAAAAATAAAAATAAATAAATAAAAAATCCAAGTAAAGCAATAAGTTTaaa                                                                                                            >  2:248407/1‑150 (MQ=255)
                                    ttttCGGGGAAGGATATAGAGATTAAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATATTTAAAAATAATAAAAGAAATAAAAACAAAATACATGTAAAACACAAAATTTaaa                                                                                                            >  1:1823814/1‑150 (MQ=255)
                                        cGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgc                                                                                                        <  2:882861/150‑1 (MQ=255)
                                              aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCa                                                                                                  >  2:158320/1‑150 (MQ=255)
                                                     agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAAACAAAAAAATAGCCATGTAAATCACTAAATTTAAATTTCGCCCCAAACAAgg                                                                                           >  2:1128101/1‑150 (MQ=255)
                                                       agaTTTAAGCAAGCATTTATTTGGAAGGTATTTTTCGTTCGACCTTTGCTTATACATACGAGTCGTTTTTTTAATTTATAAAATAATTCAAAGGAATCAAGACAATAGCGATGTTAAGGCCTAAGTTTAAATTGCGCGCCAATCATGgcg                                                                                         >  2:157088/1‑150 (MQ=255)
                                                            tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAAATCAAAGAAAAAAAGAAAAAATAGATGTAAAACACAAAATTTAAAATGAGCGCCAAAAATGTCGCGcac                                                                                    >  2:1430111/1‑150 (MQ=255)
                                                             aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAAATCAAATAAAACAAGACAAAAGCGATGTAAAGAAAGAAGTTTAAAATTCGCGCCAAACATGGCGCGcaca                                                                                   >  2:2295102/1‑150 (MQ=255)
                                                              agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa                                                                                  <  2:756465/150‑1 (MQ=255)
                                                              agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa                                                                                  <  2:1096369/150‑1 (MQ=255)
                                                              agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa                                                                                  <  2:2138257/150‑1 (MQ=255)
                                                                 aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAAAAAATAAAAAAAAAAAAGAAAAAAGCGATGTTAAACACTAAGTTTAAATTTCGCGCCAATAAAATCGCGCACaaccg                                                                               >  1:789758/1‑147 (MQ=255)
                                                                     atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAATAAATCAAGACAATAGAGATGTAAAGCACTAAGTTTAAATTGCGCGCCCATCATGGCGCGCACAAGCTataa                                                                           >  2:1910300/1‑150 (MQ=255)
                                                                     atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAATCAAGACAAAAGCGATGTAAAGCACAAAGTTTAAATTTCGCGCCAATCAATGCGCGCACAAGCTataa                                                                           >  2:463901/1‑150 (MQ=255)
                                                                     atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATTAAATAAATAAAAAAAAAAAAGAAAAAATAGATGTTAAGCACAAAGTTTAAATTTAGCGCCCAACAAGGCGCGCACaaacaaaaa                                                                           >  1:803637/1‑143 (MQ=255)
                                                                          tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTAAAATAAAAAAATAAAAAATCGATTTAAAACAAGAAGTTTAAATTTGGCGCCAAACAATGCGGGCACCAACGataaaaaaa                                                                      >  1:2042493/1‑145 (MQ=255)
                                                                            tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCTATAATaacaac                                                                    >  2:121511/1‑150 (MQ=255)
                                                                              gAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCt                                                                  <  2:11527/150‑1 (MQ=255)
                                                                                 tGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAt                                                               <  2:2259583/150‑1 (MQ=255)
                                                                                 ggTATTTTTCTTTCGACCTTTGATTATACATAACTGTAGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGAAAATAGCGATTTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAt                                                               <  1:1502980/149‑1 (MQ=255)
                                                                                            ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGTCAATAGCGATTTAAATCACAAAGTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCTATAATAACAACCTAAATTctcctcctc                                                    >  2:1880441/1‑150 (MQ=255)
                                                                                            ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATAAAGACAAAAGCGATGTAAAGCACGAAGTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCAATAAAACCAACCAAATTTatccgccgc                                                    >  2:1138292/1‑148 (MQ=255)
                                                                                                 ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                               <  1:1531439/150‑1 (MQ=255)
                                                                                                 ccTTTGATTATACATAACTGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTaga                                               <  1:1434521/150‑1 (MQ=255)
                                                                                                   tttGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATACAATAGCGATGTAAATCAATAAGTTTAAATTTCGCGCCAATCAAGGCGCGCACAAGCTATAATACCAACCAAATTTATCCTCCTCTTAGAGt                                             >  2:2099779/1‑150 (MQ=255)
                                                                                                      gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                          <  1:1865078/150‑1 (MQ=255)
                                                                                                      gCTTATAAATAACAGTCTTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                          <  1:1865082/150‑1 (MQ=255)
                                                                                                        ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAATAAAAAAGCGATGTAAATCACTAAGTTTAAATTTCGCGCCAAAAATGGCGCGCACAAGCTATAATAACAAACAAATCAATACACCTCTTAGAGTGACta                                        >  2:332659/1‑150 (MQ=255)
                                                                                                        ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGAGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCGTAGAGTGACta                                        >  2:1773664/1‑150 (MQ=255)
                                                                                                        ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAAAAAAGCGAAGTAAAGCAATAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATAACAACCTAATTTCTCCTCCTCTTAGAGTGACta                                        >  2:812452/1‑150 (MQ=255)
                                                                                                        ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGAGCGCCAATCATGGCGCGCACAAGATATAATACCACACAAATTTATCCACCGCTTAGAGTTACta                                        >  2:1478427/1‑150 (MQ=255)
                                                                                                           taCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat                                     <  2:342642/150‑1 (MQ=255)
                                                                                                           taCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTatat                                     <  2:164747/150‑1 (MQ=255)
                                                                                                              aTAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAAAGCGATGTAAATCACTAAGTTTAAATTTCGCGCCAATAAAGGCGCGCACAATCTATAATAACAAACTAATTTGTACTACACATAGAGTGACTATAtctc                                  >  2:1268552/1‑150 (MQ=255)
                                                                                                                aaCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCt                                <  2:1257011/150‑1 (MQ=255)
                                                                                                                aaCAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCt                                <  2:1256939/150‑1 (MQ=255)
                                                                                                                aaCAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGAAAATATAGATTTAAAGCAATAATTTTAAATTGCGCGCCAATCATTGCGCGCACAAGCTATAAAAACAACCTAAATTCTACTCCTCTTAGTGTGAATATATCTCCt                                >  2:1988366/1‑150 (MQ=255)
                                                                                                                aaCAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAAACAAGAAAATAGCGAGGTAAATCACGAATTTTTAATTTGGCGCCAAACATGGCGCGCACAAGGGATAATAACAAACAAATTTGTGCTCCACTTTGAGTGTCAatataggcg                                >  2:1244106/1‑145 (MQ=255)
                                                                                                                    gTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAACATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATa                            <  1:841134/150‑1 (MQ=255)
                                                                                                                       gTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGaa                         >  2:530345/1‑150 (MQ=255)
                                                                                                                        tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                        >  2:1777703/1‑150 (MQ=255)
                                                                                                                           ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                     <  2:2102076/150‑1 (MQ=255)
                                                                                                                           ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                     >  1:2192831/1‑150 (MQ=255)
                                                                                                                           ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                     <  1:214790/150‑1 (MQ=255)
                                                                                                                           ttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGc                     <  2:2102075/150‑1 (MQ=255)
                                                                                                                            tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                    <  2:284734/150‑1 (MQ=255)
                                                                                                                            tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                    >  2:542142/1‑150 (MQ=255)
                                                                                                                            tttAATTTATAAAATAATTCAAAGAAATCAAGACAAAAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                    >  1:1693444/1‑150 (MQ=255)
                                                                                                                                 tttATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAtt               <  2:2386225/150‑1 (MQ=255)
                                                                                                                                   tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             <  1:546370/150‑1 (MQ=255)
                                                                                                                                   tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             <  1:193545/150‑1 (MQ=255)
                                                                                                                                   tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             <  1:740801/150‑1 (MQ=255)
                                                                                                                                   tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             <  1:1753250/150‑1 (MQ=255)
                                                                                                                                   tataAAATAATACAAACAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAACATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa             <  1:2136808/150‑1 (MQ=255)
                                                                                                                                          taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa      <  1:1430111/150‑1 (MQ=255)
                                                                                                                                            atTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCatat    >  1:1510238/1‑150 (MQ=255)
                                                                                                                                            atTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAAt    >  2:447451/1‑150 (MQ=255)
                                                                                                                                            atTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAAt    >  1:518710/1‑150 (MQ=255)
                                                                                                                                              tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc  <  1:2295102/150‑1 (MQ=255)
                                                                                                                                              tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc  <  1:652330/150‑1 (MQ=255)
                                                                                                                                              tCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGc  <  2:1328724/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
TATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGC  >  NC_000913/2470284‑2470575

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