Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 397,959 G→A 100% A253T (GCG→ACG)  phoA → bacterial alkaline phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273397,9590GA98.7% 751.4 / 0.0 229A253T (GCG→ACG) phoAbacterial alkaline phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (126/100);  total (126/103)
Fisher's exact test for biased strand distribution p-value = 8.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 247 total aligned reads displayed.

CTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAATA  >  CP009273/397815‑398103
                                                                                                                                                |                                                                                                                                                
    aaaaaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTgg                                                                                                                                           <  2:1420257/148‑1 (MQ=255)
    aaaaaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTgg                                                                                                                                           <  1:832043/148‑1 (MQ=255)
       aaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTa                                                                                                                                        <  1:3162443/148‑1 (MQ=255)
        aGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTAt                                                                                                                                       >  1:2775939/1‑148 (MQ=255)
          gCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCa                                                                                                                                     >  2:3174236/1‑148 (MQ=255)
             gAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGtt                                                                                                                                  <  2:247099/148‑1 (MQ=255)
              aaaaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTg                                                                                                                                 <  2:2183366/148‑1 (MQ=255)
              aaaaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTg                                                                                                                                 >  1:1899051/1‑148 (MQ=255)
              aaaaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTg                                                                                                                                 <  2:1108707/148‑1 (MQ=255)
              aaaaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTg                                                                                                                                 >  2:1507381/1‑148 (MQ=255)
                aaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGt                                                                                                                               >  2:2585674/1‑148 (MQ=255)
                aaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGt                                                                                                                               <  1:1174975/148‑1 (MQ=255)
                 aGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTg                                                                                                                              >  1:954113/1‑148 (MQ=255)
                  ggATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGa                                                                                                                             <  1:3461006/148‑1 (MQ=255)
                   gATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAg                                                                                                                            >  1:1271049/1‑148 (MQ=255)
                     tCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCg                                                                                                                          >  1:681092/1‑148 (MQ=255)
                     tCGATGACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCg                                                                                                                          <  1:1618798/148‑1 (MQ=255)
                      cGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAACCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGa                                                                                                                         <  2:431193/148‑1 (MQ=255)
                      cGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGa                                                                                                                         <  2:130376/148‑1 (MQ=255)
                      cGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGa                                                                                                                         <  2:1905819/148‑1 (MQ=255)
                       gATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCGGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGag                                                                                                                        >  1:3048269/1‑147 (MQ=255)
                        aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACACTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtg                                                                                                                       <  1:3083886/148‑1 (MQ=255)
                            ccGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGACGCGGCAACCGCTGGTGCCTGGCAGGGACAAACGCTGCGTGAACCGGCACAGACGCGTGGGTCTCAGTTGGTGAGCGAggctgc                                                                                                                   >  2:3362694/1‑148 (MQ=255)
                            ccGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTTAGCGAtgctgc                                                                                                                   <  2:1029080/148‑1 (MQ=255)
                               aaCAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTc                                                                                                                >  2:2527618/1‑148 (MQ=255)
                               aaCAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTGATCAGTTGGTGAGCGATGCTGCCTc                                                                                                                >  1:2668632/1‑148 (MQ=255)
                                aCAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCa                                                                                                               >  2:2909870/1‑148 (MQ=255)
                                    ctgctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGa                                                                                                           >  2:3453822/1‑148 (MQ=255)
                                    ctgctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGa                                                                                                           <  1:2901421/148‑1 (MQ=255)
                                     tgctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGaa                                                                                                          <  2:1149810/148‑1 (MQ=255)
                                         tAACGCTCGTGCCGACGTTACGCTTGGCGGCTGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGGGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCg                                                                                                      >  1:2194004/1‑148 (MQ=255)
                                         tAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCg                                                                                                      >  1:2008192/1‑148 (MQ=255)
                                         tAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCg                                                                                                      >  1:2007161/1‑148 (MQ=255)
                                         tAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCg                                                                                                      >  1:1943910/1‑148 (MQ=255)
                                          aaCGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCgg                                                                                                     <  1:2808349/148‑1 (MQ=255)
                                          aaCGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCgg                                                                                                     <  1:2431344/148‑1 (MQ=255)
                                            cGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTg                                                                                                   >  1:2796385/1‑148 (MQ=255)
                                            cGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTg                                                                                                   >  1:3361809/1‑148 (MQ=255)
                                             gCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGa                                                                                                  <  1:2906568/148‑1 (MQ=255)
                                              cTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGACACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGAc                                                                                                 >  1:1490111/1‑148 (MQ=255)
                                                cGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACgg                                                                                               <  2:332680/148‑1 (MQ=255)
                                                 gTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGa                                                                                              <  1:471357/148‑1 (MQ=255)
                                                   gCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTGTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAg                                                                                            <  2:3150584/148‑1 (MQ=255)
                                                    ccGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGc                                                                                           >  2:1264940/1‑148 (MQ=255)
                                                      gACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGa                                                                                         <  2:2950688/148‑1 (MQ=255)
                                                       aCGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGaa                                                                                        <  1:235720/148‑1 (MQ=255)
                                                        cGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAGCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAAt                                                                                       >  2:1072909/1‑148 (MQ=255)
                                                         gTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATc                                                                                      <  2:1647745/148‑1 (MQ=255)
                                                             cGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCa                                                                                  >  1:1814054/1‑148 (MQ=255)
                                                               cTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCaaa                                                                                >  1:321853/1‑148 (MQ=255)
                                                               cTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCaaa                                                                                >  1:518217/1‑148 (MQ=255)
                                                               cTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCaaa                                                                                >  1:1091035/1‑148 (MQ=255)
                                                                 tGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCaaaaa                                                                              <  1:1662403/148‑1 (MQ=255)
                                                                 tGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCaaaaa                                                                              <  2:2080513/148‑1 (MQ=255)
                                                                     ggcggcgCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAAcccc                                                                          <  1:601004/148‑1 (MQ=255)
                                                                     ggcggcgCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAAcccc                                                                          <  1:1424857/148‑1 (MQ=255)
                                                                       cggcgCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTg                                                                        <  1:2073372/148‑1 (MQ=255)
                                                                         gcgcAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCt                                                                      >  1:2016493/1‑148 (MQ=255)
                                                                          cgcAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCtt                                                                     >  1:3368652/1‑148 (MQ=255)
                                                                              aaaaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGcc                                                                 >  2:1924781/1‑148 (MQ=255)
                                                                               aaaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAACAAAAACCCCTGCTTGGCCt                                                                >  1:237179/1‑148 (MQ=255)
                                                                                 aCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGt                                                              >  2:1863248/1‑148 (MQ=255)
                                                                                  ccTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGtt                                                             <  1:2057910/148‑1 (MQ=255)
                                                                                   cTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGTACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGttt                                                            >  2:2225255/1‑148 (MQ=255)
                                                                                   cTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGttt                                                            >  2:1807163/1‑148 (MQ=255)
                                                                                   cTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGttt                                                            >  1:2717127/1‑148 (MQ=255)
                                                                                    tttGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTg                                                           >  2:203112/1‑148 (MQ=255)
                                                                                    tttGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTg                                                           <  1:687908/148‑1 (MQ=255)
                                                                                            aaCGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGc                                                   >  2:1351235/1‑148 (MQ=255)
                                                                                             aCGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCa                                                  <  1:839983/148‑1 (MQ=255)
                                                                                                gCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAata                                               >  2:1183395/1‑148 (MQ=255)
                                                                                                   aCCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGc                                            >  2:397475/1‑148 (MQ=255)
                                                                                                     cGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCa                                          <  2:1668973/148‑1 (MQ=255)
                                                                                                          gTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTgcg                                     >  1:393077/1‑148 (MQ=255)
                                                                                                          gTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTgcg                                     <  1:1709295/148‑1 (MQ=255)
                                                                                                            gAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACTCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCt                                   >  1:695032/1‑148 (MQ=255)
                                                                                                            gAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCt                                   >  2:1521044/1‑148 (MQ=255)
                                                                                                               tGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGc                                <  1:1060703/148‑1 (MQ=255)
                                                                                                               tGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGc                                <  1:2814098/148‑1 (MQ=255)
                                                                                                                 gCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTa                              <  2:2678595/148‑1 (MQ=255)
                                                                                                                  cAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAg                             >  1:1320568/1‑148 (MQ=255)
                                                                                                                    gggAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGa                           >  1:622260/1‑148 (MQ=255)
                                                                                                                    gggAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGa                           >  1:3312596/1‑148 (MQ=255)
                                                                                                                     ggAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGAc                          <  2:2731285/148‑1 (MQ=255)
                                                                                                                     ggAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGAc                          >  2:579902/1‑148 (MQ=255)
                                                                                                                            cGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGc                   >  1:3509971/1‑148 (MQ=255)
                                                                                                                                gCGTGAACAGGCACAGACGCGGGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACg               >  1:24436/1‑148 (MQ=255)
                                                                                                                                 cGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGt              <  2:138196/148‑1 (MQ=255)
                                                                                                                                   tGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCTGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGTCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTAc            <  2:1684755/148‑1 (MQ=255)
                                                                                                                                   tGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTAc            <  2:856835/148‑1 (MQ=255)
                                                                                                                                   tGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTAc            >  1:114407/1‑148 (MQ=255)
                                                                                                                                   tGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTAc            <  2:3427592/148‑1 (MQ=255)
                                                                                                                                    gAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTAcc           >  1:543493/1‑148 (MQ=255)
                                                                                                                                     aaCAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCa          >  2:2806491/1‑148 (MQ=255)
                                                                                                                                      acaggcacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTTCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCAt         <  2:635517/148‑1 (MQ=255)
                                                                                                                                      acaggcacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCAt         >  1:2253160/1‑148 (MQ=255)
                                                                                                                                        aggcacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATgg       >  1:1020529/1‑148 (MQ=255)
                                                                                                                                          acacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCa     <  2:1490111/147‑1 (MQ=255)
                                                                                                                                            acagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAat   <  2:2796373/148‑1 (MQ=255)
                                                                                                                                             cagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAata  >  2:830157/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAATA  >  CP009273/397815‑398103

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