Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 918,501 G→A 6.8% D45N (GAC→AAC)  clpS → regulatory protein for ClpA substrate specificity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273918,5010GA 6.8% 664.9 / 25.2 191D45N (GAC→AAC) clpSregulatory protein for ClpA substrate specificity
Reads supporting (aligned to +/- strand):  ref base G (81/97);  new base A (7/6);  total (88/103)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 193 total aligned reads displayed.

GAAGATGATAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTACGCGAGGGA  >  CP009273/918352‑918649
                                                                                                                                                     |                                                                                                                                                    
gAAGATGATAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTg                                                                                                                                                      >  2:362514/1‑150 (MQ=255)
gAAGATGATAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTg                                                                                                                                                      >  2:1247127/1‑150 (MQ=255)
 aaGATGATAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGa                                                                                                                                                     <  2:4566752/150‑1 (MQ=255)
  aGATGATAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGAc                                                                                                                                                    >  2:343141/1‑150 (MQ=255)
    atgatACCTGACAATGGTTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACgt                                                                                                                                                  <  1:711253/150‑1 (MQ=255)
    atgatAACTCACAATGGGAAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACgt                                                                                                                                                  <  1:711552/150‑1 (MQ=255)
     tgatAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACgtg                                                                                                                                                 >  2:3314690/1‑150 (MQ=255)
        tAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTa                                                                                                                                              >  2:2343228/1‑150 (MQ=255)
               cAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAt                                                                                                                                       >  1:1923351/1‑150 (MQ=255)
                  tGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttct                                                                                                                                    >  1:1777184/1‑150 (MQ=255)
                   gggTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttctt                                                                                                                                   >  2:3932333/1‑150 (MQ=255)
                   gggTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttctt                                                                                                                                   <  2:4130900/150‑1 (MQ=255)
                   gggTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttctt                                                                                                                                   <  2:2721882/150‑1 (MQ=255)
                     gTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttctttt                                                                                                                                 >  2:2396149/1‑150 (MQ=255)
                     gTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttctttt                                                                                                                                 >  2:2399555/1‑150 (MQ=255)
                      tAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttcttttc                                                                                                                                >  2:2270788/1‑150 (MQ=255)
                      tAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttcttttc                                                                                                                                >  2:4344500/1‑150 (MQ=255)
                        aaaCGCACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAAttcttttctt                                                                                                                              <  2:1188302/150‑1 (MQ=255)
                          acgaacgaCTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTat                                                                                                                            <  1:1393515/150‑1 (MQ=255)
                             aacgaCTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatgat                                                                                                                         <  1:2044971/150‑1 (MQ=255)
                             aacgaCTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTatgat                                                                                                                         <  1:2044961/150‑1 (MQ=255)
                               cgaCTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGt                                                                                                                       <  1:2173986/150‑1 (MQ=255)
                                      ctggACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGt                                                                                                                <  2:2010065/150‑1 (MQ=255)
                                        ggACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGc                                                                                                              >  2:316980/1‑150 (MQ=255)
                                        ggACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGc                                                                                                              >  2:514081/1‑150 (MQ=255)
                                        ggACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGAGGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGc                                                                                                              >  2:4707748/1‑150 (MQ=255)
                                           cTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAAc                                                                                                           >  1:1763254/1‑150 (MQ=255)
                                            tttGATCAACTTGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACg                                                                                                          <  2:961892/150‑1 (MQ=255)
                                            tttGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACg                                                                                                          <  2:961963/150‑1 (MQ=255)
                                            tttGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACg                                                                                                          <  2:1583674/150‑1 (MQ=255)
                                            tttGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACg                                                                                                          <  2:4083079/150‑1 (MQ=255)
                                              tGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCa                                                                                                        >  1:487059/1‑150 (MQ=255)
                                                aTCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAAt                                                                                                      <  2:3807841/150‑1 (MQ=255)
                                                 tCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAAtt                                                                                                     <  1:1193944/150‑1 (MQ=255)
                                                 tCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAAtt                                                                                                     <  1:3036902/150‑1 (MQ=255)
                                                   aaCTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAAACGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGTTGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGa                                                                                                   <  2:3016905/150‑1 (MQ=255)
                                                    aCTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGAt                                                                                                  <  1:2734681/150‑1 (MQ=255)
                                                     cTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATg                                                                                                 <  2:4080209/150‑1 (MQ=255)
                                                      tGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGc                                                                                                >  1:1780608/1‑150 (MQ=255)
                                                       ggCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCt                                                                                               >  1:4730921/1‑150 (MQ=255)
                                                        gCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTc                                                                                              >  1:92688/1‑150 (MQ=255)
                                                           gaagaaAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCt                                                                                           <  1:3603240/150‑1 (MQ=255)
                                                               aaaaaGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTc                                                                                       <  1:2743061/150‑1 (MQ=255)
                                                                 aaaGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTTTTCAc                                                                                     <  2:3845951/150‑1 (MQ=255)
                                                                 aaaGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCAc                                                                                     <  2:3761775/150‑1 (MQ=255)
                                                                  aaGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACt                                                                                    <  1:2754276/150‑1 (MQ=255)
                                                                   aGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTa                                                                                   >  2:2131515/1‑150 (MQ=255)
                                                                   aGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTa                                                                                   >  1:1472541/1‑150 (MQ=255)
                                                                    gTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTAcc                                                                                 >  2:4322969/1‑150 (MQ=255)
                                                                         cgACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAggg                                                                             <  1:3961533/150‑1 (MQ=255)
                                                                             gcgcTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAg                                                                         >  1:2432508/1‑150 (MQ=255)
                                                                              cgcTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAgg                                                                        >  1:3285973/1‑150 (MQ=255)
                                                                              cgcTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAgg                                                                        <  1:2386833/150‑1 (MQ=255)
                                                                                 tAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCa                                                                     >  2:2776922/1‑150 (MQ=255)
                                                                                 tAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCa                                                                     <  1:3143825/150‑1 (MQ=255)
                                                                                 tAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCa                                                                     <  1:3546160/150‑1 (MQ=255)
                                                                                 tAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCa                                                                     >  1:2895806/1‑150 (MQ=255)
                                                                                   aaaCCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCAtt                                                                   >  1:1099234/1‑150 (MQ=255)
                                                                                      ccGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGc                                                                <  1:3350730/150‑1 (MQ=255)
                                                                                       cGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCg                                                               <  1:2313461/150‑1 (MQ=255)
                                                                                           aTCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGt                                                           <  2:2245662/150‑1 (MQ=255)
                                                                                           aTCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGt                                                           <  2:1880151/150‑1 (MQ=255)
                                                                                            tCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTc                                                          <  1:1958168/150‑1 (MQ=255)
                                                                                             cTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCt                                                         >  1:3633560/1‑150 (MQ=255)
                                                                                              tATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCtt                                                        <  1:2572155/150‑1 (MQ=255)
                                                                                               aTGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCttt                                                       <  2:1766592/150‑1 (MQ=255)
                                                                                               aTGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCttt                                                       >  1:37903/1‑150 (MQ=255)
                                                                                                 gTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCACTTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTAc                                                     >  1:4454771/1‑150 (MQ=255)
                                                                                                       aGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGa                                               >  1:3158182/1‑150 (MQ=255)
                                                                                                           atatTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGtt                                           >  1:810660/1‑150 (MQ=255)
                                                                                                            tatTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTg                                          <  1:3484475/150‑1 (MQ=255)
                                                                                                               tAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAg                                       >  2:1159050/1‑150 (MQ=255)
                                                                                                               gaGTCAATGATGATTACACTCCGATGTAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAg                                       <  1:169533/149‑1 (MQ=255)
                                                                                                                 gTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGaa                                     >  1:4234941/1‑150 (MQ=255)
                                                                                                                    aaTGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAAcc                                  <  1:4284690/150‑1 (MQ=255)
                                                                                                                     atgatgATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCa                                 >  2:1492671/1‑150 (MQ=255)
                                                                                                                     atgatgATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCa                                 >  2:3882537/1‑150 (MQ=255)
                                                                                                                      tgatgaTTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCaa                                >  2:2107869/1‑150 (MQ=255)
                                                                                                                        atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAg                              >  1:914854/1‑150 (MQ=255)
                                                                                                                        atgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTGTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGCGGTTGCAGAAACCAAAg                              >  1:901832/1‑150 (MQ=255)
                                                                                                                        atgatTACACTCCGATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAg                              <  1:3774903/150‑1 (MQ=255)
                                                                                                                         tgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGt                             <  1:3529893/150‑1 (MQ=255)
                                                                                                                         tgatTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGt                             >  2:1212699/1‑150 (MQ=255)
                                                                                                                           atTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTgg                           >  2:3555405/1‑150 (MQ=255)
                                                                                                                                acTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATg                      >  1:1216702/1‑150 (MQ=255)
                                                                                                                                  tCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGt                    <  1:514081/150‑1 (MQ=255)
                                                                                                                                  tCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGt                    <  1:3055222/150‑1 (MQ=255)
                                                                                                                                    cGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGa                  >  1:112788/1‑150 (MQ=255)
                                                                                                                                     gATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGaa                 >  2:665316/1‑150 (MQ=255)
                                                                                                                                     gATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGaa                 >  1:403587/1‑150 (MQ=255)
                                                                                                                                     gATGGAGTTTGTTATTAACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGaa                 >  1:2779041/1‑150 (MQ=255)
                                                                                                                                       tGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACa               >  1:1007579/1‑150 (MQ=255)
                                                                                                                                       tGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACa               <  2:1651534/150‑1 (MQ=255)
                                                                                                                                       tGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACa               >  1:3391161/1‑150 (MQ=255)
                                                                                                                                         gAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAg             <  1:4344430/150‑1 (MQ=255)
                                                                                                                                           gTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTa           <  1:1247127/150‑1 (MQ=255)
                                                                                                                                              tGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTAcgc        >  1:340028/1‑150 (MQ=255)
                                                                                                                                               gTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTAcgcg       <  2:2926026/150‑1 (MQ=255)
                                                                                                                                                 tATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTACGCGAg     <  2:3158182/150‑1 (MQ=255)
                                                                                                                                                 tATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTACGCGAg     <  2:92688/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAAGATGATAACTGACAATGGGTAAAACGAACGACTGGCTGGACTTTGATCAACTGGCGGAAGAAAAAGTTCGCGACGCGCTAAAACCGCCATCTATGTATAAAGTGATATTAGTCAATGATGATTACACTCCGATGGAGTTTGTTATTGACGTGTTACAAAAATTCTTTTCTTATGATGTAGAACGTGCAACGCAATTGATGCTCGCTGTTCACTACCAGGGGAAGGCCATTTGCGGAGTCTTTACCGCCGAGGTTGCAGAAACCAAAGTGGCGATGGTGAACAAGTACGCGAGGGA  >  CP009273/918352‑918649

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: