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.5% 800.8 / 6.3 244intergenic (+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 (4/2);  new base T (127/110);  total (132/112)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 246 total aligned reads displayed.

GCCCATCAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCG  >  CP009273/888818‑889110
                                                                                                                                                  |                                                                                                                                                  
gCCCATCAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTa                                                                                                                                                 <  2:4605395/150‑1 (MQ=255)
 cccATCAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTaa                                                                                                                                                >  2:4485289/1‑150 (MQ=255)
    aTCAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGa                                                                                                                                             >  1:1081371/1‑150 (MQ=255)
    aTCAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGa                                                                                                                                             >  2:1412455/1‑150 (MQ=255)
    aTCAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGa                                                                                                                                             >  2:2486151/1‑150 (MQ=255)
      cAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGatt                                                                                                                                           >  1:3077267/1‑150 (MQ=255)
          tgctgTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTAttt                                                                                                                                       >  1:250116/1‑150 (MQ=255)
            ctgTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCa                                                                                                                                     <  1:4879962/150‑1 (MQ=255)
              gTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACt                                                                                                                                   <  2:160521/150‑1 (MQ=255)
               ttACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTCTCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACtt                                                                                                                                  <  1:3940419/150‑1 (MQ=255)
               ttACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACtt                                                                                                                                  <  1:46443/150‑1 (MQ=255)
                tACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTtc                                                                                                                                 >  2:4528349/1‑150 (MQ=255)
                 aCCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTtct                                                                                                                                <  1:3483003/150‑1 (MQ=255)
                     aCTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTtg                                                                                                                            <  1:1031855/150‑1 (MQ=255)
                     aCTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTtg                                                                                                                            <  1:2148433/150‑1 (MQ=255)
                     aCTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTtg                                                                                                                            <  1:1459273/150‑1 (MQ=255)
                      cTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTtgt                                                                                                                           <  2:3245944/150‑1 (MQ=255)
                      cTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTtgt                                                                                                                           >  2:466877/1‑150 (MQ=255)
                        gatgatgaAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTgtgt                                                                                                                         >  1:2973929/1‑150 (MQ=255)
                             tgaAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGAtt                                                                                                                    >  1:1719972/1‑150 (MQ=255)
                               aaGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTg                                                                                                                  >  2:2279720/1‑150 (MQ=255)
                               aaGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTg                                                                                                                  >  2:3863448/1‑150 (MQ=255)
                               aaGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTg                                                                                                                  <  1:3620462/150‑1 (MQ=255)
                                aGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTgg                                                                                                                 <  2:4957674/150‑1 (MQ=255)
                                aGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTgg                                                                                                                 >  2:4598910/1‑150 (MQ=255)
                                   gCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCtt                                                                                                              <  1:1753257/150‑1 (MQ=255)
                                    cAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCttt                                                                                                             >  2:4052411/1‑150 (MQ=255)
                                     aaaaCAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTa                                                                                                            >  1:2802475/1‑150 (MQ=255)
                                        acaacGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATca                                                                                                         >  1:2237964/1‑150 (MQ=255)
                                          aacGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATcaca                                                                                                       >  1:183545/1‑150 (MQ=255)
                                           acGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACCTCTCTTGTGTCGATTTGGCTTTATCACAt                                                                                                      <  2:3174548/150‑1 (MQ=255)
                                            cGTTACCGACAACTATTTCCTCCACTGATCACTCATTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCATAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATa                                                                                                     <  2:3443928/150‑1 (MQ=255)
                                             gTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATag                                                                                                    >  2:1310695/1‑150 (MQ=255)
                                               tACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATagag                                                                                                  <  1:4782873/150‑1 (MQ=255)
                                               tACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATagag                                                                                                  >  2:4854406/1‑150 (MQ=255)
                                                aCCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAACACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGc                                                                                                 <  2:2311001/150‑1 (MQ=255)
                                                  cGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCaa                                                                                               <  1:218274/150‑1 (MQ=255)
                                                       aCTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATg                                                                                          <  2:57934/150‑1 (MQ=255)
                                                        cTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGc                                                                                         >  2:476826/1‑150 (MQ=255)
                                                           tttCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATa                                                                                      <  2:3830078/150‑1 (MQ=255)
                                                             tcctccACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATaaa                                                                                    <  1:3424859/150‑1 (MQ=255)
                                                             tcctccACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATaaa                                                                                    >  2:191329/1‑150 (MQ=255)
                                                              cctccACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTGTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATaaaa                                                                                   >  2:1061946/1‑150 (MQ=255)
                                                                tccACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAAt                                                                                 >  2:519296/1‑150 (MQ=255)
                                                                 ccACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAAtt                                                                                <  1:4081343/150‑1 (MQ=255)
                                                                  cACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCACAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAAttt                                                                               >  2:964978/1‑150 (MQ=255)
                                                                      gATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTa                                                                           >  1:2264779/1‑150 (MQ=255)
                                                                       aTAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTaa                                                                          >  2:4430389/1‑150 (MQ=255)
                                                                       aTAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTaa                                                                          <  2:3665708/150‑1 (MQ=255)
                                                                        tAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTaaa                                                                         >  1:4503607/1‑150 (MQ=255)
                                                                        tAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTaaa                                                                         >  2:2151139/1‑150 (MQ=255)
                                                                           cTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGAGTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATAc                                                                      >  1:3563690/1‑150 (MQ=255)
                                                                             ccTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCg                                                                    <  2:420365/150‑1 (MQ=255)
                                                                              cTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGt                                                                   <  2:659203/150‑1 (MQ=255)
                                                                               tttCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGtt                                                                  <  2:3119135/150‑1 (MQ=255)
                                                                                ttCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGttt                                                                 >  1:1911308/1‑150 (MQ=255)
                                                                                 tCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGtttt                                                                <  1:1426844/150‑1 (MQ=255)
                                                                                    aGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCACCTATGCATAAAAATTTGTTAATTACCGTTTTTTa                                                             >  2:4372687/1‑150 (MQ=255)
                                                                                    aGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTa                                                             >  1:4290504/1‑150 (MQ=255)
                                                                                      cACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAAt                                                           >  2:3536784/1‑150 (MQ=255)
                                                                                       aCGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATc                                                          >  1:3258490/1‑150 (MQ=255)
                                                                                        cGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATcc                                                         <  1:1647843/150‑1 (MQ=255)
                                                                                          cAGTCGCTGGTGCAGTGGCTGCGCGCTTCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGa                                                       <  2:1779959/150‑1 (MQ=255)
                                                                                           aGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAAGCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATGTGTTAAATACCGTTTTTTAATCCgtg                                                      >  1:4942285/1‑148 (MQ=255)
                                                                                            gTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGc                                                     >  2:4124968/1‑150 (MQ=255)
                                                                                              cGCTGGTGCAGTGGCTGCTCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCta                                                   <  1:653074/150‑1 (MQ=255)
                                                                                               gCTGGTGCAGGGGCTGCGAGCTGCAAAATTATCTGCTGTTTTTAAACTTTTTTTAAAGATTATTTCACTTCTATTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCtat                                                  <  1:2637895/150‑1 (MQ=255)
                                                                                                cTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCtcta                                                 >  2:4914249/1‑150 (MQ=255)
                                                                                                cTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCtata                                                 >  2:114384/1‑150 (MQ=255)
                                                                                                cTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCtata                                                 <  2:1537881/150‑1 (MQ=255)
                                                                                                  ggTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGt                                               >  1:386002/1‑150 (MQ=255)
                                                                                                  ggTGCAGTGGCTGCGCGCTCCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGt                                               >  1:2711047/1‑150 (MQ=255)
                                                                                                     gCAGTGGCTGCGCGCTGAAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTTGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGtctc                                            <  1:377707/150‑1 (MQ=255)
                                                                                                       aGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCaa                                          <  2:350905/150‑1 (MQ=255)
                                                                                                         tGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAAc                                        >  2:2762601/1‑150 (MQ=255)
                                                                                                           gCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAAccc                                      <  2:2568582/150‑1 (MQ=255)
                                                                                                           gCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAAccc                                      <  1:4588601/150‑1 (MQ=255)
                                                                                                            cTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCt                                     <  2:3379848/150‑1 (MQ=255)
                                                                                                            cTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCt                                     <  1:1530748/150‑1 (MQ=255)
                                                                                                             tGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTg                                    >  1:504169/1‑150 (MQ=255)
                                                                                                               cgcgcTGCAAAATTATCTGCTGTTTTTAACCTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCt                                 >  1:525518/1‑150 (MQ=255)
                                                                                                                  gcTGCAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTaa                               <  2:4697625/150‑1 (MQ=255)
                                                                                                                      cAAAATTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGtt                           >  1:209036/1‑150 (MQ=255)
                                                                                                                        aaaTTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTAt                         >  1:1108949/1‑150 (MQ=255)
                                                                                                                         aaTTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTAtt                        <  2:2496141/150‑1 (MQ=255)
                                                                                                                         aaTTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTAtt                        <  2:1666184/150‑1 (MQ=255)
                                                                                                                          aTTATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATACGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTc                       >  1:3359541/1‑150 (MQ=255)
                                                                                                                           ttATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTca                      <  1:829166/150‑2 (MQ=255)
                                                                                                                           ttATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCt                      <  1:2486151/150‑1 (MQ=255)
                                                                                                                           ttATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCt                      >  2:2379500/1‑150 (MQ=255)
                                                                                                                            tATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCtt                     >  2:2868923/1‑150 (MQ=255)
                                                                                                                            tATCTGCTGTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCtt                     <  2:771148/150‑1 (MQ=255)
                                                                                                                               ctgctgTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCg                  <  2:4661082/150‑1 (MQ=255)
                                                                                                                                 gctgTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGAt                >  2:400003/1‑150 (MQ=255)
                                                                                                                                    gTTTTTAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCt             <  2:2477777/150‑1 (MQ=255)
                                                                                                                                       tttAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCtttt          <  2:2973929/150‑1 (MQ=255)
                                                                                                                                       tttAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCtttt          >  2:896467/1‑150 (MQ=255)
                                                                                                                                       tttAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCtatt          >  2:3416701/1‑150 (MQ=255)
                                                                                                                                         tAACCTTTTTTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTtat        >  2:2714173/1‑150 (MQ=255)
                                                                                                                                               ttttttAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCg  <  2:3077267/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCCATCAACTGCTGTTACCGACTGATGATGAAGGGCAAAACAACGTTACCGAAAACTATTTCCTCCACTGATAACTCCTTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCG  >  CP009273/888818‑889110

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: