Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,459,976 G→T 100% intergenic (+114/‑48) argW → / → intS tRNA‑Arg/CPS‑53 (KpLE1) prophage; putative prophage CPS‑53 integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,459,9760GT98.7% 779.5 / 0.5 242intergenic (+114/‑48)argW/intStRNA‑Arg/CPS‑53 (KpLE1) prophage; putative prophage CPS‑53 integrase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (99/139);  total (101/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 242 total aligned reads displayed.

TAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTCCCTGTGTCAGG  >  CP009273/2459830‑2460121
                                                                                                                                                  |                                                                                                                                                 
tAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTAc                                                                                                                                                 <  2:2551458/149‑1 (MQ=255)
tAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTAc                                                                                                                                                 <  2:1369880/149‑1 (MQ=255)
 aaTTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACg                                                                                                                                                >  1:2942690/1‑149 (MQ=255)
 aaTTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACg                                                                                                                                                <  1:3565424/149‑1 (MQ=255)
        ggTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTaaaa                                                                                                                                         >  2:3903437/1‑149 (MQ=255)
        ggTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTaaaa                                                                                                                                         >  2:3142170/1‑149 (MQ=255)
         gTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTaaaaa                                                                                                                                        <  1:3531084/149‑1 (MQ=255)
         gTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTaaaaa                                                                                                                                        <  2:1353142/149‑1 (MQ=255)
         gTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTaaaaa                                                                                                                                        <  1:784185/149‑1 (MQ=255)
          ttCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAg                                                                                                                                       >  2:1005619/1‑149 (MQ=255)
           tCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGt                                                                                                                                      >  2:779255/1‑149 (MQ=255)
              aTTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGa                                                                                                                                   >  2:1960264/1‑149 (MQ=255)
              aTTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGa                                                                                                                                   >  1:2495819/1‑149 (MQ=255)
               ttCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGAt                                                                                                                                  <  2:416051/149‑1 (MQ=255)
               ttCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGAt                                                                                                                                  <  1:1181785/149‑1 (MQ=255)
               ttCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGAt                                                                                                                                  <  2:4078665/149‑1 (MQ=255)
               ttCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGAt                                                                                                                                  >  1:376755/1‑149 (MQ=255)
               ttCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGAt                                                                                                                                  >  1:1568589/1‑149 (MQ=255)
               ttCCTGCAGGGGACACCATTTATCAGGTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGAt                                                                                                                                  <  2:146025/149‑1 (MQ=255)
                tCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATa                                                                                                                                 <  1:2098938/149‑1 (MQ=255)
                tCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATa                                                                                                                                 >  1:1518282/1‑149 (MQ=255)
                  cTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATaaa                                                                                                                               <  1:460109/149‑1 (MQ=255)
                     cAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAaata                                                                                                                            <  2:3483510/149‑1 (MQ=255)
                      aGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAaataa                                                                                                                           <  2:3848191/149‑1 (MQ=255)
                        gggACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAAtt                                                                                                                         <  1:189193/149‑1 (MQ=255)
                        gggACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAAtt                                                                                                                         >  1:3209412/1‑149 (MQ=255)
                           acacCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTa                                                                                                                      >  1:2400865/1‑149 (MQ=255)
                             acCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTAcc                                                                                                                    <  1:3131546/149‑1 (MQ=255)
                             acCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTAcc                                                                                                                    <  1:3289486/149‑1 (MQ=255)
                              ccATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTAccc                                                                                                                   >  1:1438917/1‑149 (MQ=255)
                              ccATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTAccc                                                                                                                   >  1:1439041/1‑149 (MQ=255)
                                    aTCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGa                                                                                                             >  2:3670049/1‑149 (MQ=255)
                                    aTCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGa                                                                                                             >  2:3904339/1‑149 (MQ=255)
                                     tCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGAt                                                                                                            <  2:2470552/149‑1 (MQ=255)
                                      cAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGAtt                                                                                                           <  2:2628771/149‑1 (MQ=255)
                                      cAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGAtt                                                                                                           <  2:1834708/149‑1 (MQ=255)
                                      cAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGAtt                                                                                                           <  2:562328/149‑1 (MQ=255)
                                        gTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGAtttt                                                                                                         >  1:662749/1‑149 (MQ=255)
                                         ttCGCTCCCATCCGTACCAGTCCGCACAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGAttttt                                                                                                        <  1:2198239/149‑1 (MQ=255)
                                         ttCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGAttttt                                                                                                        <  2:3219695/149‑1 (MQ=255)
                                          tCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTa                                                                                                       >  1:3429080/1‑149 (MQ=255)
                                           cGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTAc                                                                                                      <  1:18840/149‑1 (MQ=255)
                                              tCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCa                                                                                                   >  1:1249930/1‑149 (MQ=255)
                                                      gTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCAcc                                                                                           <  2:2568363/149‑1 (MQ=255)
                                                         ccAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGtt                                                                                        >  2:111586/1‑149 (MQ=255)
                                                           aGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTaa                                                                                      >  2:2822964/1‑149 (MQ=255)
                                                            gTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAg                                                                                     >  2:387193/1‑149 (MQ=255)
                                                             tCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGc                                                                                    <  1:2230535/149‑1 (MQ=255)
                                                               cGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAg                                                                                  <  2:3277802/149‑1 (MQ=255)
                                                                  aaaaTCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGAtt                                                                               <  1:1305585/149‑1 (MQ=255)
                                                                  aaaaTCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGAtt                                                                               <  2:662842/149‑1 (MQ=255)
                                                                    aaTCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGa                                                                             >  1:3922325/1‑149 (MQ=255)
                                                                     aTCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGaa                                                                            >  1:1654986/1‑149 (MQ=255)
                                                                      tCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAag                                                                           <  1:941409/149‑1 (MQ=255)
                                                                       ccccTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAagc                                                                          >  1:3193186/1‑149 (MQ=255)
                                                                          cTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAagcagc                                                                       <  1:1969250/149‑1 (MQ=255)
                                                                           tGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAgcagca                                                                      <  2:3831817/149‑1 (MQ=255)
                                                                              atatCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAg                                                                   <  2:1287222/149‑1 (MQ=255)
                                                                               tatCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGc                                                                  <  2:2211963/149‑1 (MQ=255)
                                                                                   aaGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCgaa                                                              >  2:2224361/1‑149 (MQ=255)
                                                                                    aGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATAGAAGCAGCAAAGCCgaaa                                                             >  2:159942/1‑149 (MQ=255)
                                                                                     gCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCgaaag                                                            >  2:1522959/1‑149 (MQ=255)
                                                                                      cATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCgaaaga                                                           <  2:3502153/149‑1 (MQ=255)
                                                                                        ttCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCgaaagaaa                                                         >  1:1368695/1‑149 (MQ=255)
                                                                                        ttCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCAACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCgaaagaaa                                                         >  2:274291/1‑149 (MQ=255)
                                                                                            gTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAAcc                                                     >  1:1138479/1‑149 (MQ=255)
                                                                                              aGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCAt                                                   <  1:3789257/149‑1 (MQ=255)
                                                                                                aTTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATAc                                                 >  2:1894958/1‑149 (MQ=255)
                                                                                                 tttACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATCAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATAcc                                                >  2:3125054/1‑149 (MQ=255)
                                                                                                 tttACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATAcc                                                <  1:1643586/149‑1 (MQ=255)
                                                                                                 tttACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATAcc                                                <  1:2696877/149‑1 (MQ=255)
                                                                                                 tttACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATAcc                                                <  1:278370/149‑1 (MQ=255)
                                                                                                     cAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCt                                            <  1:696819/149‑1 (MQ=255)
                                                                                                        ttCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTtct                                         <  1:1116186/149‑1 (MQ=255)
                                                                                                        ttCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTtct                                         <  1:128374/149‑1 (MQ=255)
                                                                                                         tCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTtctc                                        <  1:3364475/149‑1 (MQ=255)
                                                                                                         tCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTtctc                                        <  1:2597703/149‑1 (MQ=255)
                                                                                                         tCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTtctc                                        >  1:2751995/1‑149 (MQ=255)
                                                                                                         tCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTtctc                                        <  2:3786541/149‑1 (MQ=255)
                                                                                                         tCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTtctc                                        <  2:3785953/149‑1 (MQ=255)
                                                                                                           gTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGa                                      <  1:1253544/149‑1 (MQ=255)
                                                                                                           gTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGa                                      <  2:2335800/149‑1 (MQ=255)
                                                                                                            tCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGAt                                     >  2:3706420/1‑149 (MQ=255)
                                                                                                              aTGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATgg                                   >  2:2456887/1‑149 (MQ=255)
                                                                                                                 gTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTaa                                <  2:294159/149‑1 (MQ=255)
                                                                                                                 gTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTaa                                <  2:833806/149‑1 (MQ=255)
                                                                                                                 gTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTaa                                <  2:3951926/149‑1 (MQ=255)
                                                                                                                    cGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATgg                             <  2:1087286/149‑1 (MQ=255)
                                                                                                                     gCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGc                            <  1:1786856/149‑1 (MQ=255)
                                                                                                                           gATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTAc                      >  1:3137164/1‑149 (MQ=255)
                                                                                                                                  gACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATg               >  1:2891569/1‑149 (MQ=255)
                                                                                                                                     aGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTcc            >  2:2712642/1‑149 (MQ=255)
                                                                                                                                      gCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTccc           >  2:868029/1‑149 (MQ=255)
                                                                                                                                      gCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTccc           <  2:696137/149‑1 (MQ=255)
                                                                                                                                      gCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTccc           >  2:3758221/1‑149 (MQ=255)
                                                                                                                                      gCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTccc           <  2:2983204/149‑1 (MQ=255)
                                                                                                                                         gCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTCCCtgt        >  2:2087463/1‑149 (MQ=255)
                                                                                                                                            cTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTCCCTGTGTc     <  2:1745215/149‑1 (MQ=255)
                                                                                                                                             tCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTCCCTGTGTCa    >  2:1440551/1‑149 (MQ=255)
                                                                                                                                               cATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTCCCTGTGTCAgg  <  1:3593572/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TAATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTCCCTGTGTCAGG  >  CP009273/2459830‑2460121

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