Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,904,553 T→A 94.8% 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,5530TA94.8% 534.0 / 15.2 172P111P (CCA→CCTpstSperiplasmic phosphate binding protein, high‑affinity
Reads supporting (aligned to +/- strand):  ref base T (4/5);  new base A (94/68);  total (98/74)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 172 total aligned reads displayed.

GCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAA  >  CP009273/3904409‑3904697
                                                                                                                                                |                                                                                                                                                
gcgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGa                                                                                                                                               <  2:246246/148‑1 (MQ=255)
gcgTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGa                                                                                                                                               <  1:1994040/148‑1 (MQ=255)
    actacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGt                                                                                                                                           <  1:3212705/148‑1 (MQ=255)
     ctacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGtt                                                                                                                                          >  1:359852/1‑148 (MQ=255)
      tacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTa                                                                                                                                         <  2:3308136/148‑1 (MQ=255)
      tacAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTa                                                                                                                                         >  2:610097/1‑148 (MQ=255)
       acAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTaa                                                                                                                                        <  2:2775069/148‑1 (MQ=255)
           cAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGc                                                                                                                                    <  2:3359371/148‑1 (MQ=255)
                 tttGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGcac                                                                                                                              <  1:705199/148‑1 (MQ=255)
                  ttGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGcacc                                                                                                                             >  1:3291827/1‑148 (MQ=255)
                      gAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCAcgc                                                                                                                         >  2:3320843/1‑148 (MQ=255)
                      gAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCAcgc                                                                                                                         >  1:3237483/1‑148 (MQ=255)
                         ggCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCAcgccgc                                                                                                                      >  2:241034/1‑148 (MQ=255)
                         ggCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCAcgccgc                                                                                                                      <  1:1893434/148‑1 (MQ=255)
                         ggCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCAcgccgc                                                                                                                      >  2:837153/1‑148 (MQ=255)
                           cAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCa                                                                                                                    <  1:1140299/148‑1 (MQ=255)
                           cAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCa                                                                                                                    >  1:1495069/1‑148 (MQ=255)
                               ttCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCa                                                                                                                <  2:311299/148‑1 (MQ=255)
                                  aGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACgg                                                                                                             <  2:1736565/148‑1 (MQ=255)
                                    aCCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTc                                                                                                           >  1:2318592/1‑148 (MQ=255)
                                       cGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCggg                                                                                                        >  1:1902864/1‑148 (MQ=255)
                                        ggATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGa                                                                                                       <  2:157972/148‑1 (MQ=255)
                                          aTTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAAc                                                                                                     <  2:1902864/148‑1 (MQ=255)
                                             cAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTCGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTgg                                                                                                  >  2:1889419/1‑148 (MQ=255)
                                             cAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTgg                                                                                                  >  1:2315519/1‑148 (MQ=255)
                                               gTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGaa                                                                                                <  1:2057256/148‑1 (MQ=255)
                                               gTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGaa                                                                                                <  1:1415930/148‑1 (MQ=255)
                                                tttGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAAc                                                                                               >  2:679940/1‑148 (MQ=255)
                                                     cGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCATTTCGCCATACTTCAGCCCATGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGAcc                                                                                          <  1:35080/148‑1 (MQ=255)
                                                     cGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGAcc                                                                                          <  2:2315519/148‑1 (MQ=255)
                                                         ggCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTcc                                                                                      <  1:405844/148‑1 (MQ=255)
                                                          gCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCt                                                                                     >  1:3030991/1‑148 (MQ=255)
                                                          gCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCt                                                                                     >  1:3031027/1‑148 (MQ=255)
                                                           cTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTg                                                                                    >  1:1879869/1‑148 (MQ=255)
                                                             tcatcaTCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAg                                                                                  <  2:3026314/148‑1 (MQ=255)
                                                             tcatcaTCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAg                                                                                  >  2:2182073/1‑148 (MQ=255)
                                                             tcatcaTCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAg                                                                                  >  1:566159/1‑148 (MQ=255)
                                                             tcatcaTCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAg                                                                                  >  2:2182076/1‑148 (MQ=255)
                                                                 catcCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAg                                                                              >  2:3334147/1‑148 (MQ=255)
                                                                   tcCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGtt                                                                            >  1:588349/1‑148 (MQ=255)
                                                                    cccACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGttt                                                                           >  1:695414/1‑148 (MQ=255)
                                                                    cccACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGttt                                                                           <  2:434898/148‑1 (MQ=255)
                                                                    cccACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGttt                                                                           <  2:434920/148‑1 (MQ=255)
                                                                     ccACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGtttt                                                                          >  2:2436695/1‑148 (MQ=255)
                                                                      cACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGttttt                                                                         >  2:955575/1‑148 (MQ=255)
                                                                         ttcttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGt                                                                      >  1:946062/1‑148 (MQ=255)
                                                                          tcttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTc                                                                     >  1:3310996/1‑148 (MQ=255)
                                                                          tcttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTc                                                                     <  1:1526232/148‑1 (MQ=255)
                                                                            ttGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTTAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAg                                                                   <  2:1227550/148‑1 (MQ=255)
                                                                             tGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGa                                                                  <  1:2182073/148‑1 (MQ=255)
                                                                             tGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGa                                                                  >  1:838181/1‑148 (MQ=255)
                                                                               aTTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACa                                                                >  2:131458/1‑148 (MQ=255)
                                                                                ttttGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAg                                                               <  1:1244944/148‑1 (MQ=255)
                                                                                 tttGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAgc                                                              >  1:2974765/1‑148 (MQ=255)
                                                                                       cAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGcgcg                                                        >  1:1081399/1‑148 (MQ=255)
                                                                                       cAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGcgcg                                                        <  1:253586/148‑1 (MQ=255)
                                                                                         ggTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTc                                                      <  2:2129218/148‑1 (MQ=255)
                                                                                         ggTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTc                                                      >  1:1618486/1‑148 (MQ=255)
                                                                                          gTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCa                                                     >  1:2648014/1‑148 (MQ=255)
                                                                                           tAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCag                                                    >  1:2434591/1‑148 (MQ=255)
                                                                                           tAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCag                                                    >  1:124755/1‑148 (MQ=255)
                                                                                            aGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCaga                                                   <  1:1637137/148‑1 (MQ=255)
                                                                                              aTGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAgg                                                 >  1:2354828/1‑148 (MQ=255)
                                                                                              aTGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAgg                                                 >  1:3377466/1‑148 (MQ=255)
                                                                                                gTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCa                                               >  1:1715039/1‑148 (MQ=255)
                                                                                                     cGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCaaa                                          >  1:1052118/1‑148 (MQ=255)
                                                                                                       aGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAAt                                        >  2:1090733/1‑148 (MQ=255)
                                                                                                         ggTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCa                                      <  1:1926368/148‑1 (MQ=255)
                                                                                                         ggTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCa                                      <  1:2167175/148‑1 (MQ=255)
                                                                                                           ttttACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAAc                                    >  2:2690265/1‑148 (MQ=255)
                                                                                                            tttACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACg                                   <  1:653250/148‑1 (MQ=255)
                                                                                                             ttACCATCCAGCACCAGTTCGCCAGACTTCAGCCCGGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACgg                                  <  2:2549274/148‑1 (MQ=255)
                                                                                                             ttACCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGCCAGAGGCACCAAAATCAACgg                                  >  2:1038371/1‑148 (MQ=255)
                                                                                                               aCCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTa                                >  2:2208728/1‑148 (MQ=255)
                                                                                                               aCCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTa                                >  2:551400/1‑148 (MQ=255)
                                                                                                               aCCATCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTa                                <  2:1213948/148‑1 (MQ=255)
                                                                                                                  aTCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTa                             >  2:1990984/1‑148 (MQ=255)
                                                                                                                   tCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAg                            <  1:720333/148‑1 (MQ=255)
                                                                                                                   tCCAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAg                            <  1:720315/148‑1 (MQ=255)
                                                                                                                     cAGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCg                          <  1:2581844/148‑1 (MQ=255)
                                                                                                                      aGCACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGa                         >  1:2820394/1‑148 (MQ=255)
                                                                                                                        cACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATa                       >  2:937721/1‑148 (MQ=255)
                                                                                                                        cACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATa                       >  2:33554/1‑148 (MQ=255)
                                                                                                                        cACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATa                       <  2:250346/148‑1 (MQ=255)
                                                                                                                        cACCAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATa                       <  2:1401409/148‑1 (MQ=255)
                                                                                                                          ccAGTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAAt                     >  1:1729855/1‑148 (MQ=255)
                                                                                                                             gTTCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTg                  >  2:2468786/1‑148 (MQ=255)
                                                                                                                              ttCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGt                 <  1:993133/148‑1 (MQ=255)
                                                                                                                              ttCGCCAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGt                 <  1:2133476/148‑1 (MQ=255)
                                                                                                                                   cAGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTAc            >  2:3097808/1‑148 (MQ=255)
                                                                                                                                    aGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACg           >  1:3370443/1‑148 (MQ=255)
                                                                                                                                    aGACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACg           >  1:3144317/1‑148 (MQ=255)
                                                                                                                                     gACTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAACTCAACGGTATTAGCGATATTCTGTTTTACGc          >  2:2964820/1‑148 (MQ=255)
                                                                                                                                      aCTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGcc         <  1:2850456/148‑1 (MQ=255)
                                                                                                                                      aCTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGcc         <  1:1542103/148‑1 (MQ=255)
                                                                                                                                      aCTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGcc         <  1:1542006/148‑1 (MQ=255)
                                                                                                                                      aCTTCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGcc         <  1:1446402/148‑1 (MQ=255)
                                                                                                                                        gtCAGCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCAc       <  2:161859/147‑1 (MQ=255)
                                                                                                                                            gCCCAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGga   >  1:1075923/1‑148 (MQ=255)
                                                                                                                                             cccAGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGgaa  <  1:3450860/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GCGTACTACAGCAATGTTTTGTGAAGGCAGTTTCAGACCCGGATTCAGTTTGGCGATGGCTTCATCATCCCACTTCTTGATTTTGCCCAGGTAGATGTCGCCGAGGGTTTTACCATCCAGCACCAGTTCGCCAGACTTCAGCCCTGGAATGTTAACCGCCAGCACCACGCCGCCAATCACGGTCGGGAACTGGAACAGACCTTCCTGAGCCAGTTTTTCGTCAGACAGCGGCGCGTCAGAGGCACCAAAATCAACGGTATTAGCGATAATCTGTTTTACGCCACCGGAA  >  CP009273/3904409‑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: