Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,137,749 C→A 12.4% L335I (CTC→ATC)  yegE → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,137,7490CA12.4% 1067.8 / 78.7 303L335I (CTC→ATC) yegEputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (139/126);  new base A (19/19);  total (158/145)
Fisher's exact test for biased strand distribution p-value = 8.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 321 total aligned reads displayed.

TGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGCTACTACAACCGCAAT  >  CP009273/2137600‑2137898
                                                                                                                                                     |                                                                                                                                                     
tGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtc                                                                                                                                                       <  1:3482577/150‑1 (MQ=255)
tGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtc                                                                                                                                                       <  1:223503/150‑1 (MQ=255)
 ggTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtct                                                                                                                                                      >  1:309069/1‑150 (MQ=255)
 ggTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtct                                                                                                                                                      <  2:233372/150‑1 (MQ=255)
  gTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtctc                                                                                                                                                     <  1:962880/150‑1 (MQ=255)
  gTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTtctc                                                                                                                                                     <  1:112618/150‑1 (MQ=255)
     aTGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCggg                                                                                                                                                  >  2:2980223/1‑150 (MQ=255)
      tGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGt                                                                                                                                                 <  2:3237845/150‑1 (MQ=255)
       ggaTGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTa                                                                                                                                                >  2:1249367/3‑150 (MQ=255)
        tATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTAc                                                                                                                                               >  2:1980418/1‑150 (MQ=255)
         aTGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACa                                                                                                                                              >  1:3521998/1‑150 (MQ=255)
          tGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGACTATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAg                                                                                                                                             >  2:1824503/1‑150 (MQ=255)
          tGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAg                                                                                                                                             <  1:1091353/150‑1 (MQ=255)
             cTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCa                                                                                                                                          <  2:3149022/150‑1 (MQ=255)
              tttCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACGTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAg                                                                                                                                         >  1:4430465/1‑150 (MQ=255)
                 cGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGaa                                                                                                                                      >  2:4375940/1‑150 (MQ=255)
                     cGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGc                                                                                                                                  >  2:595514/1‑150 (MQ=255)
                      ggAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCt                                                                                                                                 <  2:935596/150‑1 (MQ=255)
                      ggAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCt                                                                                                                                 <  1:3651129/150‑1 (MQ=255)
                       gAACGCAAACACATTTCCGAAAGCGAAACCCACTTTAGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTg                                                                                                                                >  2:2920601/1‑150 (MQ=255)
                         aCGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTgcg                                                                                                                              >  2:2161536/1‑150 (MQ=255)
                         aCGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTgcg                                                                                                                              >  2:2353207/1‑150 (MQ=255)
                          cGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTgcgc                                                                                                                             >  2:360917/1‑150 (MQ=255)
                          cGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTgcgc                                                                                                                             >  2:1436550/1‑150 (MQ=255)
                           gCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGcgcg                                                                                                                            >  1:626934/1‑150 (MQ=255)
                           gCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGAATGGCGTTAGTGGGCACCGAGGGACAAGGGCTGCAAACCACCAAAGCGCTCTGCCAGAATCTCGCGTACAGTCAGGAAGAGCTCcgcg                                                                                                                            >  2:372057/1‑150 (MQ=255)
                            cAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCgg                                                                                                                           >  2:965281/1‑150 (MQ=255)
                            cAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCgg                                                                                                                           <  1:1999781/150‑1 (MQ=255)
                             aaaCACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGa                                                                                                                          >  1:2315773/1‑150 (MQ=255)
                             aaaCACAATTCCGAAAGCGAAAACCACTTTCGGAACGCGAATGAATAATCCGATATCGGTATGACGTTAGTGGGAACAGAGGGACAATGGCAGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGTGTAAAGTCAGCAAGAGCTCCGAGGa                                                                                                                          >  2:1423095/1‑150 (MQ=255)
                               acacATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACt                                                                                                                        <  2:4191714/150‑1 (MQ=255)
                               acacATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACt                                                                                                                        <  1:2591348/150‑1 (MQ=255)
                                  caTTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCaa                                                                                                                     >  1:736652/1‑149 (MQ=255)
                                   aTTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCAcc                                                                                                                    >  1:52193/1‑150 (MQ=255)
                                    tttCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCt                                                                                                                   >  2:3307778/1‑150 (MQ=255)
                                     ttCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCtt                                                                                                                  >  2:3770437/1‑150 (MQ=255)
                                         gAAAGTGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAg                                                                                                              >  2:1072216/1‑150 (MQ=255)
                                          aaaGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGc                                                                                                             <  1:3991557/150‑1 (MQ=255)
                                           aaGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCa                                                                                                            <  2:899431/150‑1 (MQ=255)
                                              cGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACt                                                                                                         <  2:4236142/150‑1 (MQ=255)
                                              cGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACt                                                                                                         >  2:2285613/1‑150 (MQ=255)
                                               gAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTg                                                                                                        >  1:1533074/1‑150 (MQ=255)
                                                aaaCCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGa                                                                                                       >  1:1420223/1‑150 (MQ=255)
                                                 aaCCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGAc                                                                                                      >  2:1725235/1‑150 (MQ=255)
                                                 aaCCCACTTTCGGAACGCGATGGAATATTCCGCTATAGCTATCGCGTAAGTGGGCACCCAGGGACAATTCCAGCAAACCAACAAAGAGCTCAGCCAGTATCTCGGGTACAGTCAGGAAGATCATCGCGGACACACCTTTCAGCAACTGAc                                                                                                      >  1:2804254/1‑150 (MQ=255)
                                                   cccACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCt                                                                                                    >  2:3927246/1‑150 (MQ=255)
                                                    ccACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTg                                                                                                   <  2:3864312/150‑1 (MQ=255)
                                                     cACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTgg                                                                                                  >  2:1682420/1‑150 (MQ=255)
                                                      aCTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGc                                                                                                 >  2:2567956/1‑150 (MQ=255)
                                                      aCTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGc                                                                                                 >  1:4257546/1‑150 (MQ=255)
                                                       cTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGcc                                                                                                >  1:3322164/1‑150 (MQ=255)
                                                       cTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGcc                                                                                                >  1:1000539/1‑150 (MQ=255)
                                                          tCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCgga                                                                                             <  1:3410558/150‑1 (MQ=255)
                                                              aaCGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGAt                                                                                         >  1:3780973/1‑150 (MQ=255)
                                                              aaCGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGAt                                                                                         <  1:2344462/150‑1 (MQ=255)
                                                               aCGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGAtc                                                                                        >  1:1774559/1‑150 (MQ=255)
                                                                 tcgATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGAtctc                                                                                      <  1:2503613/149‑1 (MQ=255)
                                                                 gcgATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGAtctc                                                                                      <  1:2917632/150‑1 (MQ=255)
                                                                 gcgATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGAtctc                                                                                      <  1:896792/150‑1 (MQ=255)
                                                                 gcgATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGAtctc                                                                                      <  1:2438861/150‑1 (MQ=255)
                                                                         atatTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGa                                                                              <  1:4633619/150‑1 (MQ=255)
                                                                         atatTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGa                                                                              <  2:2807907/150‑1 (MQ=255)
                                                                             tCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGAtctc                                                                          >  1:2154003/1‑150 (MQ=255)
                                                                                 cTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAAc                                                                      >  2:4498892/1‑150 (MQ=255)
                                                                                  tATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACa                                                                     <  1:2375467/150‑1 (MQ=255)
                                                                                   aTCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAg                                                                    <  2:2085061/150‑1 (MQ=255)
                                                                                       gTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTg                                                                <  2:3216603/150‑1 (MQ=255)
                                                                                             cGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGc                                                          <  2:1491756/150‑1 (MQ=255)
                                                                                              gTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCt                                                         >  2:4021354/1‑150 (MQ=255)
                                                                                              gTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCt                                                         <  2:176636/150‑1 (MQ=255)
                                                                                                   tGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATaa                                                    <  1:4413693/150‑1 (MQ=255)
                                                                                                       cACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCgg                                                <  2:3485112/150‑1 (MQ=255)
                                                                                                        aCCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGt                                               >  2:1897889/1‑150 (MQ=255)
                                                                                                        aCCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGt                                               <  2:2755396/150‑1 (MQ=255)
                                                                                                         ccGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTg                                              <  1:1436551/150‑1 (MQ=255)
                                                                                                         ccGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTg                                              <  1:360916/150‑1 (MQ=255)
                                                                                                           gAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGaa                                            <  2:1009804/150‑1 (MQ=255)
                                                                                                                 cAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCATCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAac                                      <  2:3389313/150‑1 (MQ=255)
                                                                                                                  aaTGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAaca                                     >  1:117142/1‑150 (MQ=255)
                                                                                                                   aTGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTATCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAacac                                    <  1:1886863/150‑1 (MQ=255)
                                                                                                                    tGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACAcc                                   >  2:677763/1‑150 (MQ=255)
                                                                                                                      gCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTa                                 <  2:3076965/150‑1 (MQ=255)
                                                                                                                      gCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTa                                 <  2:2522956/150‑1 (MQ=255)
                                                                                                                       cTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTAt                                <  1:4319590/150‑1 (MQ=255)
                                                                                                                        tGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTAtt                               <  1:3995273/150‑1 (MQ=255)
                                                                                                                        tGCAAACCAACAAAGCGCTCTGCCAGTTTATCTGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTAtt                               <  1:83667/150‑1 (MQ=255)
                                                                                                                                aacaaAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGaa                       >  2:4614266/1‑150 (MQ=255)
                                                                                                                                aacaaAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGaa                       <  1:1261707/150‑1 (MQ=255)
                                                                                                                                aacaaAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGaa                       >  1:3791819/1‑150 (MQ=255)
                                                                                                                                aacaaAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGaa                       >  1:2769054/1‑150 (MQ=255)
                                                                                                                                  caaAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGaaaa                     <  2:3817150/150‑1 (MQ=255)
                                                                                                                                  caaAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGaaaa                     >  2:91067/1‑150 (MQ=255)
                                                                                                                                  caaAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGaaaa                     >  2:1439202/1‑150 (MQ=255)
                                                                                                                                   aaaGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGaaaaa                    <  1:2920601/150‑1 (MQ=255)
                                                                                                                                     aGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACg                  >  1:266988/1‑150 (MQ=255)
                                                                                                                                       cgcTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCAGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGct                >  2:4059511/1‑150 (MQ=255)
                                                                                                                                          tctGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctact             <  2:3521998/150‑1 (MQ=255)
                                                                                                                                            tGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGCtactac           <  2:4519863/150‑1 (MQ=255)
                                                                                                                                                   ttCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGCTACTACAACCGCa    >  2:286891/1‑150 (MQ=255)
                                                                                                                                                     ctcGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGCTACTACAACCGCAAt  <  2:1470941/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTGATGTATGCCTTTCGTGCGGAACGCAAACACATTTCCGAAAGCGAAACCCACTTTCGGAACGCGATGGAATATTCCGCTATCGGTATGGCGTTAGTGGGCACCGAGGGACAATGGCTGCAAACCAACAAAGCGCTCTGCCAGTTTCTCGGGTACAGTCAGGAAGAGCTGCGCGGACTCACCTTTCAGCAACTGACCTGGCCGGAGGATCTCAATAAAGATCTCCAACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGCTACTACAACCGCAAT  >  CP009273/2137600‑2137898

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: