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,9590GA97.7% 846.9 / 6.5 263A253T (GCG→ACG) phoAbacterial alkaline phosphatase
Reads supporting (aligned to +/- strand):  ref base G (6/0);  new base A (151/104);  total (158/105)
Fisher's exact test for biased strand distribution p-value = 8.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-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 281 total aligned reads displayed.

GCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAATAT  >  CP009273/397812‑398104
                                                                                                                                                   |                                                                                                                                                 
 ctctGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTAAACGCTCGTGCCGACGTTACTCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGacgc                                                                                                                                                >  1:1181896/1‑150 (MQ=255)
  tctGGATAAAGGCGGAAAAGGATCGATTACCGATCTGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAATCGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGAcgcg                                                                                                                                               <  2:1901329/150‑1 (MQ=255)
   ctGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGt                                                                                                                                              <  2:4676631/150‑1 (MQ=255)
   ctGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGt                                                                                                                                              >  2:1772325/1‑150 (MQ=255)
     ggAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCGGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTgg                                                                                                                                            >  2:2349147/1‑150 (MQ=255)
       aaaaaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGtt                                                                                                                                          >  2:133144/1‑150 (MQ=255)
         aaaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTAt                                                                                                                                        >  2:4255421/1‑150 (MQ=255)
         aaaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTAt                                                                                                                                        >  2:4195138/1‑150 (MQ=255)
          aaGGCGGAAAAGGATCGATTACCGACCAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATc                                                                                                                                       <  1:132363/150‑1 (MQ=255)
          aaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGACACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATc                                                                                                                                       >  2:4867060/1‑150 (MQ=255)
          aaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATc                                                                                                                                       >  1:4937390/1‑150 (MQ=255)
          aaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTATGCTGAAACGGCAAACGCTGGTGAAATGCAGGGAAAAACGCTGCGTGAACATGCACATACTCGTGGTAATc                                                                                                                                       >  1:4870895/1‑150 (MQ=255)
           aGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCa                                                                                                                                      >  1:1068654/1‑150 (MQ=255)
              cGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCTTGCCGACGTTACGCTTGGCGGCGGAGCAAAAACCTTTTCTGAAACGGCAACCGCTGTTGAATGGCAGGGAAAAACGCTGCTTGAACAGGCACAGACGCTTGGTTATCAGtt                                                                                                                                   >  1:1958109/1‑150 (MQ=255)
              cGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGtt                                                                                                                                   <  1:732704/150‑1 (MQ=255)
               ggAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTg                                                                                                                                  >  1:2354805/1‑150 (MQ=255)
               ggAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTg                                                                                                                                  >  2:5085440/1‑150 (MQ=255)
                 aaaaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGt                                                                                                                                <  2:517328/150‑1 (MQ=255)
                 aaaaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGt                                                                                                                                >  1:1237351/1‑150 (MQ=255)
                  aaaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTg                                                                                                                               >  1:3015861/1‑150 (MQ=255)
                     ggATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGc                                                                                                                            >  1:101288/1‑150 (MQ=255)
                        tCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAt                                                                                                                         >  2:3931373/1‑150 (MQ=255)
                          gATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGAtgc                                                                                                                       >  2:3277265/1‑150 (MQ=255)
                           aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGTGCGAtgct                                                                                                                      <  2:5303702/150‑1 (MQ=255)
                           aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtgct                                                                                                                      >  2:3816961/1‑150 (MQ=255)
                           aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtgct                                                                                                                      <  1:3872501/150‑1 (MQ=255)
                           aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtgct                                                                                                                      <  1:1822962/150‑1 (MQ=255)
                               ccGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCt                                                                                                                  >  1:1459851/1‑150 (MQ=255)
                                cGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTc                                                                                                                 >  1:815596/1‑150 (MQ=255)
                                 gAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCa                                                                                                                >  2:4653780/1‑150 (MQ=255)
                                  aaCAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCATTGATGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCAc                                                                                                               >  1:3796948/1‑150 (MQ=255)
                                   aCAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACt                                                                                                              <  1:153201/150‑1 (MQ=255)
                                     aGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGa                                                                                                            <  2:2030259/150‑1 (MQ=255)
                                       ctgctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAAt                                                                                                          >  2:2637251/1‑150 (MQ=255)
                                        tgctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAAtt                                                                                                         <  1:2248270/150‑1 (MQ=255)
                                        tgctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAAtt                                                                                                         >  2:5256815/1‑150 (MQ=255)
                                         gtttAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTc                                                                                                        <  1:1917482/148‑1 (MQ=255)
                                         gctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTc                                                                                                        <  2:3652787/150‑1 (MQ=255)
                                            tAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGt                                                                                                     >  1:32081/1‑150 (MQ=255)
                                             aaCGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTg                                                                                                    <  2:2532943/150‑1 (MQ=255)
                                              aCGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGa                                                                                                   <  2:1073408/150‑1 (MQ=255)
                                                gCTCGTGCCGACGTTACGCTTGGCGGTGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGTTGCCTCACTGAATTCGGTGACg                                                                                                 <  2:241461/150‑1 (MQ=255)
                                                    gTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAg                                                                                             <  1:2093182/150‑1 (MQ=255)
                                                        cGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGaa                                                                                         <  2:893558/150‑1 (MQ=255)
                                                        cGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGaa                                                                                         >  1:2453696/1‑150 (MQ=255)
                                                          aCGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATc                                                                                       <  2:4623621/150‑1 (MQ=255)
                                                            gTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAg                                                                                     <  1:4389042/150‑1 (MQ=255)
                                                            gTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAg                                                                                     >  2:1270207/1‑150 (MQ=255)
                                                              tACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCa                                                                                   >  1:2045014/1‑150 (MQ=255)
                                                              tACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCa                                                                                   <  2:4339719/150‑1 (MQ=255)
                                                                cGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACTGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCaaa                                                                                 <  2:37184/150‑1 (MQ=255)
                                                                cGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCaaa                                                                                 <  1:4572388/150‑1 (MQ=255)
                                                                    tGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAAcc                                                                             >  1:167130/1‑150 (MQ=255)
                                                                    tGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAAcc                                                                             >  1:5328023/1‑150 (MQ=255)
                                                                        ggcggcgCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTg                                                                         <  1:3618590/150‑1 (MQ=255)
                                                                        ggcggcgCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTg                                                                         >  2:169065/1‑150 (MQ=255)
                                                                         gcggcgCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGACGCGAATCAGCAAAAACCCCTGc                                                                        >  2:1805783/1‑150 (MQ=255)
                                                                          cggcgCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACAGAATTCGGTGACGGAAGAGAATCAGCAAAAACCCCTGCt                                                                       >  2:4037336/1‑150 (MQ=255)
                                                                            gcgcAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTg                                                                     >  1:421980/1‑150 (MQ=255)
                                                                            gcgcAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTg                                                                     >  2:5203470/1‑150 (MQ=255)
                                                                            gcgcAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTg                                                                     <  1:402594/150‑1 (MQ=255)
                                                                                 aaaaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTg                                                                >  2:556862/1‑150 (MQ=255)
                                                                                 aaaaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTg                                                                >  2:2376275/1‑150 (MQ=255)
                                                                                  aaaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGt                                                               >  2:3157688/1‑150 (MQ=255)
                                                                                   aaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAACACCCCTGCTTGGCCTGtt                                                              <  2:3451309/150‑1 (MQ=255)
                                                                                   aaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGtt                                                              >  1:1418015/1‑150 (MQ=255)
                                                                                    aCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGttt                                                             <  1:5021397/150‑1 (MQ=255)
                                                                                     ccTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTg                                                            >  2:3538857/1‑150 (MQ=255)
                                                                                     ccTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTg                                                            <  1:3305162/150‑1 (MQ=255)
                                                                                      cTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGc                                                           >  2:1945027/1‑150 (MQ=255)
                                                                                      cTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGc                                                           >  2:2877450/1‑150 (MQ=255)
                                                                                         tGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGa                                                        <  1:3129226/150‑1 (MQ=255)
                                                                                            tGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACgg                                                     <  2:4506092/150‑1 (MQ=255)
                                                                                               aaCGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCaa                                                  >  2:2746389/1‑150 (MQ=255)
                                                                                               aaCGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCaa                                                  <  2:2542736/150‑1 (MQ=255)
                                                                                                aCGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAat                                                 >  2:4892238/1‑150 (MQ=255)
                                                                                                    cAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGc                                             >  1:2383367/1‑150 (MQ=255)
                                                                                                         gCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTg                                        <  1:3912697/150‑1 (MQ=255)
                                                                                                             gTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCt                                    <  2:619225/150‑1 (MQ=255)
                                                                                                               gAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTgg                                  >  1:3669724/1‑150 (MQ=255)
                                                                                                                aaTGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGc                                 >  2:1218016/1‑150 (MQ=255)
                                                                                                                      aGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGAc                           >  1:3316184/1‑150 (MQ=255)
                                                                                                                       gggAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGAcc                          <  1:349559/150‑1 (MQ=255)
                                                                                                                         gAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGa                        <  2:2924667/150‑1 (MQ=255)
                                                                                                                          aaaaaCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGaa                       >  2:4422619/1‑150 (MQ=255)
                                                                                                                           aaaaCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGaaa                      >  1:5135633/1‑150 (MQ=255)
                                                                                                                           aaaaCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGaaa                      >  2:2682529/1‑150 (MQ=255)
                                                                                                                            aaaCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAg                     <  2:2129320/150‑1 (MQ=255)
                                                                                                                              aCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAGTATGCCAGTGCGCTGGCTAGGACCGAAAGCa                   <  2:4998414/150‑1 (MQ=255)
                                                                                                                               cGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCaa                  <  2:3421252/150‑1 (MQ=255)
                                                                                                                                 cTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACg                <  2:4937361/150‑1 (MQ=255)
                                                                                                                                 cTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACg                >  1:2635084/1‑150 (MQ=255)
                                                                                                                                 cTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACg                <  2:4942845/150‑1 (MQ=255)
                                                                                                                                 cTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACg                >  1:933926/1‑150 (MQ=255)
                                                                                                                                     gTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCATGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTAcc            >  2:2391373/1‑150 (MQ=255)
                                                                                                                                      tGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCa           <  1:3439930/150‑1 (MQ=255)
                                                                                                                                        aaCAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATg         <  2:3034587/150‑1 (MQ=255)
                                                                                                                                           aggcacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCa      <  2:281631/150‑1 (MQ=255)
                                                                                                                                               acagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAatat  <  1:278899/150‑1 (MQ=255)
                                                                                                                                               acagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAatat  <  1:2936515/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCTCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAATAT  >  CP009273/397812‑398104

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