Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 726,087 G→T 100% R350L (CGT→CTT)  rhsC → Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273726,0870GT99.4% 502.9 / ‑4.0 156R350L (CGT→CTT) rhsCRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (79/76);  total (79/77)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-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 156 total aligned reads displayed.

TACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAA  >  CP009273/725942‑726229
                                                                                                                                                 |                                                                                                                                              
tACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTc                                                                                                                                              <  1:2550859/148‑1 (MQ=255)
   ccGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTcaca                                                                                                                                           <  1:2371659/148‑1 (MQ=255)
      gagaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACgg                                                                                                                                        <  1:295988/148‑1 (MQ=255)
      gagaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACgg                                                                                                                                        <  1:295983/148‑1 (MQ=255)
      gagaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACgg                                                                                                                                        >  1:669490/1‑148 (MQ=255)
          aTTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCg                                                                                                                                    >  1:1758078/1‑148 (MQ=255)
                 tGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGaa                                                                                                                             >  1:3610680/1‑148 (MQ=255)
                 tGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGaa                                                                                                                             <  1:2310148/148‑1 (MQ=255)
                         cgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTa                                                                                                                     <  2:1342975/148‑1 (MQ=255)
                         cgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTa                                                                                                                     <  1:1423816/148‑1 (MQ=255)
                         cgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTa                                                                                                                     <  2:1342972/148‑1 (MQ=255)
                            tGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCg                                                                                                                  >  2:2031130/1‑148 (MQ=255)
                            tGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCg                                                                                                                  >  2:2030814/1‑148 (MQ=255)
                            tGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCg                                                                                                                  <  2:3013568/148‑1 (MQ=255)
                              gTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGtt                                                                                                                >  1:3049717/1‑148 (MQ=255)
                              gTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGtt                                                                                                                <  2:2846714/148‑1 (MQ=255)
                              gTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGtt                                                                                                                >  1:3049254/1‑148 (MQ=255)
                                    tATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACa                                                                                                          <  1:1215306/148‑1 (MQ=255)
                                    tATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACa                                                                                                          <  1:1215198/148‑1 (MQ=255)
                                      tGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGc                                                                                                        >  1:707259/1‑148 (MQ=255)
                                        gCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGa                                                                                                      >  2:520654/1‑148 (MQ=255)
                                            gACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGACATCTGTTACCGTTACGACAGCGAcggg                                                                                                  >  2:3227208/1‑148 (MQ=255)
                                            gACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGAcggg                                                                                                  >  2:19246/1‑148 (MQ=255)
                                            gACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGAcggg                                                                                                  >  1:2920041/1‑148 (MQ=255)
                                            gACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGAcggg                                                                                                  >  1:2091467/1‑148 (MQ=255)
                                            gACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGAcggg                                                                                                  <  1:1405470/148‑1 (MQ=255)
                                             aCGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGAcgggc                                                                                                 >  2:1720773/1‑148 (MQ=255)
                                                 cccGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGt                                                                                             <  2:2229171/148‑1 (MQ=255)
                                                 cccGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGt                                                                                             <  2:2229182/148‑1 (MQ=255)
                                                 cccGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGt                                                                                             <  2:88362/148‑1 (MQ=255)
                                                  ccGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTg                                                                                            >  1:2407026/1‑148 (MQ=255)
                                                     cggcgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACa                                                                                         >  1:1845044/1‑148 (MQ=255)
                                                     cggcgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACa                                                                                         >  1:76871/1‑148 (MQ=255)
                                                      ggcgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAg                                                                                        >  1:3023157/1‑148 (MQ=255)
                                                        cgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGaa                                                                                      >  1:2762981/1‑148 (MQ=255)
                                                           aCTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAg                                                                                   <  2:2552709/148‑1 (MQ=255)
                                                           aCTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAg                                                                                   <  1:307883/148‑1 (MQ=255)
                                                              ggcggcGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTa                                                                                <  2:3614263/148‑1 (MQ=255)
                                                                cggcggTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTaaa                                                                              >  1:866803/1‑148 (MQ=255)
                                                                 ggcggTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAAc                                                                             <  2:1468443/148‑1 (MQ=255)
                                                                  gcggTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAAcc                                                                            >  2:2901089/1‑148 (MQ=255)
                                                                      tgtATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGc                                                                        >  2:3398418/1‑148 (MQ=255)
                                                                      tgtATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGc                                                                        >  2:1075512/1‑148 (MQ=255)
                                                                       gtATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCTACGGGCGGGTGACAGAACAGCTAAACCCGGCa                                                                       <  2:3695505/148‑1 (MQ=255)
                                                                       gtATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCa                                                                       <  1:140211/148‑1 (MQ=255)
                                                                       gtATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCa                                                                       >  2:1877041/1‑148 (MQ=255)
                                                                       gtATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCa                                                                       <  2:1258407/148‑1 (MQ=255)
                                                                        tATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAg                                                                      >  2:1464944/1‑148 (MQ=255)
                                                                        tATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAg                                                                      >  1:736860/1‑148 (MQ=255)
                                                                        tATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAg                                                                      >  1:736902/1‑148 (MQ=255)
                                                                         aTGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAgg                                                                     >  2:112838/1‑148 (MQ=255)
                                                                          tGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGGTAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGc                                                                    >  2:2360973/1‑148 (MQ=255)
                                                                          tGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGc                                                                    >  2:2753911/1‑148 (MQ=255)
                                                                          tGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGc                                                                    <  1:1964961/148‑1 (MQ=255)
                                                                           gACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCt                                                                   <  2:2760005/148‑1 (MQ=255)
                                                                            aCCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCtt                                                                  <  2:2689688/148‑1 (MQ=255)
                                                                                cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGt                                                              <  1:1920770/148‑1 (MQ=255)
                                                                                cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGt                                                              >  1:3695934/1‑148 (MQ=255)
                                                                                cagcaATACACAGGTGCGCAGATTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGt                                                              <  1:147191/148‑1 (MQ=255)
                                                                                 agcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGtt                                                             <  1:1783716/148‑1 (MQ=255)
                                                                                   caATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTac                                                           <  2:917272/148‑1 (MQ=255)
                                                                                    aaTACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTaca                                                          <  2:1757613/148‑1 (MQ=255)
                                                                                      tACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACg                                                        >  1:1983303/1‑148 (MQ=255)
                                                                                       acacAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTCAACCCGGCAGGCTTCAGTTACACGt                                                       >  2:1643731/1‑148 (MQ=255)
                                                                                         acaGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAt                                                     >  2:935866/1‑148 (MQ=37)
                                                                                         acaGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAt                                                     <  1:1342698/148‑1 (MQ=37)
                                                                                          caGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATc                                                    >  1:3386436/1‑148 (MQ=37)
                                                                                          caGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATc                                                    >  2:3544602/1‑148 (MQ=37)
                                                                                          caGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATc                                                    <  2:2091467/148‑1 (MQ=37)
                                                                                            ggTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAg                                                  >  2:2461469/1‑148 (MQ=37)
                                                                                             gTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGt                                                 <  2:1415544/148‑1 (MQ=37)
                                                                                              tGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTa                                                <  1:3773811/148‑1 (MQ=37)
                                                                                              tGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTa                                                <  1:1882017/148‑1 (MQ=37)
                                                                                                cgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATg                                              <  1:576804/148‑1 (MQ=37)
                                                                                                cgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATg                                              <  1:1551992/148‑1 (MQ=37)
                                                                                                    gCTTTACTTACGATGGTAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGaa                                          <  1:2360973/148‑1 (MQ=255)
                                                                                                    gCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGaa                                          <  1:2753911/148‑1 (MQ=37)
                                                                                                    gCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGaa                                          <  1:1083419/148‑1 (MQ=37)
                                                                                                    gCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGaa                                          <  1:2753936/148‑1 (MQ=37)
                                                                                                    gCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGaa                                          <  2:2554177/148‑1 (MQ=37)
                                                                                                     cTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAg                                         >  1:1103764/1‑148 (MQ=38)
                                                                                                     cTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAg                                         >  2:3806531/1‑148 (MQ=38)
                                                                                                       ttacttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGa                                       >  2:146534/1‑148 (MQ=255)
                                                                                                       ttacttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGa                                       <  2:3037971/148‑1 (MQ=255)
                                                                                                            tacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGca                                  >  2:1752103/1‑148 (MQ=255)
                                                                                                               gatgatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatca                               <  1:1814187/148‑1 (MQ=255)
                                                                                                               gatgatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatca                               <  1:3357130/148‑1 (MQ=255)
                                                                                                                atgatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcac                              >  1:2417519/1‑148 (MQ=255)
                                                                                                                 tgatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacc                             >  1:3164029/1‑148 (MQ=255)
                                                                                                                  gatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacca                            >  1:2079423/1‑148 (MQ=255)
                                                                                                                   atAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccat                           >  1:3227285/1‑148 (MQ=255)
                                                                                                                   atAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccat                           >  1:2227788/1‑148 (MQ=255)
                                                                                                                        tACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACCCGTATCAGTATGAGAAGGACCGCATCACCATCACCg                      <  1:1537310/148‑1 (MQ=255)
                                                                                                                              ggCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGcc                <  2:3326542/148‑1 (MQ=255)
                                                                                                                                  ggATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGaa            <  2:707259/148‑1 (MQ=255)
                                                                                                                                    aTGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAcc          >  1:3189347/1‑148 (MQ=255)
                                                                                                                                       gtggCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAccgcc       >  1:492494/1‑148 (MQ=255)
                                                                                                                                          gcgcACCTTCACACGGGCCGACCGTAAATCTGTTACCGTTACGACAGCGACGGGAGGTTTACAGAACAGATAAACACGGCAGGCTAAAGTTACACGTATCAGTATGAGACGGACCGCATCACCATCACCGACAGACTGAAACGCCGTg    >  2:1642272/1‑148 (MQ=14)
                                                                                                                                           cgcACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGa   <  1:2901089/148‑1 (MQ=255)
                                                                                                                                            gcACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGaa  >  1:2128316/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAA  >  CP009273/725942‑726229

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