Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,424,234 C→A 83.2% S310* (TCA→TAA)  stfR → Rac prophage; putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,424,2340CA83.2% 637.4 / 145.1 273S310* (TCA→TAA) stfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base C (27/19);  new base A (125/102);  total (152/121)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 281 total aligned reads displayed.

AGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTCCGCAGCG  >  CP009273/1424085‑1424379
                                                                                                                                                     |                                                                                                                                                 
   ttCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACTAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCTTCTGCCAGTGCAGCGTcaac                                                                                                                                                >  1:1446483/1‑150 (MQ=255)
     ccTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTaaacaa                                                                                                                                              >  1:295158/1‑150 (MQ=255)
     ccTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTaaacaa                                                                                                                                              >  2:445521/1‑150 (MQ=255)
     ccTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTaaacaa                                                                                                                                              >  1:1734241/1‑150 (MQ=255)
      cTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAg                                                                                                                                             >  2:1959098/1‑150 (MQ=255)
      cTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAg                                                                                                                                             >  2:4046445/1‑150 (MQ=255)
      cTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAg                                                                                                                                             >  2:3606396/1‑150 (MQ=255)
           cAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCg                                                                                                                                        <  2:1523530/150‑1 (MQ=255)
           cAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCg                                                                                                                                        <  2:3768898/150‑1 (MQ=255)
            aaCGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCgg                                                                                                                                       >  2:3426937/1‑150 (MQ=255)
                gcggcAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAg                                                                                                                                   <  2:1364680/150‑1 (MQ=255)
                gcggcAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAg                                                                                                                                   <  2:2190258/150‑1 (MQ=255)
                    cAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCt                                                                                                                               <  1:1839763/150‑1 (MQ=255)
                     aGGAAATTCCGCGATGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTc                                                                                                                              <  1:3721114/150‑1 (MQ=255)
                      ggAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCa                                                                                                                             >  1:3513201/1‑150 (MQ=255)
                       gAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAg                                                                                                                            >  1:2652701/1‑150 (MQ=255)
                       gAAAGTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAg                                                                                                                            <  2:3956271/150‑1 (MQ=255)
                        aaaTTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGc                                                                                                                           >  1:794358/1‑150 (MQ=255)
                         aaTTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGcc                                                                                                                          >  2:2502481/1‑150 (MQ=255)
                         aaTTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGcc                                                                                                                          >  1:1065009/1‑150 (MQ=255)
                             ccGCGACGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTg                                                                                                                      <  1:2392503/150‑1 (MQ=255)
                              cgcgAAGGCGGCAAAACCGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGc                                                                                                                     <  1:1636014/150‑1 (MQ=255)
                                    agcggcAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCAccgc                                                                                                               >  2:1583383/2‑150 (MQ=255)
                                       ggcAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCAccgccgc                                                                                                            >  2:1952034/1‑150 (MQ=255)
                                         cAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCAccgccgccg                                                                                                          >  1:3680245/1‑150 (MQ=255)
                                          aaaaaCGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCgg                                                                                                         <  1:2498608/150‑1 (MQ=255)
                                               cGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAAt                                                                                                    >  2:2915605/1‑150 (MQ=255)
                                                gTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACCGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATc                                                                                                   <  1:3415987/150‑1 (MQ=255)
                                                gTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGACAATc                                                                                                   >  2:1450556/1‑150 (MQ=255)
                                                gTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATc                                                                                                   >  2:1392063/1‑150 (MQ=255)
                                                 tCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCg                                                                                                  <  1:4047184/150‑1 (MQ=255)
                                                   cGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGc                                                                                                >  1:564967/1‑150 (MQ=255)
                                                   cGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGc                                                                                                >  2:665439/1‑150 (MQ=255)
                                                   cGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGc                                                                                                <  2:188214/150‑1 (MQ=255)
                                                    gagaCGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCa                                                                                               <  2:1429108/150‑1 (MQ=255)
                                                    gagaCGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCa                                                                                               >  2:4512284/1‑150 (MQ=255)
                                                     agaCGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAg                                                                                              >  1:2618672/1‑150 (MQ=255)
                                                      gaCGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGa                                                                                             >  1:430496/1‑150 (MQ=255)
                                                         gAACGCCTGGTCTTCTGAATCGGCAGCGGGACTGAGCGCCTCGGCTGCGGCTGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCTTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGGAGAAAg                                                                                          <  2:124214/150‑1 (MQ=255)
                                                         gAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAg                                                                                          <  1:1009177/150‑1 (MQ=255)
                                                         gAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAAAAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAg                                                                                          >  2:298828/1‑150 (MQ=255)
                                                           aCGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAgcg                                                                                        <  1:109056/150‑1 (MQ=255)
                                                           aCGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAgcg                                                                                        <  2:27022/150‑1 (MQ=255)
                                                           aCGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAgcg                                                                                        <  2:2645396/150‑1 (MQ=255)
                                                             gCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGcgcc                                                                                      >  1:4183336/1‑150 (MQ=255)
                                                               cAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGcgccgc                                                                                    >  1:4267076/1‑150 (MQ=255)
                                                                aGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCa                                                                                   >  1:2064956/1‑150 (MQ=255)
                                                                 ggTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCAt                                                                                  >  2:723329/1‑150 (MQ=255)
                                                                  gTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATc                                                                                 >  2:4287844/1‑150 (MQ=255)
                                                                  gTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATc                                                                                 >  2:2014038/1‑150 (MQ=255)
                                                                         gAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCt                                                                          <  2:1729081/150‑1 (MQ=255)
                                                                         gAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCt                                                                          >  2:3714558/1‑150 (MQ=255)
                                                                          aaaCGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCtt                                                                         >  2:101197/1‑150 (MQ=255)
                                                                             cGGCTGGGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACCGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCaa                                                                      <  1:4341298/150‑1 (MQ=255)
                                                                              ggCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAAc                                                                     <  2:3649217/150‑1 (MQ=255)
                                                                                    gggACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCcac                                                               <  2:1037849/150‑1 (MQ=255)
                                                                                    gggACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCcac                                                               >  1:997559/1‑150 (MQ=255)
                                                                                    gggACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCcac                                                               >  1:2409743/1‑150 (MQ=255)
                                                                                    gggACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCcac                                                               <  1:4415314/150‑1 (MQ=255)
                                                                                    gggACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCcac                                                               >  1:2560781/1‑150 (MQ=255)
                                                                                     ggACAGCGCGCCTCGGCTGCGGCTGGCTCAAAAACCGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCGGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCcaca                                                              <  1:2683241/150‑1 (MQ=255)
                                                                                     ggACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCcaca                                                              >  1:1314487/1‑150 (MQ=255)
                                                                                     ggACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCcaca                                                              >  1:3599323/1‑150 (MQ=255)
                                                                                         agagCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGa                                                          >  1:2097296/1‑150 (MQ=255)
                                                                                         agagCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGa                                                          >  2:3361061/1‑150 (MQ=255)
                                                                                             cgcCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGc                                                      >  2:2512411/1‑150 (MQ=255)
                                                                                               ccTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTg                                                    >  1:357556/1‑150 (MQ=255)
                                                                                                 tCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCACCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGc                                                  >  1:3607953/1‑150 (MQ=255)
                                                                                                        cGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCa                                           >  1:149518/1‑150 (MQ=255)
                                                                                                         ggCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCAc                                          >  2:2560479/1‑150 (MQ=255)
                                                                                                              gCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAAc                                     <  1:3189773/150‑1 (MQ=255)
                                                                                                                 cAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAgg                                  <  2:3177579/150‑1 (MQ=255)
                                                                                                                 cAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAgg                                  <  2:3263783/150‑1 (MQ=255)
                                                                                                                 cAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAgg                                  <  2:1158414/150‑1 (MQ=255)
                                                                                                                  aaaaaCAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGc                                 >  2:82199/1‑150 (MQ=255)
                                                                                                                   aaaaCAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGcc                                >  2:2866276/1‑150 (MQ=255)
                                                                                                                   aaaaCAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGcc                                <  1:2736578/150‑1 (MQ=255)
                                                                                                                    aaaCAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCa                               <  1:1141683/150‑1 (MQ=255)
                                                                                                                     aaCAGCGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAg                              >  1:2228297/1‑150 (MQ=255)
                                                                                                                          cGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCgcag                         >  1:1483284/1‑150 (MQ=255)
                                                                                                                          cGGCTGCGTCGTCTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCgcag                         >  1:3929037/1‑150 (MQ=255)
                                                                                                                             cTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGcagcag                      >  2:4428068/1‑150 (MQ=255)
                                                                                                                                 gtcgtcTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAg                  <  2:223444/150‑1 (MQ=255)
                                                                                                                                 gtcgtcTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAg                  >  1:2240427/1‑150 (MQ=255)
                                                                                                                                   cgtcTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGt                >  2:4098678/1‑150 (MQ=255)
                                                                                                                                   cgtcTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGt                >  2:204538/1‑150 (MQ=255)
                                                                                                                                    gtcTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTc               >  2:427482/1‑150 (MQ=255)
                                                                                                                                    gtcTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTc               >  1:145371/1‑150 (MQ=255)
                                                                                                                                     tcTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCt              >  2:3829463/1‑150 (MQ=255)
                                                                                                                                     tcTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCt              <  1:3832935/150‑1 (MQ=255)
                                                                                                                                      cTGCCAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTg             <  2:4237658/150‑1 (MQ=255)
                                                                                                                                          cAGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTc         <  1:3866677/150‑1 (MQ=255)
                                                                                                                                           aGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTcc        >  1:3050476/1‑150 (MQ=255)
                                                                                                                                           aGTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTcc        <  1:3606396/150‑1 (MQ=255)
                                                                                                                                            gTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGACAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTCCg       <  2:1628249/150‑1 (MQ=255)
                                                                                                                                            gTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTCCg       <  1:378947/150‑1 (MQ=255)
                                                                                                                                            gTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTCCg       <  2:1476799/150‑1 (MQ=255)
                                                                                                                                            gTGCAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTCCg       <  2:3102050/150‑1 (MQ=255)
                                                                                                                                               cAGCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTCCGCAg    <  2:1382404/150‑1 (MQ=255)
                                                                                                                                                 gCGTAAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTCCGCAGCg  >  1:2150942/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGTCTGCCAGTGCAGCGTCAACAAGTGCCGGGCAGGCCTCAGCCAGTGCCACCGCCGCCGGAAAATCGGCAGAAAGCGCCGCATCGTCTGCTTCAACAGCCACAACGAAGGCTGGCGAAGCCACTGAACAGGCCAGCGCAGCAGCGAGGTCTGCTTCCGCAGCG  >  CP009273/1424085‑1424379

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: ‑