Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,904,553 T→A 100% P111P (CCA→CCT pstS ← periplasmic phosphate binding protein, high‑affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,904,5530TA98.8% 832.2 / ‑0.5 249P111P (CCA→CCTpstSperiplasmic phosphate binding protein, high‑affinity
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base A (130/116);  total (133/116)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 249 total aligned reads displayed.

GCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAG  >  CP009273/3904406‑3904698
                                                                                                                                                   |                                                                                                                                                 
  ggcgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGaa                                                                                                                                               >  1:1328178/1‑150 (MQ=255)
  ggcgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGaa                                                                                                                                               >  2:3728040/1‑150 (MQ=255)
   gcgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAAt                                                                                                                                              <  2:3839233/150‑1 (MQ=255)
    cgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATg                                                                                                                                             <  2:2022942/150‑1 (MQ=255)
    cgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATg                                                                                                                                             <  1:4972501/150‑1 (MQ=255)
    cgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATg                                                                                                                                             <  2:181145/150‑1 (MQ=255)
       actacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTa                                                                                                                                          >  2:629070/1‑150 (MQ=255)
       actacAGAAATGTTTTTTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGACGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTa                                                                                                                                          >  2:645663/1‑150 (MQ=255)
        ctacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTaa                                                                                                                                         >  1:2153046/1‑150 (MQ=255)
         tacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGAAGATGTCGCCGAGGGTATTACCATCCAGCACAAGATCTCCAGACTTCAGCCCAGGAATGTTAAc                                                                                                                                        >  1:4410596/1‑150 (MQ=255)
           cAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCg                                                                                                                                      >  2:4973406/1‑150 (MQ=255)
            aGCAATGTTTTGTGAAGGCAGGTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTAGATTTTGCCCAGGTATATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACATCAGCCCAGGAAAGTTAACCGc                                                                                                                                     >  2:1738240/1‑150 (MQ=255)
               aaTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAg                                                                                                                                  <  2:2927553/150‑1 (MQ=255)
               aaTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAg                                                                                                                                  >  2:764408/1‑150 (MQ=255)
                aTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGc                                                                                                                                 >  1:1086553/1‑150 (MQ=255)
                aTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGc                                                                                                                                 <  1:3689209/150‑1 (MQ=255)
                 tGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGca                                                                                                                                <  1:1579414/150‑1 (MQ=255)
                    tttGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGcacca                                                                                                                             >  2:1538757/1‑150 (MQ=255)
                          aaGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCAcgccgc                                                                                                                       >  1:4883148/1‑150 (MQ=255)
                          aaGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCAcgccgc                                                                                                                       >  1:3180632/1‑150 (MQ=255)
                           aGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACgccgcc                                                                                                                      >  1:4748738/1‑150 (MQ=255)
                             gCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCaa                                                                                                                    <  2:304541/150‑1 (MQ=255)
                              cAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAAt                                                                                                                   >  1:2808785/1‑150 (MQ=255)
                               aGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATc                                                                                                                  >  1:593380/1‑150 (MQ=255)
                               aGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCACCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATc                                                                                                                  <  2:2704970/150‑1 (MQ=255)
                                gTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCa                                                                                                                 <  2:3094751/150‑1 (MQ=255)
                                gTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCa                                                                                                                 <  1:322930/150‑1 (MQ=255)
                                gTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCa                                                                                                                 <  1:2472076/150‑1 (MQ=255)
                                       aCCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCgg                                                                                                          <  2:3251289/150‑1 (MQ=255)
                                       aCCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCgg                                                                                                          <  1:2512056/150‑1 (MQ=255)
                                            gATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACt                                                                                                     <  1:3505972/150‑1 (MQ=255)
                                                cAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGaa                                                                                                 <  1:5153534/150‑1 (MQ=255)
                                                 aGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAAc                                                                                                >  2:1836894/1‑150 (MQ=255)
                                                    ttGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGa                                                                                             >  1:4843962/1‑150 (MQ=255)
                                                        cGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCtt                                                                                         >  2:450521/1‑150 (MQ=255)
                                                        cGAGGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCtt                                                                                         <  2:5046547/150‑1 (MQ=255)
                                                         gATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTc                                                                                        <  2:1174354/150‑1 (MQ=255)
                                                         gATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGAAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTc                                                                                        >  2:2602440/1‑150 (MQ=255)
                                                          aTGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTcc                                                                                       >  1:5049050/1‑150 (MQ=255)
                                                            ggCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTg                                                                                     <  2:4212788/150‑1 (MQ=255)
                                                              cTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAg                                                                                   >  2:318291/1‑150 (MQ=255)
                                                              cTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAg                                                                                   >  1:1032026/1‑150 (MQ=255)
                                                                 catcatCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCa                                                                                <  1:2292567/150‑1 (MQ=255)
                                                                  atcatcCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACAGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAg                                                                               >  2:3493830/1‑150 (MQ=255)
                                                                    catcCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGtt                                                                             >  1:673529/1‑150 (MQ=255)
                                                                        ccACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTc                                                                         <  1:3728040/150‑1 (MQ=255)
                                                                              cttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGa                                                                   >  1:3503301/1‑150 (MQ=255)
                                                                                tgcccgtGCCCAGGTAGAAGTCGCCGAGGGTTTTACCATCCAGCACCAGGTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACa                                                                 <  1:2683768/144‑1 (MQ=255)
                                                                                tGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACa                                                                 <  2:50154/150‑1 (MQ=255)
                                                                                tGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACa                                                                 <  1:703282/150‑1 (MQ=255)
                                                                                tGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACa                                                                 <  1:86391/150‑1 (MQ=255)
                                                                                 gATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAg                                                                <  2:3071286/150‑1 (MQ=255)
                                                                                     ttGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAgcggc                                                            <  1:199518/150‑1 (MQ=255)
                                                                                         ccAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGt                                                        >  1:1200266/1‑150 (MQ=255)
                                                                                           aGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCa                                                      >  1:1265519/1‑150 (MQ=255)
                                                                                           aGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCa                                                      >  1:1664964/1‑150 (MQ=255)
                                                                                           aGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCa                                                      <  1:15239/150‑1 (MQ=255)
                                                                                           aGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCa                                                      >  1:2328323/1‑150 (MQ=255)
                                                                                              tAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCagag                                                   <  1:1892867/150‑1 (MQ=255)
                                                                                               aGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAgg                                                  >  2:4230629/1‑150 (MQ=255)
                                                                                                gATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGc                                                 <  2:2271455/150‑1 (MQ=255)
                                                                                                 aTGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGACCTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCa                                                >  2:847698/1‑150 (MQ=255)
                                                                                                 aTGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCa                                                >  2:4772030/1‑150 (MQ=255)
                                                                                                  tGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCAc                                               >  2:1155642/1‑150 (MQ=255)
                                                                                                   gTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCAcc                                              <  1:296435/150‑1 (MQ=255)
                                                                                                      gCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCaaa                                           <  2:4207178/150‑1 (MQ=255)
                                                                                                        cGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAAt                                         <  1:1873223/150‑1 (MQ=255)
                                                                                                        cGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAAt                                         >  1:4175014/1‑150 (MQ=255)
                                                                                                        cGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAAt                                         >  1:69057/1‑150 (MQ=255)
                                                                                                         gAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATc                                        <  2:1327914/150‑1 (MQ=255)
                                                                                                              ttttACCATCCAGCACCAGTTCGCCAGACTTCAGCCCATGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACgg                                   >  1:2575764/1‑150 (MQ=255)
                                                                                                              ttttACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACgg                                   >  1:4349812/1‑150 (MQ=255)
                                                                                                              ttttACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACgg                                   >  1:1294852/1‑150 (MQ=255)
                                                                                                                    cATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAg                             <  1:3594870/150‑1 (MQ=255)
                                                                                                                    cATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAg                             <  2:1260061/150‑1 (MQ=255)
                                                                                                                      tCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCg                           <  2:369713/150‑1 (MQ=255)
                                                                                                                       ccAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGa                          >  1:4119297/1‑150 (MQ=255)
                                                                                                                             ccAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCt                    >  1:5015297/1‑150 (MQ=255)
                                                                                                                             ccAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCt                    <  1:1538757/150‑1 (MQ=255)
                                                                                                                             ccAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCt                    <  2:4299907/150‑1 (MQ=255)
                                                                                                                             ccAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCt                    >  2:1103157/1‑150 (MQ=255)
                                                                                                                              cAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTg                   >  2:4183743/1‑150 (MQ=255)
                                                                                                                               aGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGt                  <  1:3301082/150‑1 (MQ=255)
                                                                                                                                gTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGtt                 <  1:3455918/150‑1 (MQ=255)
                                                                                                                                 ttCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGttt                <  1:5075749/150‑1 (MQ=255)
                                                                                                                                  tCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGtttt               <  2:4616018/150‑1 (MQ=255)
                                                                                                                                   cGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTa              <  1:3906011/150‑1 (MQ=255)
                                                                                                                                    gCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTAc             <  2:3414841/150‑1 (MQ=255)
                                                                                                                                     ccAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACg            >  2:3210418/1‑150 (MQ=255)
                                                                                                                                     ccAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACg            >  2:2788986/1‑150 (MQ=255)
                                                                                                                                      cAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGc           >  2:1020864/1‑150 (MQ=255)
                                                                                                                                      cAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGc           >  2:2310348/1‑150 (MQ=255)
                                                                                                                                        gACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCa         >  1:4474523/1‑150 (MQ=255)
                                                                                                                                        gACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCa         >  2:3590394/1‑150 (MQ=255)
                                                                                                                                        gACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCa         >  1:606141/1‑150 (MQ=255)
                                                                                                                                        gACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCa         >  2:2487322/1‑150 (MQ=255)
                                                                                                                                         aCTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCAc        <  2:2381310/150‑1 (MQ=255)
                                                                                                                                           ttCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCg      <  2:4410596/150‑1 (MQ=255)
                                                                                                                                             cAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGga    >  1:1423188/1‑150 (MQ=255)
                                                                                                                                             cAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGga    <  2:1919464/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCGGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAAG  >  CP009273/3904406‑3904698

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: