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.3% 569.9 / 0.3 174A253T (GCG→ACG) phoAbacterial alkaline phosphatase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (81/90);  total (82/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-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 186 total aligned reads displayed.

TCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAATA  >  CP009273/397814‑398103
                                                                                                                                                 |                                                                                                                                                
tctGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGacg                                                                                                                                                <  2:3425282/148‑1 (MQ=255)
  tGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGAcgcg                                                                                                                                              >  2:50122/1‑148 (MQ=255)
  tGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGAcgcg                                                                                                                                              <  2:1375748/148‑1 (MQ=255)
  tGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGAcgcg                                                                                                                                              <  1:2345801/148‑1 (MQ=255)
   ggAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGt                                                                                                                                             >  1:2942895/1‑148 (MQ=255)
     aaaaaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTgg                                                                                                                                           <  1:1884053/148‑1 (MQ=255)
     aaaaaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTgg                                                                                                                                           <  1:604680/148‑1 (MQ=255)
     aaaaaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTgg                                                                                                                                           >  2:2920247/1‑148 (MQ=255)
        aaGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTa                                                                                                                                        >  1:907111/1‑148 (MQ=255)
         aGGCGGAAACGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTAt                                                                                                                                       <  1:333815/148‑1 (MQ=255)
            cGGAAACGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAg                                                                                                                                    <  1:652955/148‑1 (MQ=255)
            cGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAg                                                                                                                                    <  1:2742948/148‑1 (MQ=255)
            cGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAg                                                                                                                                    <  1:945216/148‑1 (MQ=255)
             ggAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGt                                                                                                                                   <  2:1662404/148‑1 (MQ=255)
                aaaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTgg                                                                                                                                >  1:3211228/1‑148 (MQ=255)
                aaaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTgg                                                                                                                                >  1:453558/1‑148 (MQ=255)
                 aaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGt                                                                                                                               <  2:860085/148‑1 (MQ=255)
                 aaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGt                                                                                                                               >  2:1560348/1‑148 (MQ=255)
                 aaGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGt                                                                                                                               <  1:1281668/148‑1 (MQ=255)
                  aGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTg                                                                                                                              <  1:3457351/148‑1 (MQ=255)
                   ggATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGa                                                                                                                             >  2:847019/1‑148 (MQ=255)
                    gATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAg                                                                                                                            >  2:1369283/1‑148 (MQ=255)
                     aTCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGc                                                                                                                           >  2:657766/1‑148 (MQ=255)
                     aTCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGc                                                                                                                           >  1:3101075/1‑148 (MQ=255)
                     aTCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGc                                                                                                                           >  1:1400974/1‑148 (MQ=255)
                         atgaCCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtg                                                                                                                       <  1:1618505/145‑1 (MQ=255)
                         aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGTCGCAAAAACCTGTGCTGAAACGTCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtg                                                                                                                       <  2:2131852/148‑1 (MQ=255)
                         aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtg                                                                                                                       >  2:3065020/1‑148 (MQ=255)
                         aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtg                                                                                                                       >  2:3065250/1‑148 (MQ=255)
                         aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtg                                                                                                                       <  1:3476161/148‑1 (MQ=255)
                         aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtg                                                                                                                       <  1:3232292/148‑1 (MQ=255)
                         aTTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtg                                                                                                                       >  2:2435148/1‑148 (MQ=255)
                          ttACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtgc                                                                                                                      <  1:2506463/148‑1 (MQ=255)
                          ttACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtgc                                                                                                                      <  1:717633/148‑1 (MQ=255)
                            aCCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtgctg                                                                                                                    >  2:2940253/1‑148 (MQ=255)
                             ccGAAGAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCCTTTCTGATACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATAAGTTGGTGAGCGAtgctgc                                                                                                                   <  1:2034766/148‑1 (MQ=255)
                             ccGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGAtgctgc                                                                                                                   <  1:2226200/148‑1 (MQ=255)
                                aaCAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTc                                                                                                                >  2:776559/1‑148 (MQ=255)
                                  cAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCAc                                                                                                              >  2:132897/1‑148 (MQ=255)
                                     ctgctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGa                                                                                                           >  1:1120092/1‑148 (MQ=255)
                                      tgctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGaa                                                                                                          <  2:1924576/148‑1 (MQ=255)
                                       gctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAAt                                                                                                         >  1:749562/1‑148 (MQ=255)
                                        ctTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAAtt                                                                                                        >  2:3174369/1‑148 (MQ=255)
                                         ttAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTc                                                                                                       <  1:2454402/148‑1 (MQ=255)
                                             cGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCTGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTg                                                                                                   <  1:592441/148‑1 (MQ=255)
                                                 cGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACgg                                                                                               >  1:858991/1‑148 (MQ=255)
                                                      cGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCg                                                                                          >  1:2941845/1‑148 (MQ=255)
                                                      cGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCg                                                                                          >  1:2480154/1‑148 (MQ=255)
                                                       gACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGa                                                                                         >  1:1987927/1‑148 (MQ=255)
                                                        aCGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGaa                                                                                        <  2:855593/148‑1 (MQ=255)
                                                        aCGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGaa                                                                                        <  1:1913550/148‑1 (MQ=255)
                                                         cGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAAt                                                                                       <  2:3328068/148‑1 (MQ=255)
                                                          gTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATc                                                                                      <  1:1683271/148‑1 (MQ=255)
                                                                  tGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCaaaaa                                                                              >  1:695002/1‑148 (MQ=255)
                                                                  tGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCaaaaa                                                                              >  1:2709320/1‑148 (MQ=255)
                                                                   ggcggcggcGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAAc                                                                             <  2:1689675/148‑1 (MQ=255)
                                                                     cggcggcgCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAAccc                                                                           <  1:2303581/148‑1 (MQ=255)
                                                                             cAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTgg                                                                   >  1:413556/1‑148 (MQ=255)
                                                                                 aaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTg                                                               <  2:2103825/148‑1 (MQ=255)
                                                                                 aaCCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTg                                                               >  2:2987308/1‑148 (MQ=255)
                                                                                    cTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGttt                                                            >  1:1566684/1‑148 (MQ=255)
                                                                                      ttGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGc                                                          <  1:998802/148‑1 (MQ=255)
                                                                                      ttGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGc                                                          >  1:163502/1‑148 (MQ=255)
                                                                                             aaCGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGc                                                   >  2:2759313/1‑148 (MQ=255)
                                                                                              aCGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCa                                                  >  1:2805673/1‑148 (MQ=255)
                                                                                              aCGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCa                                                  >  1:2804863/1‑148 (MQ=255)
                                                                                                  cAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAatat                                              >  2:1448081/1‑148 (MQ=255)
                                                                                                   aaCCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATg                                             <  1:1411282/148‑1 (MQ=255)
                                                                                                     ccGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGcc                                           >  1:3046977/1‑148 (MQ=255)
                                                                                                      cGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCTATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCa                                          >  1:1822478/1‑148 (MQ=255)
                                                                                                      cGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCa                                          <  1:117845/148‑1 (MQ=255)
                                                                                                        cTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGt                                        <  2:2697917/148‑1 (MQ=255)
                                                                                                           gTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTgcg                                     >  2:9728/1‑148 (MQ=255)
                                                                                                            tGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTgcgc                                    <  2:722194/148‑1 (MQ=255)
                                                                                                            tGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTgcgc                                    <  2:722210/148‑1 (MQ=255)
                                                                                                            tGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTgcgc                                    >  2:1131600/1‑148 (MQ=255)
                                                                                                            tGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTgcgc                                    <  2:2509732/148‑1 (MQ=255)
                                                                                                              aaTGGCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTg                                  >  2:2062142/1‑148 (MQ=255)
                                                                                                                 ggCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCt                               <  2:2241425/148‑1 (MQ=255)
                                                                                                                 ggCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCt                               <  2:1806196/148‑1 (MQ=255)
                                                                                                                  gCAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTa                              >  2:362660/1‑148 (MQ=255)
                                                                                                                   cAGGGAAAAACGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAg                             >  1:860408/1‑148 (MQ=255)
                                                                                                                           aaCGCTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGaaa                     <  1:3511013/148‑1 (MQ=255)
                                                                                                                               cTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCaa                 <  1:1445816/148‑1 (MQ=255)
                                                                                                                               cTGCGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCaa                 <  1:50122/148‑1 (MQ=255)
                                                                                                                                  cGTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGt              >  2:1240505/1‑148 (MQ=255)
                                                                                                                                   gTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTa             <  2:3233842/148‑1 (MQ=255)
                                                                                                                                   gTGAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTa             <  2:269343/148‑1 (MQ=255)
                                                                                                                                     gAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTAcc           >  1:1973948/1‑148 (MQ=255)
                                                                                                                                     gAACAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTAcc           >  1:2683108/1‑148 (MQ=255)
                                                                                                                                      aaCAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCa          >  2:3316573/1‑148 (MQ=255)
                                                                                                                                      aaCAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCa          >  2:2774986/1‑148 (MQ=255)
                                                                                                                                      aaCAGGCACAGACGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATACGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCa          >  2:993867/1‑148 (MQ=255)
                                                                                                                                         aggcacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATgg       >  2:1748554/1‑148 (MQ=255)
                                                                                                                                         aggcacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATgg       >  2:1704905/1‑148 (MQ=255)
                                                                                                                                          ggcacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGc      >  2:3516944/1‑148 (MQ=255)
                                                                                                                                           gcacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCa     >  1:897387/1‑148 (MQ=255)
                                                                                                                                           gcacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCa     >  1:2594/1‑148 (MQ=255)
                                                                                                                                            cacagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCaa    >  1:390007/1‑148 (MQ=255)
                                                                                                                                              cagacGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAata  <  2:1800514/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCTGGAAAAAGGCGGAAAAGGATCGATTACCGAACAGCTGCTTAACGCTCGTGCCGACGTTACGCTTGGCGGCGGCGCAAAAACCTTTGCTGAAACGGCAACCGCTGGTGAATGGCAGGGAAAAACGCTGCGTGAACAGGCACAGGCGCGTGGTTATCAGTTGGTGAGCGATGCTGCCTCACTGAATTCGGTGACGGAAGCGAATCAGCAAAAACCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAATA  >  CP009273/397814‑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: ‑