Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,470,4360CA25.3% 216.8 / 12.8 87Q436K (CAA→AAA) yfdICPS‑53 (KpLE1) prophage; serotype specific glucosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (14/51);  new base A (22/0);  total (36/51)
Fisher's exact test for biased strand distribution p-value = 1.67e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATT  >  NC_000913/2470287‑2470582
                                                                                                                                                     |                                                                                                                                                  
ctctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                                                                                                    <  2:133512/150‑1 (MQ=255)
  ctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                                                                                                                  <  2:1997963/150‑1 (MQ=255)
   tttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                                                                                                                 <  2:718239/150‑1 (MQ=255)
   tttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTAAAAGAAATAAAg                                                                                                                                                 >  2:443438/1‑150 (MQ=255)
    ttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                                                                                                                <  1:1919113/150‑1 (MQ=255)
    ttttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                                                                                                                <  2:2568605/150‑1 (MQ=255)
     tttGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                                                                                                               <  2:1392489/150‑1 (MQ=255)
       tGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACaa                                                                                                                                             <  1:1631592/150‑1 (MQ=255)
         gAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATa                                                                                                                                           <  1:1595036/150‑1 (MQ=255)
          aaTATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAGAAATAAAAAAAATat                                                                                                                                          >  2:1313135/1‑149 (MQ=255)
              tctcAATTTTATGTGTCATTTTTCGGGGAAGGATATATAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAAAAAATAAAAACAATATCGAt                                                                                                                                      >  2:706110/1‑150 (MQ=255)
               ctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                                                                                                     <  1:2487463/150‑1 (MQ=255)
               ctcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATAAAAAAATTAAAAGAAATCAAGACAAAAGCGATg                                                                                                                                     >  1:2187755/1‑150 (MQ=255)
                tcAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAAAAAAAAACGATGt                                                                                                                                    >  1:1099685/1‑150 (MQ=255)
                 cAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAGACAATAGCGATTTa                                                                                                                                   >  2:2626108/1‑150 (MQ=255)
                    ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                                                                                                                <  1:1794768/150‑1 (MQ=255)
                     tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                                                                                                               <  1:1708168/150‑1 (MQ=255)
                      ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAATAAAAAAATAAAAAATAGATGTAAAGCa                                                                                                                              >  2:831716/1‑150 (MQ=255)
                           tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAATACAATATCGATGTTAAGCAATaat                                                                                                                         >  2:2951828/1‑149 (MQ=255)
                           tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAATAGCGATTTAAATCACTaat                                                                                                                         >  2:2951871/1‑149 (MQ=255)
                            gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTTTAAAAAAATTAAAATAAATAAATAAAATATCGAGGTAAATCACGAAGt                                                                                                                        >  1:2912581/1‑150 (MQ=255)
                             tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTATAAAATAATTAAAAAAAAAAAAAAAAAAATCGATTTAAATAACTAAAtt                                                                                                                       >  1:1526773/1‑150 (MQ=255)
                               aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAAAACAAAAGCGATTTAAAGCACTAAGTTTa                                                                                                                     >  1:2708160/1‑150 (MQ=255)
                               aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAAAAATAAATACAATATAGATGTAAATCACTAAGTTTa                                                                                                                     >  2:2462936/1‑150 (MQ=255)
                                tttttCGGGGAATGATATAGATATTTAATCAAGCATTTATTTGGAATGTATTTTTCTTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTTAGTTTaa                                                                                                                    <  1:2231051/150‑1 (MQ=255)
                                       gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                                                                                                             <  1:2270935/150‑1 (MQ=255)
                                            ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTTTAAAAAAAATAAAATAAAAAAAAACAAAAGCGATGTAAAGCACTAAGTTTTAATTTAGCGCCaa                                                                                                        >  1:440744/1‑150 (MQ=255)
                                                ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAAAAAATACAATATAGATGTAAAGCACTAAGTTTAAATTTAGCGCCAATCAt                                                                                                    >  1:1837839/1‑150 (MQ=255)
                                                ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTAATTTTTAAAAAAATTAAAAAAAATAAATAAAAAAACGATGTTAAGACCAAATTTTATATTTAGGGCCAATCaa                                                                                                    >  1:1777824/1‑149 (MQ=255)
                                                 taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                                                                                                   <  1:2520661/150‑1 (MQ=255)
                                                 taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                                                                                                   <  1:2736504/150‑1 (MQ=255)
                                                  agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg                                                                                                  <  2:1160753/150‑1 (MQ=255)
                                                   gagaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAATAAATAAATACAATAGCGATGTAAAGCACAAATTTTAAATTTCGCGCCAAACATGgc                                                                                                 >  2:2483398/1‑150 (MQ=255)
                                                        ttAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca                                                                                            <  2:2829744/150‑1 (MQ=255)
                                                           agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa                                                                                         <  1:2419065/150‑1 (MQ=255)
                                                            gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTAAAAGAAATAAAAACAAAAACGGTGTAAAGCACTAATTTTAAATATAGCGCCAATCAATGCGCGCACAAg                                                                                        >  1:2679725/1‑150 (MQ=255)
                                                            gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAATACAATATCGATTTAAATCAATAAGTTTAAATTTCGCGCCAATAATGGCGCGCACaaa                                                                                        >  2:2021624/1‑149 (MQ=255)
                                                            gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAg                                                                                        <  2:2940154/150‑1 (MQ=255)
                                                            gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATAAAGACAAAATCGATGTAAATCACTAAGTTTAAATTTCGCGCCAAAAATGTCGCGCACAAg                                                                                        >  1:726582/1‑150 (MQ=255)
                                                             caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc                                                                                       <  2:968323/150‑1 (MQ=255)
                                                               agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                                     <  1:1157962/150‑1 (MQ=255)
                                                               agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                                     <  1:1158064/150‑1 (MQ=255)
                                                                    ttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTataata                                                                                <  1:2626108/150‑1 (MQ=255)
                                                                        ttGGAATGTATTTTTCGTTCGACCTTTGCTTTACATAACAGTCGTTTTTTTAAATTATAAAATAATTCAAAGAAATAAATAAAAAATCGATGTAAAGCAATAATTTTAAATTTCGCGCCAATCATTGCGCGCACAAACTATAATACCAAc                                                                           >  2:1575589/1‑150 (MQ=255)
                                                                         tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAAc                                                                           <  2:1223273/150‑1 (MQ=255)
                                                                              tGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAAt                                                                      <  1:2220730/150‑1 (MQ=255)
                                                                              tGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAAATAAAATAAATAAATAAAAAAGCGATGTAAAACACAAAATTTAAAATTCGCGCCCATCAAGGCGCGCACAAACAAAAAAAACAACCAAAt                                                                      >  1:474514/1‑150 (MQ=255)
                                                                                  tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAAAAAAAAAAAAAAAATCGATGTAAATCACTAATTTTAAATTTCGCGCCAATAATGGCGCGCACAAACTATAATAACAAACAAAtttat                                                                  >  2:563376/1‑148 (MQ=255)
                                                                                  tttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAAAAATTCAAATAAATAAAGACAAAATCGTTGTAAAGCACTAAGTTTAAATTTCGCGCCAATAATGGCGCGCCCCAGCAATAAAAACAACCAAAtttat                                                                  >  2:563328/1‑148 (MQ=255)
                                                                                   ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTtctc                                                                 <  1:1852164/150‑1 (MQ=255)
                                                                                    tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTctcc                                                                <  2:2791757/150‑1 (MQ=255)
                                                                                          tcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCctct                                                          <  1:138711/150‑1 (MQ=255)
                                                                                            taCCTTTGATTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTa                                                        <  1:2293350/149‑1 (MQ=255)
                                                                                             aCCTTTGCTTATACATAACAGTCGTGTTTTTTATTTTTAAAATAAATCAAAGTAATAAAGACAAAAGCGATGTTAATCACGAAGTTTTAATTGAGCGCCAATCATGGCGCGCCCAAGCAATAATAACAACCAAATTTATCCACCGCTTat                                                       >  1:1332234/1‑149 (MQ=255)
                                                                                              ccTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAAAATCGATGTAAATCACGAAGTTTAAATTGAGCGCCCAAAATGGCGCGCACAAGCGATAATAACAAACTAATTTATCCTCCCCGTaga                                                      >  2:401785/1‑150 (MQ=255)
                                                                                                  tGCTTATACATAACAGTCGTTTTTTTAATTTTTAAAAAAATTCAAAGAAAACAAGAAAATATCGATGTTAATCACTAATTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCAAATTTCTGCTCCACATAGAGTGa                                                  >  2:2829255/1‑150 (MQ=255)
                                                                                                  tGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGa                                                  <  1:2483398/150‑1 (MQ=255)
                                                                                                   gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGAc                                                 <  2:2912581/150‑1 (MQ=255)
                                                                                                     ttATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACta                                               <  2:133287/150‑1 (MQ=255)
                                                                                                      tataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACtat                                              <  1:2951828/150‑1 (MQ=255)
                                                                                                      tataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACtat                                              <  1:2951871/150‑1 (MQ=255)
                                                                                                      tataAATAACAGTCTATTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCCACCTAATTTCTCCTCCTCTTAGAGTGACtat                                              <  2:2749495/150‑1 (MQ=255)
                                                                                                               cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATATCGATGTAAAGCTCTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                                     <  2:1526596/150‑1 (MQ=255)
                                                                                                               cAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGa                                     <  2:1526773/150‑1 (MQ=255)
                                                                                                                aGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAATACAATATCGATGTAAATCACTAATTTTAAATTTCGCGCCAAACATGGCGCGCACAAACGATAATAACAACCAAATTTATACTACTAATATAGTGACAATATATCCTGAt                                    >  1:2881454/1‑150 (MQ=255)
                                                                                                                aGTCGTTTTTTTAATTTATAAAATAATTCAAATAAATCAATACAATATCGATGTAAATCACTAAGTTTAAATTTCGCGCCAATCATGGCGCGCACAAGCTATAATAACAAACAAATTTATCCTCCTCTTAGATTTACTATATCTACTGAt                                    >  1:2881667/1‑150 (MQ=255)
                                                                                                                aGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAt                                    <  2:1504281/150‑1 (MQ=255)
                                                                                                                aGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAt                                    <  1:685766/150‑1 (MQ=255)
                                                                                                                  tCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAg                                  >  1:361762/1‑150 (MQ=255)
                                                                                                                     tttttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAAt                               >  1:381941/1‑150 (MQ=255)
                                                                                                                       tttttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCTCTACAATCATGTCGCGCACAAGCAATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTACTGATAGAATTg                             >  2:743705/1‑150 (MQ=255)
                                                                                                                         tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                           <  2:1335321/150‑1 (MQ=255)
                                                                                                                         tttAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCg                           <  1:1313135/150‑1 (MQ=255)
                                                                                                                                tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                    <  1:805902/150‑1 (MQ=255)
                                                                                                                                tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                    <  1:2259016/150‑1 (MQ=255)
                                                                                                                                tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa                    <  1:2258940/150‑1 (MQ=255)
                                                                                                                                 ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                   <  2:474514/150‑1 (MQ=255)
                                                                                                                                 ataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGAc                   >  1:22705/1‑150 (MQ=255)
                                                                                                                                    aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTCTTGACTAt                >  1:2657570/1‑150 (MQ=255)
                                                                                                                                     aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc               <  2:94108/150‑1 (MQ=255)
                                                                                                                                        aatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaaa            >  2:2959377/1‑150 (MQ=255)
                                                                                                                                         atTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAAt           <  1:285638/150‑1 (MQ=255)
                                                                                                                                          ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg          >  1:1364799/1‑150 (MQ=255)
                                                                                                                                            cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc        <  1:1021164/150‑1 (MQ=255)
                                                                                                                                            cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc        <  1:1021089/150‑1 (MQ=255)
                                                                                                                                                  aaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAtt  <  1:2462936/150‑1 (MQ=255)
                                                                                                                                                  aaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAtt  <  1:2845119/150‑1 (MQ=255)
                                                                                                                                                  aaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAtt  <  1:322630/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATT  >  NC_000913/2470287‑2470582

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