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,9590GA96.7% 683.1 / 10.5 210A253T (GCG→ACG) phoAbacterial alkaline phosphatase
Reads supporting (aligned to +/- strand):  ref base G (1/6);  new base A (103/100);  total (104/106)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
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 229 total aligned reads displayed.

TGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAATA  >  CP009273/397816‑398103
                                                                                                                                               |                                                                                                                                                
tGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGAcgcg                                                                                                                                              <  2:3732035/148‑1 (MQ=255)
tGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGATTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGAcgcg                                                                                                                                              <  1:2888613/148‑1 (MQ=255)
 ggAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGt                                                                                                                                             >  1:829154/1‑148 (MQ=255)
 ggAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGt                                                                                                                                             >  1:829631/1‑148 (MQ=255)
  gAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTg                                                                                                                                            >  2:542486/1‑148 (MQ=255)
      aaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTa                                                                                                                                        >  2:1706400/1‑148 (MQ=255)
         gCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCa                                                                                                                                     >  1:2602810/1‑148 (MQ=255)
         gCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCa                                                                                                                                     >  1:2601698/1‑148 (MQ=255)
          cGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAg                                                                                                                                    <  2:2020462/148‑1 (MQ=255)
           ggAAACGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGt                                                                                                                                   <  2:1100212/148‑1 (MQ=255)
           ggAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGt                                                                                                                                   <  1:1248891/148‑1 (MQ=255)
            gAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGtt                                                                                                                                  <  1:1486264/148‑1 (MQ=255)
            gAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGtt                                                                                                                                  <  1:3387799/148‑1 (MQ=255)
             aaaaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTg                                                                                                                                 <  1:2059416/148‑1 (MQ=255)
               aaGGATCGATTACAGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGt                                                                                                                               >  1:2769611/1‑148 (MQ=255)
                 ggATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGa                                                                                                                             >  1:1095594/1‑148 (MQ=255)
                    tCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCg                                                                                                                          <  2:424094/148‑1 (MQ=255)
                    tCGAGTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCg                                                                                                                          <  2:2346753/148‑1 (MQ=255)
                       aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtg                                                                                                                       >  1:1571821/1‑148 (MQ=255)
                       aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtg                                                                                                                       >  1:1406328/1‑148 (MQ=255)
                       aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtg                                                                                                                       >  2:1267339/1‑148 (MQ=255)
                        ttACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtgc                                                                                                                      <  2:1779363/148‑1 (MQ=255)
                        ttACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtgc                                                                                                                      >  1:1955061/1‑148 (MQ=255)
                          aCCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtgctg                                                                                                                    >  1:975361/1‑148 (MQ=255)
                              aaCAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTc                                                                                                                >  2:2244374/1‑148 (MQ=255)
                               aCAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCa                                                                                                               >  1:3280312/1‑148 (MQ=255)
                                  gctgctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTg                                                                                                            >  1:3545963/1‑148 (MQ=255)
                                   ctgctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGa                                                                                                           >  2:3669888/1‑148 (MQ=255)
                                     gctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAAt                                                                                                         >  1:2650821/1‑148 (MQ=255)
                                     gctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAAt                                                                                                         <  1:315938/148‑1 (MQ=255)
                                     gctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAAt                                                                                                         >  1:900553/1‑148 (MQ=255)
                                      ctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAAtt                                                                                                        <  1:1689052/148‑1 (MQ=255)
                                      ctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAAtt                                                                                                        <  2:3132643/148‑1 (MQ=255)
                                      ctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAAtt                                                                                                        >  1:3620240/1‑148 (MQ=255)
                                        tAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCg                                                                                                      <  2:1476847/148‑1 (MQ=255)
                                        tAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCg                                                                                                      <  1:2297119/148‑1 (MQ=255)
                                        tAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCg                                                                                                      <  1:373399/148‑1 (MQ=255)
                                          aCGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGt                                                                                                    <  2:3059278/148‑1 (MQ=255)
                                             cTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGAc                                                                                                 >  2:31707/1‑148 (MQ=255)
                                              tCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACg                                                                                                >  1:1132094/1‑148 (MQ=255)
                                                  gCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTTACGGAAg                                                                                            <  2:3079887/148‑1 (MQ=255)
                                                       cGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAAt                                                                                       <  1:2760948/148‑1 (MQ=255)
                                                         ttACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGATTCa                                                                                     <  2:1674772/148‑1 (MQ=255)
                                                         ttACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCa                                                                                     >  1:172525/1‑148 (MQ=255)
                                                           aCGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGc                                                                                   >  2:2165030/1‑148 (MQ=255)
                                                              cTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCaaa                                                                                >  2:2400152/1‑148 (MQ=255)
                                                              cTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCaaa                                                                                >  2:1570765/1‑148 (MQ=255)
                                                              cTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCaaa                                                                                >  2:2890611/1‑148 (MQ=255)
                                                                   cggcggcgCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAAccc                                                                           <  2:2754618/148‑1 (MQ=255)
                                                                             aaaaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGcc                                                                 >  1:3497524/1‑148 (MQ=255)
                                                                             aaaaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGcc                                                                 <  2:1406328/148‑1 (MQ=255)
                                                                             aaaaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGcc                                                                 >  2:3013610/1‑148 (MQ=255)
                                                                             aaaaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGcc                                                                 >  2:2718710/1‑148 (MQ=255)
                                                                              aaaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCt                                                                >  1:1032103/1‑148 (MQ=255)
                                                                              aaaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCt                                                                >  2:3637044/1‑148 (MQ=255)
                                                                              aaaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCt                                                                <  2:1361446/148‑1 (MQ=255)
                                                                              aaaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCt                                                                <  2:2650821/148‑1 (MQ=255)
                                                                               aaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTg                                                               <  2:1576266/148‑1 (MQ=255)
                                                                               aaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGC‑TGGCCTGt                                                              <  2:643230/148‑1 (MQ=255)
                                                                                aCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGt                                                              <  2:358972/148‑1 (MQ=255)
                                                                                aCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGt                                                              <  2:465707/148‑1 (MQ=255)
                                                                                 ccTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGtt                                                             >  2:2852152/1‑148 (MQ=255)
                                                                                    ttGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGc                                                          <  1:2718504/148‑1 (MQ=255)
                                                                                    ttGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGc                                                          >  2:1825920/1‑148 (MQ=255)
                                                                                    ttGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGc                                                          >  2:1792556/1‑148 (MQ=255)
                                                                                             cGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCaa                                                 <  1:3324816/148‑1 (MQ=255)
                                                                                              ggCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAat                                                <  1:1179557/148‑1 (MQ=255)
                                                                                               gCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAata                                               <  2:2520517/148‑1 (MQ=255)
                                                                                                cAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAatat                                              >  2:3280241/1‑148 (MQ=255)
                                                                                                   ccGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGcc                                           <  2:1954122/148‑1 (MQ=255)
                                                                                                   ccGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGcc                                           >  2:2756423/1‑148 (MQ=255)
                                                                                                   ccGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGcc                                           <  1:3120693/148‑1 (MQ=255)
                                                                                                       tGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTg                                       <  1:748861/148‑1 (MQ=255)
                                                                                                       tGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTg                                       <  1:619905/148‑1 (MQ=255)
                                                                                                          tGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTgcgc                                    <  1:1307902/148‑1 (MQ=255)
                                                                                                             aTGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTgg                                 >  2:1893341/1‑148 (MQ=255)
                                                                                                             aTGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTgg                                 <  2:3462406/148‑1 (MQ=255)
                                                                                                                gCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGATa                              >  2:2929572/1‑148 (MQ=255)
                                                                                                                    ggAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGAc                          >  1:917710/1‑148 (MQ=255)
                                                                                                                       aaaaCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGa                       >  1:1333521/1‑148 (MQ=255)
                                                                                                                        aaaCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGaa                      >  1:594180/1‑148 (MQ=255)
                                                                                                                        aaaCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGaa                      <  1:184221/148‑1 (MQ=255)
                                                                                                                        aaaCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGaa                      <  1:2093017/148‑1 (MQ=255)
                                                                                                                        aaaCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGaa                      <  1:1976514/148‑1 (MQ=255)
                                                                                                                         aaCGATGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGaaa                     <  2:2904901/148‑1 (MQ=255)
                                                                                                                          aCGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAg                    <  1:2405730/148‑1 (MQ=255)
                                                                                                                          aCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAg                    >  1:805998/1‑148 (MQ=255)
                                                                                                                          aCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAg                    >  1:1384165/1‑148 (MQ=255)
                                                                                                                           cGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGc                   <  2:3128490/148‑1 (MQ=255)
                                                                                                                           cGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGc                   >  1:2073910/1‑148 (MQ=255)
                                                                                                                            gCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCa                  <  2:707015/148‑1 (MQ=255)
                                                                                                                                cGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGt              >  1:3217567/1‑148 (MQ=255)
                                                                                                                                 gTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTa             <  1:2165030/148‑1 (MQ=255)
                                                                                                                                 gTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTa             <  1:450238/148‑1 (MQ=255)
                                                                                                                                     acaggcacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCAt         >  2:898141/1‑148 (MQ=255)
                                                                                                                                      caggcacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATg        >  2:314438/1‑148 (MQ=255)
                                                                                                                                         gcacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCa     >  1:2556839/1‑148 (MQ=255)
                                                                                                                                          cacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCACCGTACCATGGCaa    >  2:416645/1‑148 (MQ=255)
                                                                                                                                           acagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAat   <  1:2270664/148‑1 (MQ=255)
                                                                                                                                            cagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAata  <  1:1458609/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAATA  >  CP009273/397816‑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: ‑