Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,385,884 G→T 11.8% S1245* (TCG→TAG)  yhdP ← DUF3971‑AsmA2 domains protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,385,8840GT11.8% 642.1 / 43.8 194S1245* (TCG→TAG) yhdPDUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base G (87/84);  new base T (16/7);  total (103/91)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 195 total aligned reads displayed.

CCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGCGACAACCGCTT  >  CP009273/3385738‑3386031
                                                                                                                                                  |                                                                                                                                                   
ccTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTAAGCGGTTGGCGCAACACTTCGTTGATTTGAGGATCGTCAAGCGGACCCGa                                                                                                                                                    >  2:2830235/1‑148 (MQ=255)
 cTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGaa                                                                                                                                                   <  1:1209266/148‑1 (MQ=255)
 cTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGaa                                                                                                                                                   <  1:3059528/148‑1 (MQ=255)
  tACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGaaa                                                                                                                                                  >  1:2394251/1‑148 (MQ=255)
    cGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAAtg                                                                                                                                                <  2:2893636/148‑1 (MQ=255)
       ttATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGa                                                                                                                                             >  1:3036259/1‑148 (MQ=255)
         aTGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCTAAATGTGATa                                                                                                                                           <  2:371066/148‑1 (MQ=255)
            ttCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAgcg                                                                                                                                        >  2:2974964/1‑148 (MQ=255)
             tCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAgcgc                                                                                                                                       >  1:461337/1‑148 (MQ=255)
                ggCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCaaa                                                                                                                                    <  1:695044/148‑1 (MQ=255)
                 gCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCaaaa                                                                                                                                   <  2:3213665/148‑1 (MQ=255)
                  cTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAAt                                                                                                                                  >  1:2828845/1‑148 (MQ=255)
                  cTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAAt                                                                                                                                  >  1:2073571/1‑148 (MQ=255)
                   tCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATg                                                                                                                                 <  2:2394537/148‑1 (MQ=255)
                    cGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATgg                                                                                                                                <  2:913499/148‑1 (MQ=255)
                      ttGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGgag                                                                                                                              <  1:1449767/148‑1 (MQ=255)
                      ttGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGgag                                                                                                                              <  2:2237937/148‑1 (MQ=255)
                      ttGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGgag                                                                                                                              <  2:590222/148‑1 (MQ=255)
                       tGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGgaga                                                                                                                             <  2:2014078/148‑1 (MQ=255)
                       tGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGgaga                                                                                                                             >  1:1435803/1‑148 (MQ=255)
                       tGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGgaga                                                                                                                             <  2:2013765/148‑1 (MQ=255)
                         gCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTAAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCTAAATGTGATAGCGCAAAATGGAGACt                                                                                                                           <  1:1273479/148‑1 (MQ=255)
                         gCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTAAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCTAAATGTGATAGCGCAAAATGGAGACt                                                                                                                           <  1:2256112/148‑1 (MQ=255)
                            tttGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTg                                                                                                                        <  1:2150361/148‑1 (MQ=255)
                             ttGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGc                                                                                                                       >  1:1962134/1‑148 (MQ=255)
                             ttGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGc                                                                                                                       >  2:2396970/1‑148 (MQ=255)
                             ttGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTAAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCTAAATGTGATAGCGCAAAATGGAGACTTTGc                                                                                                                       >  2:777731/1‑148 (MQ=255)
                                gCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTcc                                                                                                                    >  2:2398003/1‑148 (MQ=255)
                                     gATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGc                                                                                                               <  2:420552/148‑1 (MQ=255)
                                       ttATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCgg                                                                                                             <  1:1729533/148‑1 (MQ=255)
                                          tCCTATGAGATTGGGGCAATTACGCGCCCTCGTAAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCTAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGccc                                                                                                          >  1:2860204/1‑148 (MQ=255)
                                             tATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAg                                                                                                       >  2:34381/1‑148 (MQ=255)
                                              aTGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGc                                                                                                      >  2:2709710/1‑148 (MQ=255)
                                              aTGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGc                                                                                                      <  1:1611334/148‑1 (MQ=255)
                                                 agagTGGGGCAATTCCGCGCCCTCGTCAAATCATTGCGCTGTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACt                                                                                                   <  2:1747847/148‑1 (MQ=255)
                                                  gaTTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACtt                                                                                                  <  2:1171909/148‑1 (MQ=255)
                                                    ttGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTa                                                                                                >  1:2863005/1‑148 (MQ=255)
                                                     tGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTAc                                                                                               >  2:3156119/1‑148 (MQ=255)
                                                       gggCAATTACGCGCCCTCGTAAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCTAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTg                                                                                             >  1:3030459/1‑148 (MQ=255)
                                                       gggCAATTACGCGCCCTCGTAAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCTAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTg                                                                                             >  2:1161686/1‑148 (MQ=255)
                                                          cAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTggcg                                                                                          <  2:318818/148‑1 (MQ=255)
                                                          cAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTggcg                                                                                          >  2:1366810/1‑148 (MQ=255)
                                                             ttACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCa                                                                                       >  1:3218665/1‑148 (MQ=255)
                                                              tACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCaa                                                                                      <  2:1279950/148‑1 (MQ=255)
                                                              tACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCaa                                                                                      >  2:2307269/1‑148 (MQ=255)
                                                                    cccTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTg                                                                                <  1:3208234/148‑1 (MQ=255)
                                                                       tagTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGACCCAGAACTTTACTGGCGGCAAACACTGCcg                                                                             >  2:2488122/3‑148 (MQ=255)
                                                                       tCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCcg                                                                             >  1:184043/1‑148 (MQ=255)
                                                                       tCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTAATGGCGGCAAACACTGcct                                                                             >  1:689401/1‑147 (MQ=255)
                                                                        cGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCcgc                                                                            >  2:1228837/1‑148 (MQ=255)
                                                                        cGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCcgc                                                                            >  2:2259850/1‑148 (MQ=255)
                                                                          tCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCgcgc                                                                          <  2:1681482/148‑1 (MQ=255)
                                                                            aaaTCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCca                                                                        >  2:2916532/1‑148 (MQ=255)
                                                                             aaTCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCcaa                                                                       <  2:2081656/148‑1 (MQ=255)
                                                                               tCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCcaaca                                                                     <  1:319174/148‑1 (MQ=255)
                                                                                 aTTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAAt                                                                   <  2:1328471/148‑1 (MQ=255)
                                                                                 aTTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAAt                                                                   <  2:461337/148‑1 (MQ=255)
                                                                                 aTTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAAt                                                                   <  2:103691/148‑1 (MQ=255)
                                                                                  ttGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATg                                                                  >  1:3118474/1‑148 (MQ=255)
                                                                                   tGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCTAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATgg                                                                 >  2:70090/1‑148 (MQ=255)
                                                                                             cTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGc                                                       >  2:304869/1‑148 (MQ=255)
                                                                                                  cgcgGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAg                                                  <  1:184312/148‑1 (MQ=255)
                                                                                                   gcgGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGc                                                 <  2:2868262/148‑1 (MQ=255)
                                                                                                    cgGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGcc                                                >  2:2977218/1‑148 (MQ=255)
                                                                                                        tGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAg                                            >  1:1093088/1‑148 (MQ=255)
                                                                                                              aaCACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGc                                      <  1:1008295/148‑1 (MQ=255)
                                                                                                               acacTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGcc                                     >  2:1896890/1‑148 (MQ=255)
                                                                                                               acacTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGcc                                     <  1:1952488/148‑1 (MQ=255)
                                                                                                                 acTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCa                                   >  1:2748103/1‑148 (MQ=255)
                                                                                                                   ttCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACTGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCAcc                                 >  2:2367408/1‑148 (MQ=255)
                                                                                                                     cGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGt                               <  1:1654852/148‑1 (MQ=255)
                                                                                                                       ttGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCg                             >  2:1611837/1‑148 (MQ=255)
                                                                                                                        tGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGc                            >  2:554098/1‑148 (MQ=255)
                                                                                                                        tGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGc                            >  2:1571626/1‑148 (MQ=255)
                                                                                                                         gATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCa                           <  1:816415/148‑1 (MQ=255)
                                                                                                                            ttGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGaa                        >  2:437472/1‑148 (MQ=255)
                                                                                                                            ttGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGaa                        >  2:2100357/1‑148 (MQ=255)
                                                                                                                            ttGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGaa                        >  1:3157059/1‑148 (MQ=255)
                                                                                                                             tGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGaaa                       <  1:1436449/148‑1 (MQ=255)
                                                                                                                              gCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAAt                      <  1:2398003/148‑1 (MQ=255)
                                                                                                                                 gATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAAtctc                   <  2:613321/148‑1 (MQ=255)
                                                                                                                                 gATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAAtctc                   <  1:1366810/148‑1 (MQ=255)
                                                                                                                                  aTCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATctct                  >  1:1883484/1‑148 (MQ=255)
                                                                                                                                     gTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGt               <  1:2238449/148‑1 (MQ=255)
                                                                                                                                     gTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGt               <  1:2238497/148‑1 (MQ=255)
                                                                                                                                       ccAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGc             >  2:290673/1‑148 (MQ=255)
                                                                                                                                       ccAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGc             >  1:293550/1‑148 (MQ=255)
                                                                                                                                       ccAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGc             >  2:2646601/1‑148 (MQ=255)
                                                                                                                                       ccAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGc             >  1:2407041/1‑148 (MQ=255)
                                                                                                                                         aGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGCGa           <  1:2448534/148‑1 (MQ=255)
                                                                                                                                           cGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGCGACa         <  1:272691/148‑1 (MQ=255)
                                                                                                                                           cGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGCGACa         <  1:2916532/148‑1 (MQ=255)
                                                                                                                                             gACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGCGACAAc       >  2:191371/1‑148 (MQ=255)
                                                                                                                                               cccTAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGCGACAACCg     >  2:668672/1‑148 (MQ=255)
                                                                                                                                               cccTAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGCGACAACCg     >  2:668807/1‑148 (MQ=255)
                                                                                                                                               cccGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGCGACAACCg     <  1:2972704/148‑1 (MQ=255)
                                                                                                                                                ccGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGCGACAACCGc    <  2:2146534/148‑1 (MQ=255)
                                                                                                                                                ccGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGCGACAACCGc    <  1:1026203/148‑1 (MQ=255)
                                                                                                                                                 cGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGCGACAACCGCt   <  1:1309040/148‑1 (MQ=255)
                                                                                                                                                  gAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGCGACAACCGCtt  >  2:56641/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTTCGTTGATTTGCGGATCGTCCAGCGGACCCGAAATGTGATAGCGCAAAATGGAGACTTTGCTCCACAGCGGCCCCAGCACTTTACTGGCGGCAAACACTGCCGCGCCAACAATGGGGTTAACCGCAAAAGCCGCAGCCACGCCCACCGTCGCAGAAATCTCTGGTGCGACAACCGCTT  >  CP009273/3385738‑3386031

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: