Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,885,611 G→T 100% R628S (CGT→AGT)  yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,885,6110GT97.3% 615.5 / 6.4 203R628S (CGT→AGT) yoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base T (121/76);  total (124/79)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 203 total aligned reads displayed.

TCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACC  >  CP009273/1885466‑1885753
                                                                                                                                                 |                                                                                                                                              
tcgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTAGCTGGAGGATGGTATCGCAAGGAAACTaa                                                                                                                                              >  1:1529064/1‑148 (MQ=255)
tcgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTaa                                                                                                                                              >  1:5694/1‑148 (MQ=255)
tcgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTaa                                                                                                                                              >  1:2286175/1‑148 (MQ=255)
tcgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTaa                                                                                                                                              >  1:1530397/1‑148 (MQ=255)
   ttCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCg                                                                                                                                           <  2:1672496/148‑1 (MQ=255)
             cGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAAt                                                                                                                                 >  1:3480120/1‑148 (MQ=255)
             cGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAAt                                                                                                                                 >  1:1717055/1‑148 (MQ=255)
                cAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTc                                                                                                                              <  1:3509082/148‑1 (MQ=255)
                 agagCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCa                                                                                                                             <  1:3170362/148‑1 (MQ=255)
                 agagCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCa                                                                                                                             >  2:45356/1‑148 (MQ=255)
                 agagCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCa                                                                                                                             <  1:1396700/148‑1 (MQ=255)
                       ggCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTg                                                                                                                       >  2:806350/1‑148 (MQ=255)
                        gCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGc                                                                                                                      >  2:902998/1‑148 (MQ=255)
                        gCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGc                                                                                                                      >  1:3270501/1‑148 (MQ=255)
                         ccTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCg                                                                                                                     <  2:2775579/148‑1 (MQ=255)
                          ctctGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGt                                                                                                                    <  1:2592856/148‑1 (MQ=255)
                          ctctGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGt                                                                                                                    >  1:3244838/1‑148 (MQ=255)
                           tctGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTg                                                                                                                   <  2:1148768/148‑1 (MQ=255)
                           tctGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTg                                                                                                                   <  1:1850379/148‑1 (MQ=255)
                           tctGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTg                                                                                                                   <  2:1151685/148‑1 (MQ=255)
                                  cgGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCgcggg                                                                                                            >  1:1320554/1‑148 (MQ=255)
                                     tATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGACACTAACCGCACGGGCAATGTCACGGGTGCGTGGCgcgggcgg                                                                                                         >  1:3378972/1‑148 (MQ=255)
                                     tATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCgcgggcgg                                                                                                         >  2:1239451/1‑148 (MQ=255)
                                     tATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCgcgggcgg                                                                                                         <  2:2502764/148‑1 (MQ=255)
                                     tATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCgcgggcgg                                                                                                         >  2:2585505/1‑148 (MQ=255)
                                      aTCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGggcggc                                                                                                        <  2:3212877/148‑1 (MQ=255)
                                        cGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAg                                                                                                      >  2:1255951/1‑148 (MQ=255)
                                        cGAGAGCCAGAAGTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCAGACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAg                                                                                                      <  2:2407960/148‑1 (MQ=255)
                                         gATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGa                                                                                                     <  1:3285632/148‑1 (MQ=255)
                                         gATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGa                                                                                                     >  1:641058/1‑148 (MQ=255)
                                          aTAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGAc                                                                                                    >  2:3081920/1‑148 (MQ=255)
                                           tAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACt                                                                                                   <  2:651347/148‑1 (MQ=255)
                                           tAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACt                                                                                                   <  2:2059224/148‑1 (MQ=255)
                                            aGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTg                                                                                                  >  2:1290634/1‑148 (MQ=255)
                                                  aaTTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCgaga                                                                                            >  1:2898545/1‑148 (MQ=255)
                                                  aaTTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCgaga                                                                                            >  2:2011940/1‑148 (MQ=255)
                                                   aTTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGaa                                                                                           >  1:195549/1‑148 (MQ=255)
                                                     tCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAAc                                                                                         <  2:3119771/148‑1 (MQ=255)
                                                      cGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACg                                                                                        <  1:1200398/148‑1 (MQ=255)
                                                          ttggTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTcgc                                                                                    >  1:2987206/3‑148 (MQ=255)
                                                          tGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTcgc                                                                                    >  2:2708506/1‑148 (MQ=255)
                                                           gggTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTcgcg                                                                                   >  1:1634935/1‑148 (MQ=255)
                                                            ggTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCgcgc                                                                                  >  2:946317/1‑148 (MQ=255)
                                                            ggTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCgcgc                                                                                  >  1:2551490/1‑148 (MQ=255)
                                                            ggTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCgcgc                                                                                  >  1:2550845/1‑148 (MQ=255)
                                                            ggTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCgcgc                                                                                  <  1:2084673/148‑1 (MQ=255)
                                                                    tGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGa                                                                          <  2:1266762/148‑1 (MQ=255)
                                                                        tAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCa                                                                      >  1:282116/1‑148 (MQ=255)
                                                                        tAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCa                                                                      >  1:2131438/1‑148 (MQ=255)
                                                                          gATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATc                                                                    >  2:1547865/1‑148 (MQ=255)
                                                                          gATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATAACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGAGCGGGCGTCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATc                                                                    >  1:142566/1‑148 (MQ=255)
                                                                            tCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACNGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCAc                                                                  >  1:1807028/1‑148 (MQ=255)
                                                                            tCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCAc                                                                  <  1:2445733/148‑1 (MQ=255)
                                                                              aaTAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCa                                                                >  1:3056143/1‑148 (MQ=255)
                                                                               aTAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCAGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACAGGCGAGAAACGTCGAGCCGTAAGGACGCATCACCAg                                                               >  2:1483997/1‑148 (MQ=255)
                                                                                tAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGc                                                              >  1:2145467/1‑148 (MQ=255)
                                                                                  aaaaGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCg                                                            >  1:2397290/1‑148 (MQ=255)
                                                                                   aaaGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGa                                                           >  2:3534031/1‑148 (MQ=255)
                                                                                       gcgcgcATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGt                                                       >  2:3425583/1‑148 (MQ=255)
                                                                                        cgcgcATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTc                                                      <  1:426382/148‑1 (MQ=255)
                                                                                              tcatACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAAt                                                >  1:2428363/1‑148 (MQ=255)
                                                                                              tcatACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGAGTGTCACAAAt                                                >  2:2262217/1‑148 (MQ=255)
                                                                                               cataccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATc                                               >  1:1772294/1‑148 (MQ=255)
                                                                                                ataccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCa                                              <  2:2916734/148‑1 (MQ=255)
                                                                                                ataccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACAAGCCGATTGTCACAAATCa                                              >  1:2909494/1‑148 (MQ=255)
                                                                                                  accatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCAcc                                            >  2:1314318/1‑148 (MQ=255)
                                                                                                     atacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGc                                         >  2:3147766/1‑148 (MQ=255)
                                                                                                      tacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCa                                        <  2:2446453/148‑1 (MQ=255)
                                                                                                      tacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCa                                        >  1:2366290/1‑148 (MQ=255)
                                                                                                        cTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACg                                      <  1:2432721/148‑1 (MQ=255)
                                                                                                          ccGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGcc                                    <  1:492495/148‑1 (MQ=255)
                                                                                                           cGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCa                                   <  1:2122277/148‑1 (MQ=255)
                                                                                                              aaCAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGa                                >  2:1125578/1‑148 (MQ=255)
                                                                                                              aaCAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGa                                <  1:899523/148‑1 (MQ=255)
                                                                                                               aCAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAt                               <  2:197340/148‑1 (MQ=255)
                                                                                                               aCAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAt                               >  2:666186/1‑148 (MQ=255)
                                                                                                                 aaaTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcg                             <  1:1226024/148‑1 (MQ=255)
                                                                                                                 aaaTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcg                             >  1:792546/1‑148 (MQ=255)
                                                                                                                  aaTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgt                            <  1:3110665/148‑1 (MQ=255)
                                                                                                                   aTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgtc                           >  2:2133929/1‑148 (MQ=255)
                                                                                                                   aTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgtc                           >  1:532138/1‑148 (MQ=255)
                                                                                                                   aTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgtc                           >  1:532115/1‑148 (MQ=255)
                                                                                                                     tACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCgg                         >  1:3000499/1‑148 (MQ=255)
                                                                                                                     tACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCgg                         >  2:2651865/1‑148 (MQ=255)
                                                                                                                          ggaggaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTcg                    <  1:2671433/148‑1 (MQ=255)
                                                                                                                             ggaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGa                 >  2:698775/1‑148 (MQ=255)
                                                                                                                              gaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGaa                >  1:2427709/1‑148 (MQ=255)
                                                                                                                              gaTGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGaa                <  1:1674763/148‑1 (MQ=255)
                                                                                                                               aTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAAt               <  1:1134237/148‑1 (MQ=255)
                                                                                                                                tGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATc              <  1:2713139/148‑1 (MQ=255)
                                                                                                                                tGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATc              >  1:1268518/1‑148 (MQ=255)
                                                                                                                                  gTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAg            <  1:58188/148‑1 (MQ=255)
                                                                                                                                   tATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGt           <  1:3160770/148‑1 (MQ=255)
                                                                                                                                    aTCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTc          >  1:2352535/1‑148 (MQ=255)
                                                                                                                                       gCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGcc       <  1:786318/148‑1 (MQ=255)
                                                                                                                                       gCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGcc       <  1:902601/148‑1 (MQ=255)
                                                                                                                                        cAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGccc      <  2:2987206/148‑1 (MQ=255)
                                                                                                                                        cAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGccc      <  1:596885/148‑1 (MQ=255)
                                                                                                                                          aGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTa    >  1:1010344/1‑148 (MQ=255)
                                                                                                                                            gAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTAcc  >  2:1169324/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACC  >  CP009273/1885466‑1885753

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