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,9760GT96.3% 647.7 / 14.8 214intergenic (+114/‑48)argW/intStRNA‑Arg/CPS‑53 (KpLE1) prophage; putative prophage CPS‑53 integrase
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base T (98/108);  total (101/113)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 214 total aligned reads displayed.

AATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTCCCTGTGTCAG  >  CP009273/2459831‑2460120
                                                                                                                                                 |                                                                                                                                                
aaTTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTAc                                                                                                                                                >  1:2406032/1‑148 (MQ=255)
 aTTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACg                                                                                                                                               <  1:1143255/148‑1 (MQ=255)
 aTTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACg                                                                                                                                               <  2:1583778/148‑1 (MQ=255)
  ttGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACgg                                                                                                                                              <  1:1336652/148‑1 (MQ=255)
  ttGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACgg                                                                                                                                              <  2:2381531/148‑1 (MQ=255)
       ggTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTaaa                                                                                                                                         >  1:3144042/1‑148 (MQ=255)
        gTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTaaaa                                                                                                                                        <  1:2216174/148‑1 (MQ=255)
        gTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTaaaa                                                                                                                                        <  2:1781878/148‑1 (MQ=255)
         ttCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTaaaaa                                                                                                                                       >  2:1983744/1‑148 (MQ=255)
            gATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTg                                                                                                                                    <  1:2050746/148‑1 (MQ=255)
              ttCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGa                                                                                                                                  >  1:555995/1‑148 (MQ=255)
                 cTGCAGGGGACACAATATATCAGTTCTCTCCCATCCGTACCACTCCGCAAAATCCCCTGAATATCAAGCATTCAGTAGATGTACAGTTCTTCCTGGTTCTCTTCAGATCTTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATaa                                                                                                                               <  2:3097581/148‑1 (MQ=255)
                  tGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATaaa                                                                                                                              <  2:2091134/148‑1 (MQ=255)
                   gCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATaaaa                                                                                                                             <  2:1699516/148‑1 (MQ=255)
                     aGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAaata                                                                                                                           <  1:3303842/148‑1 (MQ=255)
                       gggACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGGTAAaataat                                                                                                                         >  1:733980/1‑148 (MQ=255)
                       gggACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAaataat                                                                                                                         >  2:3119032/1‑148 (MQ=255)
                           cacCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTa                                                                                                                     <  2:2658830/148‑1 (MQ=255)
                              cATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTAccc                                                                                                                  >  1:1656008/1‑148 (MQ=255)
                                   aTCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTTCCCACCgg                                                                                                             >  2:3259639/1‑148 (MQ=255)
                                   aTCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCgg                                                                                                             >  2:3137100/1‑148 (MQ=255)
                                   aTCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCgg                                                                                                             <  2:3024169/148‑1 (MQ=255)
                                    tCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGa                                                                                                            <  2:3316565/148‑1 (MQ=255)
                                    tCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGa                                                                                                            >  1:288977/1‑148 (MQ=255)
                                     cAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGAt                                                                                                           <  2:533002/148‑1 (MQ=255)
                                      aGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTATATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGAtt                                                                                                          >  1:3222654/1‑148 (MQ=255)
                                           gCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCNCTTCAGAGCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTAc                                                                                                     >  1:1786005/1‑148 (MQ=255)
                                           gCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTAc                                                                                                     <  2:3250336/148‑1 (MQ=255)
                                             tCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTAccc                                                                                                   <  2:1776821/148‑1 (MQ=255)
                                                cATCCGTACCAGTCCGCAAACTCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATg                                                                                                <  1:1874437/148‑1 (MQ=255)
                                                  tCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCACACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCt                                                                                              >  2:865859/1‑148 (MQ=255)
                                                   ccGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTc                                                                                             >  1:432052/1‑148 (MQ=255)
                                                   ccGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTc                                                                                             >  1:1296340/1‑148 (MQ=255)
                                                   ccGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTc                                                                                             <  2:1226904/148‑1 (MQ=255)
                                                    cGTACCAGTCCTCAAAATCCCCTGAATATCCAGCATTCCGTAGATTTACAGTTCGTCATTGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCa                                                                                            <  1:3379756/148‑1 (MQ=255)
                                                     gTACCAGTCCGCAAAATCCCCTGAATATCAAGCAGTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCAc                                                                                           <  2:3409305/148‑1 (MQ=255)
                                                       aCCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCg                                                                                         >  2:213/1‑148 (MQ=255)
                                                       aCCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCg                                                                                         >  2:242/1‑148 (MQ=255)
                                                       aCCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCg                                                                                         >  2:3066381/1‑148 (MQ=255)
                                                       aCCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCg                                                                                         <  1:2088919/148‑1 (MQ=255)
                                                           gTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTaa                                                                                     >  1:2817337/1‑148 (MQ=255)
                                                            tCCGCAAAATCCCCTGAATATCCCGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAg                                                                                    <  1:786039/148‑1 (MQ=255)
                                                            tCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAg                                                                                    >  2:2501253/1‑148 (MQ=255)
                                                             ccGCAAAATCCCCTGAATCTCAAGCATTCCGTAGATTTACAGTTCGTCTTGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGc                                                                                   <  2:509844/148‑1 (MQ=255)
                                                              cGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCa                                                                                  <  1:1077104/148‑1 (MQ=255)
                                                               gCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAg                                                                                 >  2:1302411/1‑148 (MQ=255)
                                                                  aaaTCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGAtt                                                                              >  1:8583/1‑148 (MQ=255)
                                                                  aaaTCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTAGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGAtt                                                                              <  1:1838469/148‑1 (MQ=255)
                                                                   aaTCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTg                                                                             <  1:2172571/148‑1 (MQ=255)
                                                                    aTCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGa                                                                            <  2:2504074/148‑1 (MQ=255)
                                                                     tCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGaa                                                                           >  2:1881405/1‑148 (MQ=255)
                                                                      ccccTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGTTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAag                                                                          <  2:2796254/148‑1 (MQ=255)
                                                                      ccccTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAag                                                                          <  2:2883825/148‑1 (MQ=255)
                                                                       cccTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAagc                                                                         <  1:1003454/148‑1 (MQ=255)
                                                                         cTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAagcag                                                                       >  2:1557062/1‑148 (MQ=255)
                                                                          tGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAagcagc                                                                      >  1:888023/1‑148 (MQ=255)
                                                                          tGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAagcagc                                                                      >  2:750571/1‑148 (MQ=255)
                                                                          tGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAagcagc                                                                      >  2:174063/1‑148 (MQ=255)
                                                                           gAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAgaagca                                                                     >  2:3479566/1‑148 (MQ=255)
                                                                                   aGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCgaa                                                             <  1:1587561/148‑1 (MQ=255)
                                                                                     cATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCgaaag                                                           <  2:1380828/148‑1 (MQ=255)
                                                                                     cATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCgaaag                                                           >  2:1948976/1‑148 (MQ=255)
                                                                                      aTTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCgaaaga                                                          >  1:540625/1‑148 (MQ=255)
                                                                                      aTTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCgaaaga                                                          <  1:3098109/148‑1 (MQ=255)
                                                                                         ccGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGaaaa                                                       >  1:2160569/1‑148 (MQ=255)
                                                                                            tAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGGTAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAAcc                                                    <  2:733980/148‑1 (MQ=255)
                                                                                            tAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAAcc                                                    >  1:513276/1‑148 (MQ=255)
                                                                                            tAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAAcc                                                    <  2:268674/148‑1 (MQ=255)
                                                                                            tAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAAcc                                                    >  2:3179330/1‑148 (MQ=255)
                                                                                             aGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCa                                                   >  1:2977907/1‑148 (MQ=255)
                                                                                             aGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCa                                                   >  1:97933/1‑148 (MQ=255)
                                                                                                 ttACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATAcc                                               >  1:2480163/1‑148 (MQ=255)
                                                                                                    cAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGcc                                            <  2:1376249/148‑1 (MQ=255)
                                                                                                    cAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGcc                                            <  1:210741/148‑1 (MQ=255)
                                                                                                      gTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCtt                                          >  1:1190636/1‑148 (MQ=255)
                                                                                                      gTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCtt                                          >  1:1808893/1‑148 (MQ=255)
                                                                                                           tCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGa                                     <  1:1732484/148‑1 (MQ=255)
                                                                                                           tCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGa                                     >  2:2816593/1‑148 (MQ=255)
                                                                                                             aTGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATg                                   >  2:2349330/1‑148 (MQ=255)
                                                                                                              tGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATgg                                  <  2:1639341/148‑1 (MQ=255)
                                                                                                               ggTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGt                                 <  1:372662/148‑1 (MQ=255)
                                                                                                               ggTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGt                                 >  2:1173071/1‑148 (MQ=255)
                                                                                                                gTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTa                                <  2:3013752/148‑1 (MQ=255)
                                                                                                                gTTCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTa                                <  2:2294624/148‑1 (MQ=255)
                                                                                                                 ttCGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTaa                               >  1:3171177/1‑148 (MQ=255)
                                                                                                                   cGCTTCAGATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATg                             <  2:615582/148‑1 (MQ=255)
                                                                                                                          gATCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTa                      >  1:533864/1‑148 (MQ=255)
                                                                                                                           aTCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTAc                     >  2:2840735/1‑148 (MQ=255)
                                                                                                                           aTCGTTGACAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGAGGGTAATGGCCTGTAc                     >  2:2500944/1‑148 (MQ=255)
                                                                                                                                  aCAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATg              <  1:255029/148‑1 (MQ=255)
                                                                                                                                  aCAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATg              <  1:2501276/148‑1 (MQ=255)
                                                                                                                                   cAGCCGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGt             <  2:560239/148‑1 (MQ=255)
                                                                                                                                     gCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTcc           <  1:1948976/148‑1 (MQ=255)
                                                                                                                                      ccGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTccc          <  1:2424450/148‑1 (MQ=255)
                                                                                                                                      ccGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTccc          <  1:661718/148‑1 (MQ=255)
                                                                                                                                      ccGCACTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTccc          >  2:1674539/1‑148 (MQ=255)
                                                                                                                                      ccGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTccc          <  1:1985078/148‑1 (MQ=255)
                                                                                                                                           cTCCATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTCCCTgtgt     <  2:3421689/148‑1 (MQ=255)
                                                                                                                                              cATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTCCCTGTGTCAg  <  1:1202341/148‑1 (MQ=255)
                                                                                                                                              cATTACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTCCCTGTGTCAg  <  2:3171153/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AATTGCAGGTTCGATTCCTGCAGGGGACACCATTTATCAGTTCGCTCCCATCCGTACCAGTCCGCAAAATCCCCTGAATATCAAGCATTCCGTAGATTTACAGTTCGTCATGGTTCGCTTCAGATCGTTGACAGCCGCACTCCATGACGGGTAAAAAGTGGATAAAATAATTTTACCCACCGGATTTTTACCCATGCTCACCGTTAAGCAGATTGAAGCAGCAAAGCCGAAAGAAAAACCATACCGCCTTCTCGATGGTAATGGCCTGTACCTTTATGTCCCTGTGTCAG  >  CP009273/2459831‑2460120

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: