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

TGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  NC_000913/2470279‑2470576
                                                                                                                                                    |                                                                                                                                                     
tGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCa                                                                                                                                                      <  1:1682611/150‑1 (MQ=255)
tGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCa                                                                                                                                                      <  2:879061/150‑1 (MQ=255)
 gtagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTaaa                                                                                                                                                     >  2:1518744/1‑150 (MQ=255)
 gtagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTaaa                                                                                                                                                     >  2:138711/1‑150 (MQ=255)
 gtagtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCaa                                                                                                                                                     >  1:2791757/1‑150 (MQ=255)
    gtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTTATTTATAAAATAATTAAAAGa                                                                                                                                                  >  2:1150902/1‑150 (MQ=255)
    gtaTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGa                                                                                                                                                  <  2:2516381/150‑1 (MQ=255)
      aTCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                                                                                                <  2:2761952/150‑1 (MQ=255)
        ctctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                                                                                              <  2:133512/150‑1 (MQ=255)
        ctctTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTGTTTTTTATTTATAAAATAAATaaaaaaaaaa                                                                                                                                              >  1:786949/1‑148 (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:1158064/150‑1 (MQ=255)
                                                                       agcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta                                                                               <  1:1157962/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:2951871/150‑1 (MQ=255)
                                                                                                              tataCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACtat                                        <  1:2951828/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                              <  1:685766/150‑1 (MQ=255)
                                                                                                                        aGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGAt                              <  2:1504281/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:2258940/150‑1 (MQ=255)
                                                                                                                                        tataAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGa              <  1:2259016/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:1021089/150‑1 (MQ=255)
                                                                                                                                                    cAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGcc  <  1:1021164/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGTAGTATCTCTTTTTGGAATATCTCAATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCC  >  NC_000913/2470279‑2470576

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: