Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,467,970 G→T 47.6% D154Y (GAC→TAC)  yfdQ → CPS‑53 (KpLE1) prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,467,9700GT47.6% 110.2 / 263.4 216D154Y (GAC→TAC) yfdQCPS‑53 (KpLE1) prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (60/53);  new base T (53/50);  total (113/103)
Fisher's exact test for biased strand distribution p-value = 8.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 216 total aligned reads displayed.

CAGATAACACCGCCACTCTCAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTCGAAGCGAAGACCAAAGACATTATGCC  >  CP009273/2467824‑2468115
                                                                                                                                                   |                                                                                                                                                 
     aaCACCGCCACTCTCAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCa                                                                                                                                             >  2:2906738/1‑148 (MQ=255)
      acacCGCCACTCTCAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCAt                                                                                                                                            <  1:344124/148‑1 (MQ=255)
      acacCGCCACTCTCAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCAt                                                                                                                                            >  1:3508568/1‑148 (MQ=255)
      acacCGCCACTCTCAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCAt                                                                                                                                            >  1:1744747/1‑148 (MQ=255)
          cGCCACTCTCAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAg                                                                                                                                        >  2:28567/1‑148 (MQ=255)
          cGCCACTCTCAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAg                                                                                                                                        >  2:30798/1‑148 (MQ=255)
             cACTCTCAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTAGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGca                                                                                                                                     <  2:1251238/148‑1 (MQ=255)
               ctctCAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGcaac                                                                                                                                   >  1:2036642/1‑148 (MQ=255)
               ctctCAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGcaac                                                                                                                                   >  1:591644/1‑148 (MQ=255)
                tctcAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGcaaca                                                                                                                                  <  1:1021402/148‑1 (MQ=255)
                  tcAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACaaa                                                                                                                                >  1:1262693/1‑148 (MQ=255)
                   cAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACaaaa                                                                                                                               <  2:368365/148‑1 (MQ=255)
                     aaCTGAAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCGGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAg                                                                                                                              <  1:1435297/148‑1 (MQ=255)
                       cTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCgg                                                                                                                           <  1:485222/148‑1 (MQ=255)
                        tGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGc                                                                                                                          <  1:442930/148‑1 (MQ=255)
                            aaaaGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGc                                                                                                                       >  2:1569998/1‑148 (MQ=255)
                             aaaGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCg                                                                                                                      >  2:232789/1‑148 (MQ=255)
                                gACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCa                                                                                                                   >  2:1236203/1‑148 (MQ=255)
                                gACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCa                                                                                                                   >  2:1533620/1‑148 (MQ=255)
                                 aCAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCaa                                                                                                                  >  1:1016014/1‑148 (MQ=255)
                                  cAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAAt                                                                                                                 >  1:1364487/1‑148 (MQ=255)
                                    gCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATcc                                                                                                               >  2:3415349/1‑148 (MQ=255)
                                     cACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCg                                                                                                              >  1:1322302/1‑148 (MQ=255)
                                      aCCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGt                                                                                                             <  2:3116120/148‑1 (MQ=255)
                                       ccGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTa                                                                                                            >  1:2233500/1‑148 (MQ=255)
                                        cGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTaa                                                                                                           >  1:3388977/1‑148 (MQ=255)
                                         gTTCTCTGCGGTGTTGTCTGTTAACGGCGCGCGTAACTCCCATAAATCACGGGCAGATTGGATTGAAGACTGGGCCGACTACCTGGTGGGCTTTGATGCTATTGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTaaa                                                                                                          <  1:703807/148‑1 (MQ=255)
                                          ttCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTaaaa                                                                                                         >  2:2192210/1‑148 (MQ=255)
                                           tctcTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAAt                                                                                                        >  1:76889/1‑148 (MQ=255)
                                            ctctGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGTATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTATTGGTTACGCCTTTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATc                                                                                                       <  1:2599617/148‑1 (MQ=255)
                                            ctctGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGAAAAAAAAGCGGCTGCGGCAATCCGTAAAATc                                                                                                       <  1:812393/148‑1 (MQ=255)
                                             tctGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCa                                                                                                      <  2:3167279/148‑1 (MQ=255)
                                              ctGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCAc                                                                                                     >  1:1231177/1‑148 (MQ=255)
                                                gCCCTGTTGTCTTTTACCGGCGAGCGTAACTCCAATAAATCACTGGCATATTGGATTGAAGAGGGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCAATTCAGGAAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGa                                                                                                   <  1:2716508/148‑1 (MQ=255)
                                                gCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGa                                                                                                   <  2:1986047/148‑1 (MQ=255)
                                                 cccTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGAt                                                                                                  <  2:2033419/148‑1 (MQ=255)
                                                   cTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTg                                                                                                >  1:903696/1‑148 (MQ=255)
                                                      ttgtCTGTTAACGGCGAGCGTAACTCCCCGAAATCCCTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAg                                                                                             <  2:2383757/148‑1 (MQ=255)
                                                      ttgtCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAg                                                                                             >  1:2441647/1‑148 (MQ=255)
                                                      ttgtCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAg                                                                                             <  2:1381791/148‑1 (MQ=255)
                                                       tgtCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGc                                                                                            >  2:1365772/1‑148 (MQ=255)
                                                         ttttttACCGGCGAGCGTAACTCCCATAAATCACTGGAATATTGGATTGAAGAGTGGGCCGACTACATTGGGGGATTTGATGCTAATGTTGACGCAATTCAGGAAAAAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCaa                                                                                          <  1:1528683/144‑1 (MQ=255)
                                                            gTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAAcc                                                                                       <  2:628783/148‑1 (MQ=255)
                                                             ttAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCa                                                                                      >  2:2606463/1‑148 (MQ=255)
                                                             ttAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCa                                                                                      >  1:172180/1‑148 (MQ=255)
                                                             ttAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCa                                                                                      <  2:1351044/148‑1 (MQ=255)
                                                                aCGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGAc                                                                                   <  2:865300/148‑1 (MQ=255)
                                                                  ggCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCg                                                                                 <  2:251384/148‑1 (MQ=255)
                                                                  ggCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCg                                                                                 >  1:1000956/1‑148 (MQ=255)
                                                                   gCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGc                                                                                <  2:577621/148‑1 (MQ=255)
                                                                     gAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTg                                                                              <  2:379966/148‑1 (MQ=255)
                                                                      aGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGa                                                                             <  2:1846471/148‑1 (MQ=255)
                                                                       gCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGAt                                                                            >  1:2835277/1‑148 (MQ=255)
                                                                       gCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGAt                                                                            >  2:1775804/1‑148 (MQ=255)
                                                                        cGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGAtt                                                                           >  1:3331976/1‑148 (MQ=255)
                                                                           aaCTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTCCCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTg                                                                        <  1:1125078/148‑1 (MQ=255)
                                                                            aCTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGa                                                                       >  1:2814009/1‑148 (MQ=255)
                                                                             cTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGaa                                                                      >  2:1015872/1‑148 (MQ=255)
                                                                               cccAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGa                                                                    >  2:3199929/1‑148 (MQ=255)
                                                                                ccAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGAt                                                                   >  2:2652770/1‑148 (MQ=255)
                                                                                      aTCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGAc                                                             >  1:2083279/1‑148 (MQ=255)
                                                                                      aTCACGGGCAGATGTGTTTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGATAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGAc                                                             <  2:903696/148‑1 (MQ=255)
                                                                                       tCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACt                                                            >  2:1351551/1‑148 (MQ=255)
                                                                                        cACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACtt                                                           >  2:1832718/1‑148 (MQ=255)
                                                                                        cACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACtt                                                           <  1:2192210/148‑1 (MQ=255)
                                                                                         aCTGGCATAATGGATTGAAGACTGGGCCGACTACGTTGTGGGCTTTGTTGATATTGGTTACGCAATTCAGGAAAAAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTc                                                          <  1:308783/148‑1 (MQ=255)
                                                                                             gCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCg                                                      >  1:3116909/1‑148 (MQ=255)
                                                                                              cAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCgg                                                     <  2:3495987/148‑1 (MQ=255)
                                                                                                gAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCa                                                   >  2:2657902/1‑148 (MQ=255)
                                                                                                 aaTGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAGACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCaa                                                  >  2:958155/1‑148 (MQ=255)
                                                                                                 aaTGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCac                                                  >  2:2064302/1‑147 (MQ=255)
                                                                                                 aaTGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCaa                                                  <  1:1789516/148‑1 (MQ=255)
                                                                                                 aaTGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCaa                                                  <  2:591644/148‑1 (MQ=255)
                                                                                                  aTGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCaaa                                                 <  2:3341974/148‑1 (MQ=255)
                                                                                                    ggATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACg                                               <  2:3208459/148‑1 (MQ=255)
                                                                                                    ggATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACg                                               >  2:996221/1‑148 (MQ=255)
                                                                                                    ggATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACg                                               >  2:2779581/1‑148 (MQ=255)
                                                                                                        tGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTcc                                           >  1:2321096/1‑148 (MQ=255)
                                                                                                         gAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTccc                                          <  2:2096898/148‑1 (MQ=255)
                                                                                                           aGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTg                                        <  1:1775792/148‑1 (MQ=255)
                                                                                                           aGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTg                                        <  1:1775836/148‑1 (MQ=255)
                                                                                                           aGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTg                                        <  1:3131101/148‑1 (MQ=255)
                                                                                                                gggCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGa                                   <  1:168784/148‑1 (MQ=255)
                                                                                                                     gACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTg                              <  2:564422/148‑1 (MQ=255)
                                                                                                                     gACTACCTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTg                              >  1:2892236/1‑148 (MQ=255)
                                                                                                                     gACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTg                              >  1:2502508/1‑148 (MQ=255)
                                                                                                                          cgttGTGGGCTTTGAGGATATTGTTTCCGCCATTCCGCCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTcgaa                         <  1:1527308/146‑1 (MQ=255)
                                                                                                                          ccTTGTGGGCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTcgaa                         <  1:123347/148‑1 (MQ=255)
                                                                                                                            ttGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTcgaagc                       >  1:2220911/1‑148 (MQ=255)
                                                                                                                              gtgGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTcgaagcga                     <  1:1365841/148‑1 (MQ=255)
                                                                                                                                 ggCTTTGATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTCGAAGCGAAGa                  >  2:2627150/1‑148 (MQ=255)
                                                                                                                                 ggCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTCGAAGCGAAGa                  >  2:2980111/1‑148 (MQ=255)
                                                                                                                                       gATGCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTCGAAGCGAAGACCAAAg            <  2:76889/148‑1 (MQ=255)
                                                                                                                                       gATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTCGAAGCGAAGACCAAAg            <  2:172180/148‑1 (MQ=255)
                                                                                                                                          gCTAATGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTCGAAGCGAAGACCAAAGACa         <  2:2996474/148‑1 (MQ=255)
                                                                                                                                             attGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTCGAAGCGAAGACCAAAGACATTa      <  1:3204091/146‑1 (MQ=255)
                                                                                                                                             aaTGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTCGAAGCGAAGACCAAAGACATTa      <  2:3155170/148‑1 (MQ=255)
                                                                                                                                             aaTGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCAGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTCGAAGCGAAGACCAAAGACATTa      <  2:1778471/148‑1 (MQ=255)
                                                                                                                                               tGGTTACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTCGAAGCGAAGACCAAAGACATTATg    >  2:1304892/1‑148 (MQ=255)
                                                                                                                                                 gTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTCGAAGCGAAGACCAAAGACATTATGcc  <  2:2468457/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAGATAACACCGCCACTCTCAAACTGAAAAAGACAGCACCGTTCTCTGCCCTGTTGTCTGTTAACGGCGAGCGTAACTCCCAGAAATCACTGGCAGAATGGATTGAAGACTGGGCCGACTACCTTGTGGGCTTTGATGCTAATGGTGACGCCATTCAGGCAACAAAAGCGGCTGCGGCAATCCGTAAAATCACGATTGAAGCAAACCAGACCGCTGATTTTGAAGATAATGACTTCAGCGGCAAACGCTCCCTGATGGAATCTGTCGAAGCGAAGACCAAAGACATTATGCC  >  CP009273/2467824‑2468115

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