Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,612,781 C→A 5.6% A399E (GCG→GAG)  yjjJ → putative protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,612,7810CA 5.6% 1273.2 / 30.9 335A399E (GCG→GAG) yjjJputative protein kinase
Reads supporting (aligned to +/- strand):  ref base C (157/159);  new base A (6/13);  total (163/172)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 347 total aligned reads displayed.

GTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAG  >  CP009273/4612632‑4612928
                                                                                                                                                     |                                                                                                                                                   
gTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGc                                                                                                                                                     >  2:87802/1‑150 (MQ=255)
gTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGc                                                                                                                                                     <  2:1543731/150‑1 (MQ=255)
gTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGc                                                                                                                                                     <  2:1668951/150‑1 (MQ=255)
  tttatttaTCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGa                                                                                                                                                   >  2:659836/1‑150 (MQ=255)
  tttatttaTCTGAACAGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGa                                                                                                                                                   >  2:3202877/1‑150 (MQ=255)
     attcatCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATcc                                                                                                                                                >  2:2601441/1‑150 (MQ=255)
        tatCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGc                                                                                                                                             >  2:2683749/1‑150 (MQ=255)
          tCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTg                                                                                                                                           >  2:3996398/1‑150 (MQ=255)
             gAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGgcg                                                                                                                                        >  1:4205219/1‑150 (MQ=255)
               aCCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCgca                                                                                                                                      >  1:2800525/1‑150 (MQ=255)
                ccGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCgcag                                                                                                                                     <  1:762569/150‑1 (MQ=255)
                ccGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCgcag                                                                                                                                     >  2:3872193/1‑150 (MQ=255)
                  gCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCgcagca                                                                                                                                   <  1:1711756/150‑1 (MQ=255)
                  gCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCgcagca                                                                                                                                   >  1:3237013/1‑150 (MQ=255)
                     aTTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAGAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGtt                                                                                                                                >  1:633129/1‑150 (MQ=255)
                      tttGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTc                                                                                                                               >  2:55530/1‑150 (MQ=255)
                      tttGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTc                                                                                                                               <  1:2706812/150‑1 (MQ=255)
                       ttGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCt                                                                                                                              <  2:1852599/150‑1 (MQ=255)
                       ttGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCt                                                                                                                              <  1:502684/150‑1 (MQ=255)
                              gACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTGGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAAc                                                                                                                       >  1:4431174/1‑150 (MQ=255)
                                cGCCCGTCTACGACATGCTGCCGATGGTCTATGCTCCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACgg                                                                                                                     <  1:2911420/150‑1 (MQ=255)
                                cGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACgg                                                                                                                     <  1:3046417/150‑1 (MQ=255)
                                   ccGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCg                                                                                                                  >  1:3525937/1‑150 (MQ=255)
                                     gTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGcc                                                                                                                >  2:1587337/1‑150 (MQ=255)
                                      tCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCa                                                                                                               <  1:4373724/150‑1 (MQ=255)
                                        tACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCaga                                                                                                             <  1:2088161/150‑1 (MQ=255)
                                          cGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAgaga                                                                                                           >  2:2197289/1‑150 (MQ=255)
                                           gACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGAt                                                                                                          >  2:1038441/1‑150 (MQ=255)
                                           gACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGAt                                                                                                          <  2:2955976/150‑1 (MQ=255)
                                            aCATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATc                                                                                                         <  2:2542184/150‑1 (MQ=255)
                                            aCATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATc                                                                                                         <  1:2642759/150‑1 (MQ=255)
                                                      gATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAg                                                                                               >  1:1090995/1‑150 (MQ=255)
                                                        tGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCg                                                                                             <  1:51653/150‑1 (MQ=255)
                                                         ggTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGa                                                                                            <  1:3379445/150‑1 (MQ=255)
                                                         ggTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGa                                                                                            >  2:420985/1‑150 (MQ=255)
                                                              aTGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGt                                                                                       <  2:1923120/150‑1 (MQ=255)
                                                                 cACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTc                                                                                    >  2:2075704/1‑150 (MQ=255)
                                                                 cACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTc                                                                                    >  2:2075182/1‑150 (MQ=255)
                                                                 cACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTc                                                                                    <  2:4175021/150‑1 (MQ=255)
                                                                 cACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTc                                                                                    <  1:3296215/150‑1 (MQ=255)
                                                                    cAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGcc                                                                                 >  2:3215001/1‑150 (MQ=255)
                                                                       aCAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCcaca                                                                              >  2:757400/1‑150 (MQ=255)
                                                                         aGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACAtt                                                                            <  2:2380263/150‑1 (MQ=255)
                                                                          gcgcTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTg                                                                           <  1:3808490/150‑1 (MQ=255)
                                                                          gcgcTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTg                                                                           <  1:1324291/150‑1 (MQ=255)
                                                                           cgcTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTgc                                                                          >  1:995034/1‑150 (MQ=255)
                                                                            gcTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTgcg                                                                         >  2:3727264/1‑150 (MQ=255)
                                                                            gcTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTgcg                                                                         <  1:4115481/150‑1 (MQ=255)
                                                                            gcTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTgcg                                                                         <  1:1864935/150‑1 (MQ=255)
                                                                              tGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCa                                                                       <  2:1842121/150‑1 (MQ=255)
                                                                                gAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAg                                                                     <  1:3124591/150‑1 (MQ=255)
                                                                                gAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAg                                                                     <  2:1402597/150‑1 (MQ=255)
                                                                                 aaTGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGa                                                                    >  1:1942017/1‑150 (MQ=255)
                                                                                 aaTGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGa                                                                    >  2:3829861/1‑150 (MQ=255)
                                                                                  aTGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGaa                                                                   >  1:2297115/1‑150 (MQ=255)
                                                                                    gctgcGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAAt                                                                 <  2:852683/150‑1 (MQ=255)
                                                                                     ctgcGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATg                                                                <  2:2105990/150‑1 (MQ=255)
                                                                                     ctgcGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATg                                                                >  2:4509041/1‑150 (MQ=255)
                                                                                     ctgcGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATg                                                                <  1:3735463/150‑1 (MQ=255)
                                                                                         gTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCgg                                                            <  1:622859/150‑1 (MQ=255)
                                                                                            aTGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGaaa                                                         <  2:1154570/150‑1 (MQ=255)
                                                                                            aTGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGaaa                                                         >  1:3971731/1‑150 (MQ=255)
                                                                                               ctgcCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGaaaaaa                                                      >  2:3303880/1‑150 (MQ=255)
                                                                                                tgcCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAAt                                                     <  2:1336721/150‑1 (MQ=255)
                                                                                                tgcCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAAt                                                     <  2:2994931/150‑1 (MQ=255)
                                                                                                 gcCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATc                                                    <  1:4522553/150‑1 (MQ=255)
                                                                                                 gcCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATc                                                    >  1:1138247/1‑150 (MQ=255)
                                                                                                     ttGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCCCAGTGCGCAAGAAATGCCGGAAAAAATCCGGc                                                >  1:3343822/1‑150 (MQ=255)
                                                                                                        aGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCaaa                                             >  2:1745169/1‑150 (MQ=255)
                                                                                                        aGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCaaa                                             >  1:1411729/1‑150 (MQ=255)
                                                                                                                    aTCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAg                                 <  2:3487959/150‑1 (MQ=255)
                                                                                                                    aTCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAg                                 <  1:2167787/150‑1 (MQ=255)
                                                                                                                    aTCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAg                                 <  2:1259816/150‑1 (MQ=255)
                                                                                                                     tCTTAACGTCAGTCAAAGCTCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGt                                <  2:2915698/150‑1 (MQ=255)
                                                                                                                      cTTAACGTCAGTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGtt                               >  1:2201492/1‑150 (MQ=255)
                                                                                                                       ttAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTg                              >  1:4386220/1‑150 (MQ=255)
                                                                                                                       ttAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTg                              >  1:2428742/1‑150 (MQ=255)
                                                                                                                             tCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCa                        >  2:4262394/1‑150 (MQ=255)
                                                                                                                               aGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATg                      >  2:2144331/1‑150 (MQ=255)
                                                                                                                                gTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATgg                     <  2:4144634/150‑1 (MQ=255)
                                                                                                                                gTAAAAGCGCTTGGTTAACGGAGATCCCGCTGGCGAAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATgg                     <  2:4449621/150‑1 (MQ=255)
                                                                                                                                  aaaaGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGc                   <  2:2800522/150‑1 (MQ=255)
                                                                                                                                  aaaaGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGc                   <  2:1962547/150‑1 (MQ=255)
                                                                                                                                   aaaGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCg                  <  1:2384188/150‑1 (MQ=255)
                                                                                                                                    aaGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCgg                 >  2:2989334/1‑150 (MQ=255)
                                                                                                                                     aGCGCTTGGTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCggg                <  2:4081879/150‑1 (MQ=255)
                                                                                                                                       cgcTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTa              <  1:1197008/150‑1 (MQ=255)
                                                                                                                                       cgcTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTa              <  1:42370/150‑1 (MQ=255)
                                                                                                                                          ttGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACt           >  2:2238360/1‑150 (MQ=255)
                                                                                                                                          ttGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACt           <  1:142088/150‑1 (MQ=255)
                                                                                                                                           tGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCACATCGAAGAGAAAGTTGCCCGCATGGGCGGGTCACTa          >  2:3566131/1‑150 (MQ=255)
                                                                                                                                           tGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCTAAGAGAAAGTTGCCCGCATGGGCGGGTAACTa          >  2:3450548/1‑150 (MQ=255)
                                                                                                                                            ggTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTAc         >  1:1054456/1‑150 (MQ=255)
                                                                                                                                            ggTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTAc         <  1:868064/150‑1 (MQ=255)
                                                                                                                                            ggTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTAc         <  1:867702/150‑1 (MQ=255)
                                                                                                                                            ggTTAACGGAGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTAc         <  1:426696/150‑1 (MQ=255)
                                                                                                                                                 aCGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTAc    <  1:2075186/150‑1 (MQ=255)
                                                                                                                                                 aCGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTAc    <  1:389306/150‑1 (MQ=255)
                                                                                                                                                   ggCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAg  <  1:2683749/150‑1 (MQ=255)
                                                                                                                                                   ggCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAg  >  2:742796/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTTTTATTTATCTGAACCGCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAG  >  CP009273/4612632‑4612928

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: