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,9730CA95.7% 630.3 / 17.5 209P93P (CCG→CCT) ybgDputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base C (3/6);  new base A (80/120);  total (83/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 210 total aligned reads displayed.

CTGAACCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGA  >  CP009273/747828‑748120
                                                                                                                                                 |                                                                                                                                                   
cTGAACCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAgg                                                                                                                                                   <  2:1173137/148‑1 (MQ=255)
   aaCCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCAt                                                                                                                                                >  2:2609350/1‑148 (MQ=255)
    aCCTAATACTATGTTACCTTCATTTTTGTACATCAGTCGTACACCGCCCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAAAAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCAAACTTTGGATACAGGCAtt                                                                                                                                               <  1:2442154/148‑1 (MQ=255)
     ccTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTc                                                                                                                                              <  1:800353/148‑1 (MQ=255)
     ccTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTc                                                                                                                                              <  1:2981819/148‑1 (MQ=255)
        aaTACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGc                                                                                                                                           >  1:1211930/1‑148 (MQ=255)
        aaTACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGc                                                                                                                                           <  2:1637483/148‑1 (MQ=255)
           aCGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTAc                                                                                                                                        <  1:2431489/148‑1 (MQ=255)
           aCGATGTTACCGACATTTTTGACAATCAGTCGTACAACGACCCCAGTTCCTTCGCCTGCACTGGTGTTGATCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTAACTACAACTTTGGATACAGGCAATCAGCTAc                                                                                                                                        >  2:1923431/1‑148 (MQ=255)
              aTGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGt                                                                                                                                     <  1:3449994/148‑1 (MQ=255)
                 ttACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGt                                                                                                                                  <  2:3546280/148‑1 (MQ=255)
                 ttACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGt                                                                                                                                  <  2:3546341/148‑1 (MQ=255)
                       tCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAg                                                                                                                            <  2:2808175/148‑1 (MQ=255)
                       tCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAg                                                                                                                            >  2:1851481/1‑148 (MQ=255)
                        cATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgc                                                                                                                           <  1:1522463/148‑1 (MQ=255)
                            tttGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAgcaggc                                                                                                                       >  2:2461527/1‑148 (MQ=255)
                                tCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGAt                                                                                                                   >  1:150450/1‑148 (MQ=255)
                                  cATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATca                                                                                                                 >  2:3408579/1‑148 (MQ=255)
                                    tCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATcaca                                                                                                               >  2:1055163/1‑148 (MQ=255)
                                    tCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATcaca                                                                                                               >  2:2795597/1‑148 (MQ=255)
                                     cAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAg                                                                                                              <  2:441216/148‑1 (MQ=255)
                                           tACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCAGTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATc                                                                                                        <  1:1382644/148‑1 (MQ=255)
                                            acacCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCa                                                                                                       >  1:1072502/1‑148 (MQ=255)
                                              acCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACC‑TTGTCAGAAGCAGGCAGATCACAGTTGATCAAGc                                                                                                    <  2:842665/148‑1 (MQ=255)
                                                cGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCg                                                                                                   >  1:456903/1‑148 (MQ=255)
                                                 gACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGa                                                                                                  <  1:369567/148‑1 (MQ=255)
                                                  aCCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGaa                                                                                                 <  1:3523611/148‑1 (MQ=255)
                                                     ccAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGt                                                                                              <  2:2321426/148‑1 (MQ=255)
                                                      cAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTc                                                                                             >  1:1906048/1‑148 (MQ=255)
                                                       aGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCg                                                                                            >  1:1819967/1‑148 (MQ=255)
                                                       aGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCg                                                                                            <  1:433429/148‑1 (MQ=255)
                                                           gCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAATCGAAAGTCGACgg                                                                                        >  1:1093494/1‑148 (MQ=255)
                                                              tCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCa                                                                                     <  2:2335000/148‑1 (MQ=255)
                                                              tCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCa                                                                                     <  2:2056590/148‑1 (MQ=255)
                                                              tCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCa                                                                                     <  2:648995/148‑1 (MQ=255)
                                                                 ccTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACtt                                                                                  <  1:1620148/148‑1 (MQ=255)
                                                                  cTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACttt                                                                                 <  2:3056403/148‑1 (MQ=255)
                                                                  cTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACttt                                                                                 <  2:3094832/148‑1 (MQ=255)
                                                                     cACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTAtt                                                                              >  2:2573704/1‑148 (MQ=255)
                                                                      aCTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTg                                                                             >  1:701577/1‑148 (MQ=255)
                                                                       cTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGc                                                                            <  2:1722394/148‑1 (MQ=255)
                                                                        tGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCt                                                                           <  1:20229/148‑1 (MQ=255)
                                                                           tgtTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatg                                                                        <  2:1569220/148‑1 (MQ=255)
                                                                            gtTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatga                                                                       <  1:264297/148‑1 (MQ=255)
                                                                             ttGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatgat                                                                      >  2:1359957/1‑148 (MQ=255)
                                                                                 tCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTc                                                                  <  1:1893272/148‑1 (MQ=255)
                                                                                  caacaaAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTCGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTcc                                                                 <  2:1118905/148‑1 (MQ=255)
                                                                                     caaAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCgg                                                              >  2:3411571/1‑148 (MQ=255)
                                                                                       aaGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGtt                                                            <  1:887052/148‑1 (MQ=255)
                                                                                        aGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGttt                                                           <  2:335607/148‑1 (MQ=255)
                                                                                        aGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGttt                                                           >  2:1528029/1‑148 (MQ=255)
                                                                                               gCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATggg                                                    <  1:2249239/148‑1 (MQ=255)
                                                                                                cACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGt                                                   <  1:3466214/148‑1 (MQ=255)
                                                                                                  ccAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTg                                                 <  1:2287805/148‑1 (MQ=255)
                                                                                                   cAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGt                                                <  2:150450/148‑1 (MQ=255)
                                                                                                   cAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGt                                                >  2:704440/1‑148 (MQ=255)
                                                                                                       tGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTg                                            >  2:159760/1‑148 (MQ=255)
                                                                                                       tGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTGGTCGTg                                            >  2:1250823/1‑148 (MQ=255)
                                                                                                         tCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGAc                                          <  1:1284863/148‑1 (MQ=255)
                                                                                                         tCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGAc                                          >  1:1478652/1‑148 (MQ=255)
                                                                                                          cTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACt                                         <  2:701577/148‑1 (MQ=255)
                                                                                                              gCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTc                                     <  1:736687/148‑1 (MQ=255)
                                                                                                              gCCGTGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTc                                     <  1:1337072/148‑1 (MQ=255)
                                                                                                                cGTGCTATCGAAGGTTACGCCAACTTTGGAGACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCa                                   <  2:1824745/148‑1 (MQ=255)
                                                                                                                  tGCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAg                                 <  2:2317971/148‑1 (MQ=255)
                                                                                                                   gCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGa                                <  1:3332514/148‑1 (MQ=255)
                                                                                                                   gCTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGa                                <  1:2525510/148‑1 (MQ=255)
                                                                                                                    cTATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGAt                               >  1:917539/1‑148 (MQ=255)
                                                                                                                     tATCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATc                              <  1:3508468/148‑1 (MQ=255)
                                                                                                                      aTCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCt                             >  1:273967/1‑148 (MQ=255)
                                                                                                                      aTCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCt                             <  2:3025503/148‑1 (MQ=255)
                                                                                                                      aTCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCt                             <  2:1046068/148‑1 (MQ=255)
                                                                                                                      aTCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCt                             <  1:2461527/148‑1 (MQ=255)
                                                                                                                       tCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTa                            <  2:544074/148‑1 (MQ=255)
                                                                                                                       tCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTa                            <  1:3298356/148‑1 (MQ=255)
                                                                                                                       tCGAAGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTa                            >  1:1718877/1‑148 (MQ=255)
                                                                                                                        cGAAGGTTACGCCAACTTTCGATACAGGCATTCCGAAACCGTTGACAGAAGCACGCACATCACAGATGATCAAGCCAATGTCTACGGACACTTTATTCAAATGATGCTCACGGTTTATATGGTTTGTAGTTACTAGACCAAGATCAAt                           >  2:422805/1‑148 (MQ=255)
                                                                                                                         gAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATg                          <  1:2687431/148‑1 (MQ=255)
                                                                                                                         gAAGGTTACGCCAACTGTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATg                          <  2:3329635/148‑1 (MQ=255)
                                                                                                                          aaGGTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGt                         >  1:2199706/1‑148 (MQ=255)
                                                                                                                            ggTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGttt                       >  2:2256235/1‑148 (MQ=255)
                                                                                                                            ggTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGttt                       <  1:59547/148‑1 (MQ=255)
                                                                                                                            ggTTACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGttt                       <  2:139947/148‑1 (MQ=255)
                                                                                                                               tACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTa                    <  2:2460614/148‑1 (MQ=255)
                                                                                                                               tACGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTa                    <  2:3155668/148‑1 (MQ=255)
                                                                                                                               tACGCCAACTTTGGATACAGGCAGTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTa                    <  2:3155589/148‑1 (MQ=255)
                                                                                                                                 cGCCAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATc                  >  2:3169450/1‑148 (MQ=255)
                                                                                                                                   ccAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCaa                <  2:2545353/148‑1 (MQ=255)
                                                                                                                                    cAACTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATAGCTATGATGCTCCCGCTTTAGATGGGTTGTCGTGACTTGTCCAACATCTATGTTTTTATCAat               >  2:3211846/1‑148 (MQ=255)
                                                                                                                                     aaCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAata              >  2:1528835/1‑148 (MQ=255)
                                                                                                                                      aCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAatat             >  2:3283941/1‑148 (MQ=255)
                                                                                                                                      aCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAatat             <  2:2684772/148‑1 (MQ=255)
                                                                                                                                      aCTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAatat             >  2:1762424/1‑148 (MQ=255)
                                                                                                                                       cTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtc            <  1:2036740/148‑1 (MQ=255)
                                                                                                                                       cTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtc            >  1:2478436/1‑148 (MQ=255)
                                                                                                                                       cTTTGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtc            >  2:166630/1‑148 (MQ=255)
                                                                                                                                        tttGGATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtct           >  1:1018374/1‑148 (MQ=255)
                                                                                                                                           ggATACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtcttct        <  2:2287910/148‑1 (MQ=255)
                                                                                                                                              tACAGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGa     >  1:2691381/1‑148 (MQ=255)
                                                                                                                                                 cGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGa  >  2:933001/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTGAACCTAATACGATGTTACCGTCATTTTTGTCCATCAGTCGTACACCGACCCCAGTTGCTTCGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGA  >  CP009273/747828‑748120

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