Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 747,973 C→A 100% P93P (CCG→CCT)  ybgD ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273747,9730CA97.9% 772.5 / 4.0 240P93P (CCG→CCT) ybgDputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base A (125/110);  total (127/113)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 240 total aligned reads displayed.

CTGAACCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGAT  >  CP009273/747828‑748118
                                                                                                                                                 |                                                                                                                                                 
    aCCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTc                                                                                                                                            >  1:3324877/1‑149 (MQ=255)
     ccTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTcc                                                                                                                                           >  1:2061440/1‑149 (MQ=255)
      cTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCg                                                                                                                                          >  1:3578583/1‑149 (MQ=255)
       tAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGc                                                                                                                                         <  1:3790555/149‑1 (MQ=255)
       tAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGc                                                                                                                                         >  1:739564/1‑149 (MQ=255)
       tAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGc                                                                                                                                         >  2:3129849/1‑149 (MQ=255)
           aCGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTAcc                                                                                                                                     >  2:3034210/1‑149 (MQ=255)
               tGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGGTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTg                                                                                                                                 >  2:3775675/1‑149 (MQ=255)
                    ccGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGa                                                                                                                            <  1:3053752/149‑1 (MQ=255)
                       tCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGTGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgc                                                                                                                         <  1:3016668/149‑1 (MQ=255)
                       tCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgc                                                                                                                         >  1:492148/1‑149 (MQ=255)
                         aTTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgcag                                                                                                                       <  2:485213/149‑1 (MQ=255)
                         aTTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgcag                                                                                                                       <  1:2913697/149‑1 (MQ=255)
                         aTTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgcag                                                                                                                       <  1:1903657/149‑1 (MQ=255)
                         aTTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgcag                                                                                                                       <  2:3608103/149‑1 (MQ=255)
                          tttttGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgcagg                                                                                                                      <  1:2740289/149‑1 (MQ=255)
                          tttttGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgcagg                                                                                                                      >  1:2432050/1‑149 (MQ=255)
                           ttttGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgcaggc                                                                                                                     <  1:338056/149‑1 (MQ=255)
                               gTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGAt                                                                                                                 >  2:3474578/1‑149 (MQ=255)
                                    tCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAg                                                                                                            >  1:1470399/1‑149 (MQ=255)
                                    tCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAg                                                                                                            >  1:1139984/1‑149 (MQ=255)
                                    tCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAg                                                                                                            >  2:57163/1‑149 (MQ=255)
                                     cAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGt                                                                                                           >  2:3179332/1‑149 (MQ=255)
                                      aGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGtt                                                                                                          <  2:1661449/149‑1 (MQ=255)
                                      aGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGtt                                                                                                          >  2:3299726/1‑149 (MQ=255)
                                       gTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTg                                                                                                         >  1:1619942/1‑149 (MQ=255)
                                       gTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTg                                                                                                         >  1:1416510/1‑149 (MQ=255)
                                         cGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGAt                                                                                                       <  1:826245/149‑1 (MQ=255)
                                          gTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATc                                                                                                      <  2:1463860/149‑1 (MQ=255)
                                            acacCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCaa                                                                                                    <  1:1487521/149‑1 (MQ=255)
                                                 gACCCCCGGTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGaa                                                                                               <  1:2121011/149‑1 (MQ=255)
                                                 gACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGac                                                                                               >  1:4078453/1‑148 (MQ=255)
                                                 gACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGaa                                                                                               >  1:581307/1‑149 (MQ=255)
                                                 gACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGaa                                                                                               >  2:2199340/1‑149 (MQ=255)
                                                  aCCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAAt                                                                                              >  2:2704829/1‑149 (MQ=255)
                                                  aCCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAAt                                                                                              <  1:2305500/149‑1 (MQ=255)
                                                    cccAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGt                                                                                            <  2:2622851/149‑1 (MQ=255)
                                                    cccAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTAACGTTGTCAGAAGCAGGCAGATCACATTTGATCAAGCGAATGt                                                                                            >  1:1064081/1‑149 (MQ=255)
                                                       aGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGa                                                                                         <  2:2907870/149‑1 (MQ=255)
                                                        gTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGAc                                                                                        >  2:3725766/1‑149 (MQ=255)
                                                             ttCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCa                                                                                   >  2:124066/1‑149 (MQ=255)
                                                               cGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACt                                                                                 <  2:3346020/149‑1 (MQ=255)
                                                                gCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACtt                                                                                >  2:2536455/1‑149 (MQ=255)
                                                                gCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACtt                                                                                <  2:1494/149‑1 (MQ=255)
                                                                gCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACtt                                                                                <  2:3955862/149‑1 (MQ=255)
                                                                   tGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTAt                                                                             <  2:3924170/149‑1 (MQ=255)
                                                                   tGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTAt                                                                             <  2:3924249/149‑1 (MQ=255)
                                                                    gCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTAtt                                                                            <  2:1143871/149‑1 (MQ=255)
                                                                        tGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTa                                                                        >  1:3556293/1‑149 (MQ=255)
                                                                        tGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTa                                                                        >  1:931623/1‑149 (MQ=255)
                                                                           tgtTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatga                                                                     <  2:883637/149‑1 (MQ=255)
                                                                            gtTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCATGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatgat                                                                    <  2:257643/149‑1 (MQ=255)
                                                                            gtTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatgat                                                                    <  2:2875290/149‑1 (MQ=255)
                                                                            gtTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatgat                                                                    <  1:2416279/149‑1 (MQ=255)
                                                                                 tCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTcc                                                               <  1:2231651/149‑1 (MQ=255)
                                                                                 tCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTcc                                                               <  1:2418644/149‑1 (MQ=255)
                                                                                  caacaaAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTccc                                                              <  1:618750/149‑1 (MQ=255)
                                                                                  caacaaAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTccc                                                              <  2:2429653/149‑1 (MQ=255)
                                                                                     caaAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGt                                                           <  2:1422/149‑1 (MQ=255)
                                                                                       aaGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGttt                                                         <  2:1113014/149‑1 (MQ=255)
                                                                                        aGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAATCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTta                                                        >  2:2695304/1‑149 (MQ=255)
                                                                                        aGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTta                                                        >  2:3803708/1‑149 (MQ=255)
                                                                                         ggCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTtat                                                       >  1:3818499/1‑149 (MQ=255)
                                                                                            gTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATg                                                    >  1:3744124/1‑149 (MQ=255)
                                                                                               gCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGt                                                 >  1:1232053/1‑149 (MQ=255)
                                                                                               gCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGt                                                 <  1:3107005/149‑1 (MQ=255)
                                                                                               gCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGt                                                 >  1:2591551/1‑149 (MQ=255)
                                                                                               gCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGt                                                 <  1:2452727/149‑1 (MQ=255)
                                                                                                cACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGtt                                                <  1:3228336/149‑1 (MQ=255)
                                                                                                  ccAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGt                                              >  2:3583042/1‑149 (MQ=255)
                                                                                                  ccAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGt                                              >  1:408708/1‑149 (MQ=255)
                                                                                                  ccAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGt                                              >  1:822432/1‑149 (MQ=255)
                                                                                                    aGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCg                                            <  1:1685321/149‑1 (MQ=255)
                                                                                                     gTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGt                                           >  1:642407/1‑149 (MQ=255)
                                                                                                     gTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGt                                           <  2:4082132/149‑1 (MQ=255)
                                                                                                      ttGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTg                                          >  2:2035463/1‑149 (MQ=255)
                                                                                                       tGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGa                                         >  1:1035234/1‑149 (MQ=255)
                                                                                                       tGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGa                                         <  1:3474578/149‑1 (MQ=255)
                                                                                                           ttAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGTAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTg                                     <  2:641591/149‑1 (MQ=255)
                                                                                                                  tGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGa                              <  2:1729625/149‑1 (MQ=255)
                                                                                                                      aTCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTa                          <  2:802879/149‑1 (MQ=255)
                                                                                                                      aTCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTa                          <  1:959930/149‑1 (MQ=255)
                                                                                                                      aTCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTa                          <  1:960084/149‑1 (MQ=255)
                                                                                                                       tCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTAt                         >  2:1093693/1‑149 (MQ=255)
                                                                                                                          aaGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGtt                      <  1:3641103/149‑1 (MQ=255)
                                                                                                                              ttACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTa                  <  1:2777341/149‑1 (MQ=255)
                                                                                                                                aCGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATc                <  2:3830526/149‑1 (MQ=255)
                                                                                                                                  gCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCaa              >  1:2953870/1‑149 (MQ=255)
                                                                                                                                   ccAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAat             >  1:709542/1‑149 (MQ=255)
                                                                                                                                   ccAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAat             <  1:1498338/149‑1 (MQ=255)
                                                                                                                                   ccAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAat             <  1:868592/149‑1 (MQ=255)
                                                                                                                                   ccAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAat             >  1:3928435/1‑149 (MQ=255)
                                                                                                                                   ccAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAat             <  1:905627/149‑1 (MQ=255)
                                                                                                                                    cAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGT‑CAAGATCTATGTTTTTATCAatat           <  2:3315452/149‑1 (MQ=255)
                                                                                                                                     aaCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCCGCGCCCAGTTGATCAAGCCAAAGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAatat           >  1:4034753/1‑149 (MQ=255)
                                                                                                                                     aaCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAatat           >  2:1750283/1‑149 (MQ=255)
                                                                                                                                     aaCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAatat           >  2:601667/1‑149 (MQ=255)
                                                                                                                                      aCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtc          >  2:2434968/1‑149 (MQ=255)
                                                                                                                                        tttGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtctt        <  1:2835895/149‑1 (MQ=255)
                                                                                                                                           ggATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTg     >  2:2061814/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTGAACCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGAT  >  CP009273/747828‑748118

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