Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,617,026 A→G 11.0% T90A (ACA→GCA)  yhhH → putative immunity protein for putative polymorphic toxin RhsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,617,0260AG11.0% 468.9 / 41.1 155T90A (ACA→GCA) yhhHputative immunity protein for putative polymorphic toxin RhsB
Reads supporting (aligned to +/- strand):  ref base A (86/52);  new base G (10/7);  total (96/59)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 160 total aligned reads displayed.

CGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCT  >  CP009273/3616879‑3617174
                                                                                                                                                   |                                                                                                                                                    
 gcAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTaacaa                                                                                                                                                   <  2:2375945/150‑1 (MQ=255)
    aGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAAcaacaa                                                                                                                                                <  2:4579505/150‑1 (MQ=255)
     gTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAg                                                                                                                                               <  1:972977/150‑1 (MQ=255)
         gAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAGCAACAAGTTgg                                                                                                                                           >  1:3478167/1‑150 (MQ=255)
          aaaGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTggg                                                                                                                                          >  2:2109522/1‑150 (MQ=255)
           aaGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgt                                                                                                                                         >  1:4399849/1‑150 (MQ=255)
           aaGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgt                                                                                                                                         >  1:4400548/1‑150 (MQ=255)
             gTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAGCAACAAGTTGGgttg                                                                                                                                       <  1:1943205/150‑1 (MQ=255)
             gTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAGCAACAAGTTGGgttg                                                                                                                                       <  1:1943208/150‑1 (MQ=255)
              tGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgttgt                                                                                                                                      >  2:3346320/1‑150 (MQ=255)
               gAAGGAAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgttgtt                                                                                                                                     <  1:2268055/150‑1 (MQ=255)
               gAAGCTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgttgtt                                                                                                                                     <  1:4581417/150‑1 (MQ=255)
                aaGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGttgttg                                                                                                                                    >  2:3760713/1‑150 (MQ=255)
                 acctAAAAATACCGTAATAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGATATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGa                                                                                                                                   <  2:4355274/147‑1 (MQ=255)
                   gTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGCTGGGTTGTTGAAg                                                                                                                                 <  1:2951432/150‑1 (MQ=255)
                   ctAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGCTGGGTTGTTGAAg                                                                                                                                 <  1:2951197/149‑1 (MQ=255)
                   ctAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGCTGGGTTGTTGAAg                                                                                                                                 <  1:2344314/149‑1 (MQ=255)
                    tAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAgg                                                                                                                                <  2:1727939/150‑1 (MQ=255)
                     aaaTATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGt                                                                                                                               <  2:1291205/150‑1 (MQ=255)
                       atatACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTAc                                                                                                                             <  1:2639769/150‑1 (MQ=255)
                              gTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCAc                                                                                                                      >  1:1690145/1‑150 (MQ=255)
                                tAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTc                                                                                                                    >  2:4054152/1‑150 (MQ=255)
                                tAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTc                                                                                                                    <  1:958442/150‑1 (MQ=255)
                                tAGAAATTTGGTTAAAAATAAAGCTAGAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTc                                                                                                                    <  1:958422/150‑1 (MQ=255)
                                       ttGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCaa                                                                                                             >  1:919841/1‑150 (MQ=255)
                                        tGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCaaa                                                                                                            >  1:3495555/1‑150 (MQ=255)
                                               aaaTAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCtgg                                                                                                     <  1:3379269/150‑1 (MQ=255)
                                                aaTAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAGCAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCtggt                                                                                                    <  1:1508429/150‑1 (MQ=255)
                                                aaTAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCtggt                                                                                                    <  1:3893732/150‑1 (MQ=255)
                                                aaTAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCtggt                                                                                                    <  1:3500733/150‑1 (MQ=255)
                                                aaTAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCtggt                                                                                                    <  1:3894443/150‑1 (MQ=255)
                                                aaTAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCtggt                                                                                                    >  2:4445471/1‑150 (MQ=255)
                                                   aaaGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTggg                                                                                                 >  2:4692134/1‑150 (MQ=255)
                                                      gCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGtt                                                                                              <  2:1600575/150‑1 (MQ=255)
                                                      gCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGtt                                                                                              <  2:3504069/150‑1 (MQ=255)
                                                       cTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCTTTATGTACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGATAACAACAAGTTGGGTTGTTGAAGGTACCATACACTCAGACCAAATTGCTGGTGGGGttt                                                                                             >  2:2448976/1‑150 (MQ=255)
                                                       cTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGttt                                                                                             >  2:679642/1‑150 (MQ=255)
                                                       cTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGttt                                                                                             >  2:3049459/1‑150 (MQ=255)
                                                       cTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGttt                                                                                             >  2:3049452/1‑150 (MQ=255)
                                                        tataGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGtttt                                                                                            >  1:683969/1‑150 (MQ=255)
                                                         ataGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGttttt                                                                                           >  1:432187/1‑150 (MQ=255)
                                                         ataGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGttttt                                                                                           <  1:263944/150‑1 (MQ=255)
                                                          taGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAGCAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTtta                                                                                          <  1:1190271/150‑1 (MQ=255)
                                                              aTTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTtattat                                                                                      >  1:2522007/1‑150 (MQ=255)
                                                              aTTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTtattat                                                                                      >  2:3650645/1‑150 (MQ=255)
                                                                  gAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGaa                                                                                  >  2:3641431/1‑150 (MQ=255)
                                                                  gAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGaa                                                                                  >  2:3641445/1‑150 (MQ=255)
                                                                   aaTTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGaaa                                                                                 >  1:1794230/1‑150 (MQ=255)
                                                                    aTTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAAt                                                                                <  2:3121849/150‑1 (MQ=255)
                                                                       aGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAgg                                                                             >  2:1394568/1‑150 (MQ=255)
                                                                          tGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCaa                                                                          >  2:820166/1‑150 (MQ=255)
                                                                            agaTATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAGCAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCaaaa                                                                        >  1:1435037/1‑150 (MQ=255)
                                                                            agaTATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAGCAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCaaaa                                                                        >  1:1434997/1‑150 (MQ=255)
                                                                                    gTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATggg                                                                >  1:2795415/1‑150 (MQ=255)
                                                                                       aaGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGgaga                                                             >  2:3458745/1‑150 (MQ=255)
                                                                                         gAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAGCAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAAt                                                           <  2:3204531/150‑1 (MQ=255)
                                                                                          aaTCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAAtt                                                          >  1:4396403/1‑150 (MQ=255)
                                                                                          aaTCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAAtt                                                          >  1:1386314/1‑150 (MQ=255)
                                                                                          aaTCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAAtt                                                          >  2:1675252/1‑150 (MQ=255)
                                                                                               ttATGGACAGGATGCCTCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGaa                                                     <  2:4710728/150‑1 (MQ=255)
                                                                                                tATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAt                                                    >  2:1042794/1‑150 (MQ=255)
                                                                                                tATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAt                                                    <  2:532779/150‑1 (MQ=255)
                                                                                                tATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAt                                                    >  1:3286813/1‑150 (MQ=255)
                                                                                                tATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAt                                                    >  2:2038353/1‑150 (MQ=255)
                                                                                                            gCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGa                                        <  1:491948/150‑1 (MQ=255)
                                                                                                             ccGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGaa                                       <  1:1101863/150‑1 (MQ=255)
                                                                                                                cAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAAt                                    >  1:3708847/1‑150 (MQ=255)
                                                                                                                 ataagaagaaAAAACATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATa                                   <  2:4099364/148‑1 (MQ=255)
                                                                                                                 agaagaagaAAAACCATACGAAATCACTGAGTTAGCAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATa                                   <  2:201617/150‑1 (MQ=255)
                                                                                                                 agaagaagaAAAACCATACGAAATCACTGAGTTAGCAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGAGGGGAGAATTCTGAATTTTGGCCACGGACAATa                                   >  1:3123667/1‑150 (MQ=255)
                                                                                                                 agaagaagaAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATa                                   >  1:4439268/1‑150 (MQ=255)
                                                                                                                 agaagaagaAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATa                                   <  2:4099298/150‑1 (MQ=255)
                                                                                                                   aagaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCATACCAAATTGCTGGTGGGGTTTTTATTATAGAAATATGCAAAAATGATTGTAGAATTATTAATTTTTGACACTGAAAATaat                                 >  1:2043284/1‑149 (MQ=255)
                                                                                                                   aagaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAg                                 >  2:1349843/1‑150 (MQ=255)
                                                                                                                    agaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGt                                <  1:2404523/150‑1 (MQ=255)
                                                                                                                     gaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTa                               >  2:1351241/1‑150 (MQ=255)
                                                                                                                      aagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTAt                              >  1:1332540/1‑150 (MQ=255)
                                                                                                                      aagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTAt                              >  1:1332557/1‑150 (MQ=255)
                                                                                                                      aagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTAt                              <  2:432187/150‑1 (MQ=255)
                                                                                                                       agaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATg                             >  1:2760697/1‑150 (MQ=255)
                                                                                                                         aaaaaCCATACGAAATCACTGAGTTAGCAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTa                           >  1:3682732/1‑150 (MQ=255)
                                                                                                                         aaaaaCCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTa                           >  1:1938825/1‑150 (MQ=255)
                                                                                                                         aaaaaCCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTa                           <  2:3639962/150‑1 (MQ=255)
                                                                                                                            aaCCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAgg                        >  1:4431625/1‑150 (MQ=255)
                                                                                                                             aCCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGa                       <  1:3652991/150‑1 (MQ=255)
                                                                                                                               cATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATa                     >  2:4366380/1‑150 (MQ=255)
                                                                                                                                aTACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATaa                    >  1:1902174/1‑150 (MQ=255)
                                                                                                                                aTACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATaa                    <  1:3278865/150‑1 (MQ=255)
                                                                                                                                 tACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAAt                   >  2:3144445/1‑150 (MQ=255)
                                                                                                                                  aCGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATc                  >  2:3094304/1‑150 (MQ=255)
                                                                                                                                   cGAAATCACTGAGTTAGCAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCa                 >  1:3276731/1‑150 (MQ=255)
                                                                                                                                    gAAATCACTGAGTTAGCAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGACTTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCaa                >  1:4148895/1‑150 (MQ=255)
                                                                                                                                    gAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGTTAATCaa                >  2:3887213/1‑150 (MQ=255)
                                                                                                                                     aaaTCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCaaa               >  2:3509249/1‑150 (MQ=255)
                                                                                                                                     aaaTCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCaaa               >  2:1236266/1‑150 (MQ=255)
                                                                                                                                     aaaTCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCaaa               >  2:1120619/1‑150 (MQ=255)
                                                                                                                                        tCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACAt            >  2:4706508/1‑150 (MQ=255)
                                                                                                                                            tGAGTTAACAACAAGCTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCgtg        >  2:2853876/1‑150 (MQ=255)
                                                                                                                                             gAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCgtgt       <  1:2277126/150‑1 (MQ=255)
                                                                                                                                                  aacaacAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCt  >  2:1479470/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCT  >  CP009273/3616879‑3617174

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: