Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,107,728 C→A 13.8% E27D (GAG→GAT wcaN ← putative regulatory subunit for GalU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,107,7280CA13.8% 767.8 / 67.3 245E27D (GAG→GATwcaNputative regulatory subunit for GalU
Reads supporting (aligned to +/- strand):  ref base C (91/120);  new base A (17/17);  total (108/137)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 256 total aligned reads displayed.

GAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAATTA  >  CP009273/2107581‑2107876
                                                                                                                                                   |                                                                                                                                                    
gAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCAtct                                                                                                                                                     >  1:4088390/1‑149 (MQ=255)
 aGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCAtctc                                                                                                                                                    >  2:2287760/1‑149 (MQ=255)
 aGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCAtctc                                                                                                                                                    >  2:582060/1‑149 (MQ=255)
  gTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATctct                                                                                                                                                   <  2:3753000/149‑1 (MQ=255)
      ttCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTgg                                                                                                                                               >  1:3949948/1‑149 (MQ=255)
           aCTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATc                                                                                                                                          <  1:537494/149‑1 (MQ=255)
           aCTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATc                                                                                                                                          >  2:2573346/1‑149 (MQ=255)
               aTAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCt                                                                                                                                      <  2:3592806/149‑1 (MQ=255)
                tAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCtt                                                                                                                                     <  1:73851/149‑1 (MQ=255)
                  agagGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAg                                                                                                                                   >  1:291560/1‑149 (MQ=255)
                    agGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTg                                                                                                                                 >  2:1825496/1‑149 (MQ=255)
                      gtgtCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTggc                                                                                                                               >  1:2275334/1‑149 (MQ=255)
                       tgtCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTggca                                                                                                                              <  2:2845792/149‑1 (MQ=255)
                       tgtCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTggca                                                                                                                              <  1:4223/149‑1 (MQ=255)
                        gtCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTggcag                                                                                                                             >  2:3122503/1‑149 (MQ=255)
                        gtCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTggcag                                                                                                                             <  1:151508/149‑1 (MQ=255)
                        gtCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTggcag                                                                                                                             <  1:151444/149‑1 (MQ=255)
                         tCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTggcagg                                                                                                                            <  2:1762897/149‑1 (MQ=255)
                               tGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACat                                                                                                                      <  2:2864551/149‑1 (MQ=255)
                                     tCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCAt                                                                                                                <  1:1980613/149‑1 (MQ=255)
                                       gACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATcc                                                                                                              <  1:1639007/149‑1 (MQ=255)
                                       gACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATcc                                                                                                              <  2:1549916/149‑1 (MQ=255)
                                          cgcgTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCga                                                                                                           <  2:913513/149‑1 (MQ=255)
                                                 ttGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGc                                                                                                    >  2:856041/1‑149 (MQ=255)
                                                 ttGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGc                                                                                                    <  1:2818317/149‑1 (MQ=255)
                                                    gACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTAc                                                                                                 <  2:3886016/149‑1 (MQ=255)
                                                      cgcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAg                                                                                               <  1:997343/149‑1 (MQ=255)
                                                       gcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAgg                                                                                              >  2:941288/1‑149 (MQ=255)
                                                         gTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGaa                                                                                            <  2:2286735/149‑1 (MQ=255)
                                                         gTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGaa                                                                                            <  2:3058955/149‑1 (MQ=255)
                                                             gTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAAc                                                                                        >  2:2669392/1‑149 (MQ=255)
                                                              ttACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACt                                                                                       >  1:1901917/1‑149 (MQ=255)
                                                              ttACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACt                                                                                       <  1:2690053/149‑1 (MQ=255)
                                                                   aggaggATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCttt                                                                                  >  1:3434767/1‑149 (MQ=255)
                                                                    ggaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCtttt                                                                                 >  1:3520068/1‑149 (MQ=255)
                                                                     gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTa                                                                                >  1:3801428/1‑149 (MQ=255)
                                                                      aggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTaa                                                                               >  1:236781/1‑149 (MQ=255)
                                                                      aggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTaa                                                                               >  2:112252/1‑149 (MQ=255)
                                                                             cTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTc                                                                        <  1:2138248/149‑1 (MQ=255)
                                                                             cTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTc                                                                        >  1:3287075/1‑149 (MQ=255)
                                                                              tttGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCa                                                                       >  1:1297697/1‑149 (MQ=255)
                                                                               ttGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCat                                                                      <  1:3918519/149‑1 (MQ=255)
                                                                               ttGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCat                                                                      >  1:1444128/1‑149 (MQ=255)
                                                                                 gATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCatta                                                                    <  2:2115964/149‑1 (MQ=255)
                                                                                  aTCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCattat                                                                   >  1:3853737/1‑149 (MQ=255)
                                                                                   tCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCattatt                                                                  <  2:69299/149‑1 (MQ=255)
                                                                                   tCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCattatt                                                                  <  1:3137860/149‑1 (MQ=255)
                                                                                      cTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCa                                                               <  2:3517338/149‑1 (MQ=255)
                                                                                      cTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCa                                                               <  1:1452124/149‑1 (MQ=255)
                                                                                        gCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATc                                                             <  2:3109484/149‑1 (MQ=255)
                                                                                        gCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATc                                                             <  2:1992126/149‑1 (MQ=255)
                                                                                        gCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATc                                                             <  2:2241958/149‑1 (MQ=255)
                                                                                        gCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATc                                                             <  2:3948224/149‑1 (MQ=255)
                                                                                          aGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCa                                                           >  1:3080780/1‑149 (MQ=255)
                                                                                           gCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCAc                                                          >  1:1829830/1‑149 (MQ=255)
                                                                                              acaATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTg                                                       >  1:1904675/1‑149 (MQ=255)
                                                                                              acaATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTg                                                       >  1:3626017/1‑149 (MQ=255)
                                                                                                aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTa                                                     >  2:2993258/1‑149 (MQ=255)
                                                                                                 aTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTaa                                                    >  1:8908/1‑149 (MQ=255)
                                                                                                 aTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTaa                                                    >  1:1072108/1‑149 (MQ=255)
                                                                                                 aTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTaa                                                    >  1:486328/1‑149 (MQ=255)
                                                                                                   ctcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTaaaa                                                  <  1:3543803/149‑1 (MQ=255)
                                                                                                    tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAAt                                                 <  1:1839076/149‑1 (MQ=255)
                                                                                                       tCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGt                                              >  2:1031855/1‑149 (MQ=255)
                                                                                                            aaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTg                                         >  2:2569230/1‑149 (MQ=255)
                                                                                                            aaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTg                                         >  2:2593709/1‑149 (MQ=255)
                                                                                                             aTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGa                                        <  2:142758/149‑1 (MQ=255)
                                                                                                             aTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGa                                        <  2:546040/149‑1 (MQ=255)
                                                                                                              tGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGaa                                       >  1:487074/1‑149 (MQ=255)
                                                                                                               gTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAAt                                      <  2:1501191/149‑1 (MQ=255)
                                                                                                               gTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAAt                                      <  2:1208542/149‑1 (MQ=255)
                                                                                                               gTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAAt                                      <  1:3724590/149‑1 (MQ=255)
                                                                                                               gTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAAt                                      <  1:3726697/149‑1 (MQ=255)
                                                                                                                tACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAAtt                                     <  1:1848321/149‑1 (MQ=255)
                                                                                                                 aCTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATta                                    <  2:1225911/149‑1 (MQ=255)
                                                                                                                  cTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATtat                                   <  1:2572089/149‑1 (MQ=255)
                                                                                                                   tGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATtata                                  >  1:645393/1‑149 (MQ=255)
                                                                                                                       tCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCtt                              <  1:2923318/149‑1 (MQ=255)
                                                                                                                        cATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTg                             >  1:1220510/1‑149 (MQ=255)
                                                                                                                         aTTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGt                            >  2:2383496/1‑149 (MQ=255)
                                                                                                                            ggCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCg                         >  2:3668443/1‑149 (MQ=255)
                                                                                                                             gCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGa                        <  2:3853737/149‑1 (MQ=255)
                                                                                                                              cTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGAt                       <  1:333826/149‑1 (MQ=255)
                                                                                                                                 gTCGACGATTGGTAGCATATCTTTGGGTATCGACTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGAtttt                    >  2:3921214/1‑149 (MQ=255)
                                                                                                                                   cgacgaTTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGAtttttt                  <  2:768187/149‑1 (MQ=255)
                                                                                                                                   cgacgaTTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGAtttttt                  <  1:1740587/149‑1 (MQ=255)
                                                                                                                                   cgacgaTTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGAtttttt                  <  2:2617130/149‑1 (MQ=255)
                                                                                                                                     acgaTTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTc                >  1:169856/1‑149 (MQ=255)
                                                                                                                                      cgaTTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCg               <  2:3390369/149‑1 (MQ=255)
                                                                                                                                       gaTTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGc              <  2:407278/149‑1 (MQ=255)
                                                                                                                                       gaTTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGc              <  2:1960630/149‑1 (MQ=255)
                                                                                                                                       gaTTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGc              <  2:3996980/149‑1 (MQ=255)
                                                                                                                                         ttGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATCTGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCa            <  2:2701831/149‑1 (MQ=255)
                                                                                                                                           ggTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGc          <  1:3017270/149‑1 (MQ=255)
                                                                                                                                            gTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCa         <  1:1957858/149‑1 (MQ=255)
                                                                                                                                            gTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCa         <  2:896822/149‑1 (MQ=255)
                                                                                                                                                cATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTAGCCAGCATCaa     >  2:2672681/1‑149 (MQ=255)
                                                                                                                                                  tctcTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAAtt   >  1:1010733/1‑149 (MQ=255)
                                                                                                                                                   ctctTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAATTa  <  1:3122503/149‑1 (MQ=255)
                                                                                                                                                   ctctTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAATTa  <  1:2669392/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAATTA  >  CP009273/2107581‑2107876

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: