Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 888,964 C→T 100% intergenic (+75/‑276) ybjN → / → potF negative regulator of motility; multicopy suppressor of coaA(Ts)/putrescine transporter subunit: periplasmic‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273888,9640CT97.4% 757.6 / 7.5 233intergenic (+75/‑276)ybjN/potFnegative regulator of motility; multicopy suppressor of coaA(Ts)/putrescine transporter subunit: periplasmic‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base T (99/127);  total (102/131)
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.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 238 total aligned reads displayed.

CCATCAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAG  >  CP009273/888820‑889108
                                                                                                                                                |                                                                                                                                                
ccATCAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTa                                                                                                                                               <  1:2198863/148‑1 (MQ=255)
ccATCAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTa                                                                                                                                               <  2:628335/148‑1 (MQ=255)
  aTCAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTaaa                                                                                                                                             >  2:1790373/1‑148 (MQ=255)
    cAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGa                                                                                                                                           <  1:1763196/148‑1 (MQ=255)
    cAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGa                                                                                                                                           <  1:1298052/148‑1 (MQ=255)
    cAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGa                                                                                                                                           <  2:448748/148‑1 (MQ=255)
     aaCTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGat                                                                                                                                          <  2:3528306/148‑1 (MQ=255)
      aCTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGatt                                                                                                                                         <  2:2792336/148‑1 (MQ=255)
      aCTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGatt                                                                                                                                         <  2:813275/148‑1 (MQ=255)
       ctgctgTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGatta                                                                                                                                        <  1:1631009/148‑1 (MQ=255)
       ctgctgTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGatta                                                                                                                                        <  1:3061231/148‑1 (MQ=255)
          ctgTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTAttt                                                                                                                                     <  1:1790373/148‑1 (MQ=255)
           tgTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTc                                                                                                                                    >  2:1509540/1‑148 (MQ=255)
           tgTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTc                                                                                                                                    <  1:502903/148‑1 (MQ=255)
            gTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCa                                                                                                                                   >  1:2979980/1‑148 (MQ=255)
                  gACTGATGATGAAGGGCAACACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGGGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTtctc                                                                                                                             <  2:3088221/148‑1 (MQ=255)
                  gACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTtctc                                                                                                                             <  2:634802/148‑1 (MQ=255)
                    cTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCtt                                                                                                                           <  1:125229/148‑1 (MQ=255)
                     tgatgatgaAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTtg                                                                                                                          <  2:1769704/148‑1 (MQ=255)
                     ggatgatgaAGGGCAAAACAACGTTACCGAAAACTCTTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTtg                                                                                                                          <  2:2240490/147‑1 (MQ=255)
                        tgatgaAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTGTCTTgtgt                                                                                                                       >  1:3375105/1‑148 (MQ=255)
                        tgatgaAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTgtgt                                                                                                                       >  1:3375101/1‑148 (MQ=255)
                             aaGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGAtt                                                                                                                  >  2:772026/1‑148 (MQ=255)
                             aaGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGAtt                                                                                                                  <  2:460564/148‑1 (MQ=255)
                              aGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGAttt                                                                                                                 <  1:2111208/148‑1 (MQ=255)
                                  cAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCt                                                                                                             <  1:1384151/148‑1 (MQ=255)
                                     aacaacGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTa                                                                                                          >  1:1260177/1‑148 (MQ=255)
                                      acaacGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTAt                                                                                                         <  1:1936407/148‑1 (MQ=255)
                                       caacGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATc                                                                                                        >  2:668998/1‑148 (MQ=255)
                                       caacGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATc                                                                                                        >  2:2494913/1‑148 (MQ=255)
                                          cGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATcaca                                                                                                     >  1:2034120/1‑148 (MQ=255)
                                          cGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATcaca                                                                                                     <  1:1890789/148‑1 (MQ=255)
                                          cGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATcaca                                                                                                     <  1:904769/148‑1 (MQ=255)
                                             tACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATag                                                                                                  <  1:88737/148‑1 (MQ=255)
                                             tACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATag                                                                                                  >  2:3129093/1‑148 (MQ=255)
                                             tACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATag                                                                                                  <  2:2157429/148‑1 (MQ=255)
                                                 gAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCa                                                                                              <  1:1689764/148‑1 (MQ=255)
                                                  aaaaCTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCaa                                                                                             >  2:1579532/1‑148 (MQ=255)
                                                  aaaaCTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCaa                                                                                             >  2:2855841/1‑148 (MQ=255)
                                                   aaaCTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCaaa                                                                                            <  1:2829057/148‑1 (MQ=255)
                                                   aaaCTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCaaa                                                                                            >  2:2622957/1‑148 (MQ=255)
                                                     aCTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCACAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCCAata                                                                                          >  2:2572144/1‑148 (MQ=255)
                                                        aTTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGc                                                                                       <  2:2282044/148‑1 (MQ=255)
                                                        aTTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGc                                                                                       >  2:2999378/1‑148 (MQ=255)
                                                          ttCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCAt                                                                                     >  2:2252018/1‑148 (MQ=255)
                                                          ttCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCAt                                                                                     >  1:1003267/1‑148 (MQ=255)
                                                            cctccACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATaa                                                                                   <  2:1003267/148‑1 (MQ=255)
                                                              tccACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATaaaa                                                                                 <  1:2305409/148‑1 (MQ=255)
                                                                    gATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTGTCTTGGGTCGATTTGGCTTTATCACATATAGCAAATATGCATAAAAATTTGt                                                                           >  1:2533957/1‑148 (MQ=255)
                                                                    gATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGt                                                                           <  1:2862221/148‑1 (MQ=255)
                                                                     aTAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGtt                                                                          <  1:239033/148‑1 (MQ=255)
                                                                     aTAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGtt                                                                          >  1:3554721/1‑148 (MQ=255)
                                                                         cTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAGTTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAAt                                                                      <  1:1408755/148‑1 (MQ=255)
                                                                             tttCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCg                                                                  <  1:1809253/148‑1 (MQ=255)
                                                                             tttCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCGTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCg                                                                  >  1:2456254/1‑148 (MQ=255)
                                                                              ttCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAACTTTGTTAAATACCGt                                                                 >  2:1427977/1‑148 (MQ=255)
                                                                              ttCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGt                                                                 <  2:897738/148‑1 (MQ=255)
                                                                              ttCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGt                                                                 <  2:2979980/148‑1 (MQ=255)
                                                                              ttCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGt                                                                 <  1:1165973/148‑1 (MQ=255)
                                                                               tCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGtt                                                                <  1:1753975/148‑1 (MQ=255)
                                                                               tCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATTCATAAAAATTTGTTAAATACCttt                                                                >  1:992845/1‑148 (MQ=255)
                                                                                    cACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTTTCACATAGAGCAAATATGCATTAAAATTTGTTAAATACCGTTTTTTa                                                           >  2:1321686/1‑148 (MQ=255)
                                                                                    cACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTa                                                           >  2:1228112/1‑148 (MQ=255)
                                                                                      cGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAAt                                                         <  2:1658192/148‑1 (MQ=255)
                                                                                       gCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATc                                                        <  1:1083570/148‑1 (MQ=255)
                                                                                       gCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATc                                                        <  2:1985733/148‑1 (MQ=255)
                                                                                               tGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCtat                                                <  2:1527248/148‑1 (MQ=255)
                                                                                               tGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCtat                                                <  2:3274933/148‑1 (MQ=255)
                                                                                                 gTGCAGTGGCTGCGCGCTTCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAg                                              <  2:2725516/148‑1 (MQ=255)
                                                                                                      gTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCa                                         >  2:3496627/1‑148 (MQ=255)
                                                                                                      gTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCaa                                        <  2:520038/148‑1 (MQ=255)
                                                                                                       tGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCCAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCaa                                        >  2:721315/1‑148 (MQ=255)
                                                                                                        ggCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCaaa                                       >  1:3044937/1‑148 (MQ=255)
                                                                                                        ggCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCaaa                                       >  1:3045189/1‑148 (MQ=255)
                                                                                                          cTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAAcc                                     >  1:165695/1‑148 (MQ=255)
                                                                                                          cTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAAcc                                     <  1:316320/148‑1 (MQ=255)
                                                                                                               cgcTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGc                                <  1:1727324/148‑1 (MQ=255)
                                                                                                                gcTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCt                               <  1:3181825/148‑1 (MQ=255)
                                                                                                                     aaaaTTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACACTGGCTAAAGt                          >  2:3430708/1‑148 (MQ=255)
                                                                                                                      aaaTTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAATCCCTGGCTAAAGtt                         >  2:3135757/1‑148 (MQ=255)
                                                                                                                       aaTTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTCAAGTTa                        >  2:2207633/1‑148 (MQ=255)
                                                                                                                       aaTTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTa                        >  2:1921740/1‑148 (MQ=255)
                                                                                                                       aaTTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTa                        <  1:816187/148‑1 (MQ=255)
                                                                                                                       aaTTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCAGAGCTATAGTCTCAAACCCTGGCTAAAGTTa                        >  1:2774669/1‑148 (MQ=255)
                                                                                                                        aTTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTAt                       >  1:3281508/1‑148 (MQ=255)
                                                                                                                        aTTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTAt                       <  1:2999813/148‑1 (MQ=255)
                                                                                                                         ttATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTAtt                      >  2:276230/1‑148 (MQ=255)
                                                                                                                         ttATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTAtt                      >  2:1234618/1‑148 (MQ=255)
                                                                                                                          tATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTc                     >  2:1726527/1‑148 (MQ=255)
                                                                                                                              tgctgTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGc                 >  1:2463142/1‑148 (MQ=255)
                                                                                                                              tgctgTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGc                 >  2:1558980/1‑148 (MQ=255)
                                                                                                                                   tttttAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGc            >  1:663694/1‑148 (MQ=255)
                                                                                                                                    ttttAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCt           >  1:1994708/1‑148 (MQ=255)
                                                                                                                                     tttAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCACATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCtt          >  2:1128166/1‑148 (MQ=255)
                                                                                                                                        aaCCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTta       <  1:69672/148‑1 (MQ=255)
                                                                                                                                           cTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTtatat    >  2:2548780/1‑148 (MQ=255)
                                                                                                                                           cTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTtatat    >  2:3349961/1‑148 (MQ=255)
                                                                                                                                           cTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTtatat    >  2:2868101/1‑148 (MQ=255)
                                                                                                                                            ttttttTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTtatata   <  2:357153/148‑1 (MQ=255)
                                                                                                                                             ttttttAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAg  >  1:406398/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CCATCAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAG  >  CP009273/888820‑889108

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: