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,9730CA96.5% 695.2 / 14.2 228P93P (CCG→CCT) ybgDputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base A (95/125);  total (97/131)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 228 total aligned reads displayed.

CTGAACCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGA  >  CP009273/747828‑748117
                                                                                                                                                 |                                                                                                                                                
  gAACCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCa                                                                                                                                              <  1:1296286/148‑1 (MQ=255)
  gAACCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCa                                                                                                                                              >  2:2630831/1‑148 (MQ=255)
       tAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCg                                                                                                                                         >  1:34133/1‑148 (MQ=255)
        aaTACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGc                                                                                                                                        >  2:3156258/1‑148 (MQ=255)
        aaTACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGc                                                                                                                                        >  2:3244392/1‑148 (MQ=255)
         aTACGCTGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCt                                                                                                                                       <  2:3108149/148‑1 (MQ=255)
         aTACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCt                                                                                                                                       <  2:3400616/148‑1 (MQ=255)
           aCGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTAc                                                                                                                                     <  1:1565482/148‑1 (MQ=255)
              aTGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGt                                                                                                                                  <  2:1285724/148‑1 (MQ=255)
                 ttACCGTCATTTTTGTCCATCAGTCGAACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGt                                                                                                                               >  2:1636551/1‑148 (MQ=255)
                    ccGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAAATTTGGATACAGGCATTCACCTACCGTTGTAAg                                                                                                                            >  1:3094940/1‑148 (MQ=255)
                     cGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCTGCTACCGTTGTCAGa                                                                                                                           <  2:2770533/148‑1 (MQ=255)
                        cATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgc                                                                                                                        <  1:55275/148‑1 (MQ=255)
                        cATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgc                                                                                                                        >  1:2434519/1‑148 (MQ=255)
                        cATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgc                                                                                                                        <  1:2070779/148‑1 (MQ=255)
                        cATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgc                                                                                                                        >  1:1676149/1‑148 (MQ=255)
                         aTTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgca                                                                                                                       <  2:603872/148‑1 (MQ=255)
                           ttttGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgcagg                                                                                                                     <  2:662572/148‑1 (MQ=255)
                                   aTCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATcac                                                                                                             >  1:3514446/1‑148 (MQ=255)
                                   aTCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATcac                                                                                                             >  2:451212/1‑148 (MQ=255)
                                     cAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAg                                                                                                           >  2:1291332/1‑148 (MQ=255)
                                     cAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAg                                                                                                           <  2:652346/148‑1 (MQ=255)
                                           tACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATc                                                                                                     >  1:2990490/1‑148 (MQ=255)
                                           tACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATc                                                                                                     >  2:2205143/1‑148 (MQ=255)
                                           tACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATc                                                                                                     >  2:1699517/1‑148 (MQ=255)
                                            acacCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCa                                                                                                    >  1:514948/1‑148 (MQ=255)
                                             cacCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCaa                                                                                                   <  2:2544928/148‑1 (MQ=255)
                                              acCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAg                                                                                                  <  1:766290/148‑1 (MQ=255)
                                              acCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAg                                                                                                  >  1:1824616/1‑148 (MQ=255)
                                               ccGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGc                                                                                                 <  2:210861/148‑1 (MQ=255)
                                               ccGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGc                                                                                                 >  1:1702045/1‑148 (MQ=255)
                                               ccGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGc                                                                                                 <  2:211482/148‑1 (MQ=255)
                                                cGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTCTCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCg                                                                                                >  2:149224/1‑148 (MQ=255)
                                                cGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCg                                                                                                >  1:1846525/1‑148 (MQ=255)
                                                  aCCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGaa                                                                                              <  1:2605789/148‑1 (MQ=255)
                                                       aGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCg                                                                                         <  1:1380603/148‑1 (MQ=255)
                                                       aGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCg                                                                                         <  2:2026162/148‑1 (MQ=255)
                                                            cTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGc                                                                                    <  1:207661/148‑1 (MQ=255)
                                                            cTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGc                                                                                    <  1:1286848/148‑1 (MQ=255)
                                                                 ccTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACtt                                                                               >  1:3263691/1‑148 (MQ=255)
                                                                  cTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACttt                                                                              <  1:2488756/148‑1 (MQ=255)
                                                                  cTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACttt                                                                              >  2:1310894/1‑148 (MQ=255)
                                                                    gCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTAt                                                                            <  1:2549487/148‑1 (MQ=255)
                                                                    gCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTAt                                                                            <  2:852362/148‑1 (MQ=255)
                                                                       cTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCa                                                                                  <  2:3115293/139‑1 (MQ=255)
                                                                         ggTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTa                                                                       <  1:2016042/148‑1 (MQ=255)
                                                                         ggTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTa                                                                       <  2:1733246/148‑1 (MQ=255)
                                                                                  caacaaAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTcc                                                              >  1:68673/1‑148 (MQ=255)
                                                                                      aaaGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGt                                                          <  2:34133/148‑1 (MQ=255)
                                                                                        aGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGttt                                                        <  1:1349562/148‑1 (MQ=255)
                                                                                        aGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGttt                                                        >  1:978621/1‑148 (MQ=255)
                                                                                         ggCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTta                                                       >  1:1780219/1‑148 (MQ=255)
                                                                                          gCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTtat                                                      >  1:2144672/1‑148 (MQ=255)
                                                                                          gCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTtat                                                      >  1:2144820/1‑148 (MQ=255)
                                                                                           cGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTtata                                                     <  2:2459508/148‑1 (MQ=255)
                                                                                             tAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATg                                                   <  2:2679816/148‑1 (MQ=255)
                                                                                                cACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGt                                                >  1:1002825/1‑148 (MQ=255)
                                                                                                cACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGt                                                >  2:3166352/1‑148 (MQ=255)
                                                                                                   cAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGt                                             >  2:3486765/1‑148 (MQ=255)
                                                                                                         tCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGAc                                       <  2:1539247/148‑1 (MQ=255)
                                                                                                         tCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGAc                                       <  1:454490/148‑1 (MQ=255)
                                                                                                          cTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACt                                      <  2:2434525/148‑1 (MQ=255)
                                                                                                           ttAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACtt                                     >  1:1316135/1‑148 (MQ=255)
                                                                                                             aGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGt                                   <  2:3542783/148‑1 (MQ=255)
                                                                                                              gCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTc                                  <  1:36830/148‑1 (MQ=255)
                                                                                                               ccGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTcc                                 >  1:599412/1‑148 (MQ=255)
                                                                                                                cGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCa                                >  2:1029135/1‑148 (MQ=255)
                                                                                                                  tGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAg                              >  1:770392/1‑148 (MQ=255)
                                                                                                                   gCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGa                             <  1:1504977/148‑1 (MQ=255)
                                                                                                                    cTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGAt                            <  1:573678/148‑1 (MQ=255)
                                                                                                                    cTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATNAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGAt                            <  1:1809247/148‑1 (MQ=255)
                                                                                                                    cTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGAt                            <  1:2781939/148‑1 (MQ=255)
                                                                                                                     tATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATc                           <  2:3550076/148‑1 (MQ=255)
                                                                                                                      aTCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCt                          <  1:1722971/148‑1 (MQ=255)
                                                                                                                      aTCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCt                          <  1:1836800/148‑1 (MQ=255)
                                                                                                                           aGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGtt                     <  1:1048120/148‑1 (MQ=255)
                                                                                                                           aGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGtt                     <  1:1476047/148‑1 (MQ=255)
                                                                                                                            ggTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGttt                    <  1:2836462/148‑1 (MQ=255)
                                                                                                                               tACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTa                 <  2:2266445/148‑1 (MQ=255)
                                                                                                                               tACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTa                 <  1:3244392/148‑1 (MQ=255)
                                                                                                                                  gCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCa              <  2:2634847/148‑1 (MQ=255)
                                                                                                                                  gCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCa              <  1:2205143/148‑1 (MQ=255)
                                                                                                                                  gCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCa              <  2:1345114/148‑1 (MQ=255)
                                                                                                                                   ccAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCaa             <  2:930552/148‑1 (MQ=255)
                                                                                                                                   ccAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCaa             >  1:2889459/1‑148 (MQ=255)
                                                                                                                                   ccAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCaa             <  2:514948/148‑1 (MQ=255)
                                                                                                                                   ccAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCaa             >  1:2002401/1‑148 (MQ=255)
                                                                                                                                    cAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAat            >  1:2640024/1‑148 (MQ=255)
                                                                                                                                    cAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAat            <  2:1988394/148‑1 (MQ=255)
                                                                                                                                     aaCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAata           >  1:3369118/1‑148 (MQ=255)
                                                                                                                                     aaCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAata           <  2:2859504/148‑1 (MQ=255)
                                                                                                                                     aaCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAata           <  2:732804/148‑1 (MQ=255)
                                                                                                                                     aaCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAata           <  2:2520696/148‑1 (MQ=255)
                                                                                                                                      aCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAatat          <  1:2538895/148‑1 (MQ=255)
                                                                                                                                       cTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtc         <  1:3305874/148‑1 (MQ=255)
                                                                                                                                        tttGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtct        >  2:157972/1‑148 (MQ=255)
                                                                                                                                         ttGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtctt       <  2:1977395/148‑1 (MQ=255)
                                                                                                                                           ggATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtcttct     <  2:1067458/148‑1 (MQ=255)
                                                                                                                                           ggATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtcttct     <  1:3323802/148‑1 (MQ=255)
                                                                                                                                            gATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTg    >  2:947885/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTGAACCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGA  >  CP009273/747828‑748117

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