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.4% 600.6 / ‑0.7 183P111P (CCA→CCTpstSperiplasmic phosphate binding protein, high‑affinity
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (93/86);  total (95/88)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 183 total aligned reads displayed.

GGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAA  >  CP009273/3904408‑3904697
                                                                                                                                                 |                                                                                                                                                
ggcgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAgg                                                                                                                                                <  2:1160498/148‑1 (MQ=255)
    tactacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATg                                                                                                                                            <  2:2739151/148‑1 (MQ=255)
    tactacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATg                                                                                                                                            <  1:2998982/148‑1 (MQ=255)
    tactacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATg                                                                                                                                            >  1:1533037/1‑148 (MQ=255)
    tactacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATg                                                                                                                                            >  1:615893/1‑148 (MQ=255)
         cAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAAc                                                                                                                                       >  2:1064540/1‑148 (MQ=255)
          aGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAAcc                                                                                                                                      >  2:3181779/1‑148 (MQ=255)
              aTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCa                                                                                                                                  >  1:2464869/1‑148 (MQ=255)
              aTGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCa                                                                                                                                  >  2:339495/1‑148 (MQ=255)
                  tttGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGcac                                                                                                                              >  2:1074862/1‑148 (MQ=255)
                    tgtgAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGcacca                                                                                                                            <  2:1677850/148‑1 (MQ=255)
                    tgtgAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGcacca                                                                                                                            <  2:1753271/148‑1 (MQ=255)
                         aGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCAcgccg                                                                                                                       >  1:1310870/1‑148 (MQ=255)
                         aGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCAcgccg                                                                                                                       >  1:1765780/1‑148 (MQ=255)
                          ggCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCAcgccgc                                                                                                                      >  2:559009/1‑148 (MQ=255)
                              gTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAAt                                                                                                                  <  2:710824/148‑1 (MQ=255)
                              gTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAAt                                                                                                                  <  2:1135646/148‑1 (MQ=255)
                                ttCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCa                                                                                                                >  1:1457622/1‑148 (MQ=255)
                                ttCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCa                                                                                                                <  2:1386396/148‑1 (MQ=255)
                                    gACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGt                                                                                                            >  1:3417524/1‑148 (MQ=255)
                                    gACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGt                                                                                                            <  2:1887633/148‑1 (MQ=255)
                                     aCCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTc                                                                                                           >  1:2436077/1‑148 (MQ=255)
                                      cccGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCg                                                                                                          <  2:715793/148‑1 (MQ=255)
                                       ccGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCgg                                                                                                         >  1:1080390/1‑148 (MQ=255)
                                       ccGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCgg                                                                                                         >  2:2579579/1‑148 (MQ=255)
                                          gATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGaa                                                                                                      <  1:2653818/148‑1 (MQ=255)
                                          gATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGaa                                                                                                      <  2:2984189/148‑1 (MQ=255)
                                          gATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGaa                                                                                                      <  2:1561921/148‑1 (MQ=255)
                                          gATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGaa                                                                                                      <  1:2892650/148‑1 (MQ=255)
                                           aTTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAAc                                                                                                     >  2:2161433/1‑148 (MQ=255)
                                             tCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTg                                                                                                   >  1:1196269/1‑148 (MQ=255)
                                                gTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGaa                                                                                                <  1:2662950/148‑1 (MQ=255)
                                                 tttGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAAc                                                                                               >  2:1681372/1‑148 (MQ=255)
                                                    ggCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGa                                                                                            <  2:751264/148‑1 (MQ=255)
                                                     gCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGAc                                                                                           >  2:497829/1‑148 (MQ=255)
                                                      cGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGAcc                                                                                          >  1:3497642/1‑148 (MQ=255)
                                                        aTGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCtt                                                                                        <  2:1661230/148‑1 (MQ=255)
                                                        aTGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCtt                                                                                        >  1:760807/1‑148 (MQ=255)
                                                        aTGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCtt                                                                                        >  1:2519868/1‑148 (MQ=255)
                                                         tGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTc                                                                                       >  2:329767/1‑148 (MQ=255)
                                                          ggCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTcc                                                                                      <  2:117727/148‑1 (MQ=255)
                                                          ggCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTcc                                                                                      >  1:431624/1‑148 (MQ=255)
                                                            cTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTg                                                                                    >  2:2651017/1‑148 (MQ=255)
                                                              tcatcaTCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAg                                                                                  >  1:205239/1‑148 (MQ=255)
                                                               catcatCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGc                                                                                 <  2:1645756/148‑1 (MQ=255)
                                                               catcatCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGc                                                                                 <  2:1646019/148‑1 (MQ=255)
                                                               catcatCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGc                                                                                 >  1:3253296/1‑148 (MQ=255)
                                                               catcatCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGc                                                                                 <  1:3093526/148‑1 (MQ=255)
                                                                atcatcCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGcc                                                                                >  2:1495072/1‑148 (MQ=255)
                                                                atcatcCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGcc                                                                                >  1:3381071/1‑148 (MQ=255)
                                                                 tcatcCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCa                                                                               >  1:1927468/1‑148 (MQ=255)
                                                                    tcCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGtt                                                                            >  1:3507668/1‑148 (MQ=255)
                                                                     cccACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGttt                                                                           <  2:1320829/148‑1 (MQ=255)
                                                                         cttcttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCg                                                                       >  2:845623/1‑148 (MQ=255)
                                                                         cttcttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACAATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCg                                                                       >  2:2539401/1‑148 (MQ=255)
                                                                          ttcttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGt                                                                      <  1:543956/148‑1 (MQ=255)
                                                                               gATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGAc                                                                 <  1:274599/148‑1 (MQ=255)
                                                                                 ttttGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAg                                                               >  2:3026213/1‑148 (MQ=255)
                                                                                 ttttGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAg                                                               >  1:3460501/1‑148 (MQ=255)
                                                                                 ttttGCCCAGGTAGATGTCGCCGAGGGTTTTACCAACCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCCCCAGCACCACGCCGCCAATCACGGTCGGGACCTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAg                                                               >  2:3026200/1‑148 (MQ=255)
                                                                                  tttGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGACTGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAgc                                                              >  1:1423415/1‑148 (MQ=255)
                                                                                  tttGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAgc                                                              >  2:2372337/1‑148 (MQ=255)
                                                                                   ttGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAgcg                                                             <  1:1056591/148‑1 (MQ=255)
                                                                                     gCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAgcggc                                                           >  2:3487995/1‑148 (MQ=255)
                                                                                            tAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCag                                                    <  1:3276393/148‑1 (MQ=255)
                                                                                               aTGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAgg                                                 >  1:1417151/1‑148 (MQ=255)
                                                                                               aTGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAgg                                                 >  2:2013708/1‑148 (MQ=255)
                                                                                                 gTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCa                                               <  2:409672/148‑1 (MQ=255)
                                                                                                  tCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCAc                                              <  1:842255/148‑1 (MQ=255)
                                                                                                   cGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCAcc                                             >  2:2838266/1‑148 (MQ=255)
                                                                                                      cGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCaaa                                          <  1:737546/148‑1 (MQ=255)
                                                                                                       gAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCaaaa                                         >  1:1393736/1‑148 (MQ=255)
                                                                                                          ggTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCa                                      >  1:986813/1‑148 (MQ=255)
                                                                                                           gTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCaa                                     <  1:2101655/148‑1 (MQ=255)
                                                                                                            ttttACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAAc                                    <  1:1153653/148‑1 (MQ=255)
                                                                                                            ttttACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAAc                                    <  1:1365452/148‑1 (MQ=255)
                                                                                                             tttACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACg                                   >  1:1213196/1‑148 (MQ=255)
                                                                                                                aCCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTa                                >  1:2344018/1‑148 (MQ=255)
                                                                                                                aCCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTa                                <  2:351361/148‑1 (MQ=255)
                                                                                                                 ccATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTAt                               <  2:3058457/148‑1 (MQ=255)
                                                                                                                  cATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTAtt                              <  1:134175/148‑1 (MQ=255)
                                                                                                                         cACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATa                       >  1:3428080/1‑148 (MQ=255)
                                                                                                                            cAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATc                    <  1:1681372/148‑1 (MQ=255)
                                                                                                                            cAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATc                    >  1:2885798/1‑148 (MQ=255)
                                                                                                                               ttCGCCAGACTTCAGCCCAGGACTGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCCAAATCAACGGTATTAGCGATAATCTgg                 >  2:344435/1‑147 (MQ=255)
                                                                                                                               ttCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGt                 <  2:2231960/148‑1 (MQ=255)
                                                                                                                                tCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCGTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATGAGCGATAATCTGtt                <  2:1310870/148‑1 (MQ=255)
                                                                                                                                 cGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGttt               <  1:2161433/148‑1 (MQ=255)
                                                                                                                                  gCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGtttt              <  2:531454/148‑1 (MQ=255)
                                                                                                                                   ccAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTa             <  2:2464869/148‑1 (MQ=255)
                                                                                                                                    cAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTAc            >  2:1448116/1‑148 (MQ=255)
                                                                                                                                     aGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACg           >  1:2595522/1‑148 (MQ=255)
                                                                                                                                     aGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACg           >  1:3260104/1‑148 (MQ=255)
                                                                                                                                      gACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGc          <  1:2224639/148‑1 (MQ=255)
                                                                                                                                      gACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGc          >  1:375281/1‑148 (MQ=255)
                                                                                                                                       aCTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGcc         >  1:1872779/1‑148 (MQ=255)
                                                                                                                                        cTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCa        <  1:2961763/148‑1 (MQ=255)
                                                                                                                                             gCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGga   <  2:615893/148‑1 (MQ=255)
                                                                                                                                             gCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGga   <  2:1533037/148‑1 (MQ=255)
                                                                                                                                              cccAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGgaa  >  1:939016/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAA  >  CP009273/3904408‑3904697

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: