Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,514,754 A→T 9.5% intergenic (+67/+198) lysV → / ← xapR tRNA‑Lys/transcriptional activator of xapAB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,514,7540AT 9.5% 386.4 / 7.3 168intergenic (+67/+198)lysV/xapRtRNA‑Lys/transcriptional activator of xapAB
Reads supporting (aligned to +/- strand):  ref base A (109/33);  new base T (2/16);  total (115/53)
Fisher's exact test for biased strand distribution p-value = 7.75e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

TTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGA  >  CP009273/2514609‑2514902
                                                                                                                                                 |                                                                                                                                                    
 tCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                                                                                 <  2:259402/150‑1 (MQ=33)
 tCGGGTCGTGAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt                                                                                                                                                 <  2:213394/150‑1 (MQ=18)
             cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                                                                     <  2:3960972/150‑1 (MQ=32)
              tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                                                                                                                    >  2:3380343/1‑150 (MQ=32)
              tCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTGTTACTCTCTGCGATACTCCAAGATTCGAACCTGCAGCAGGGGTGAGGTGAgcgctg                                                                                                                                    >  2:3676980/1‑148 (MQ=11)
               cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGc                                                                                                                                   <  2:3302536/150‑1 (MQ=32)
                aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                                                                                                                  <  1:2725878/150‑1 (MQ=32)
                aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                                                                                                                  <  2:3369216/150‑1 (MQ=32)
                aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                                                                                                                  >  2:1905356/1‑150 (MQ=32)
                aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTAATATATTCGATAATCAAAGATTAGAACCTGCCGCAGGTTTGAGTTGTGCGCAGCg                                                                                                                                  >  1:822197/1‑150 (MQ=11)
                   tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAACAAGCGCCCTACAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                                                                               <  2:4212714/150‑1 (MQ=17)
                   tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                                                                               <  1:1576472/150‑1 (MQ=32)
                    ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                                                                              <  1:971087/150‑1 (MQ=32)
                    ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                                                                              >  2:292687/1‑150 (MQ=32)
                    ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                                                                              <  1:970985/150‑1 (MQ=32)
                    ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                                                                              >  1:2955003/1‑150 (MQ=32)
                     gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaaca                                                                                                                             <  2:2774010/150‑1 (MQ=17)
                     gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaaca                                                                                                                             <  2:381739/150‑1 (MQ=32)
                      tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTAGAACCTTCAGCAGGTTTGAGTTGAGCGCAGCGAaacaa                                                                                                                            >  1:4136880/1‑150 (MQ=17)
                         agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                                                                                                                         >  1:1716882/1‑150 (MQ=32)
                          gctgTTGAGTTTTAATCATTTGGTCGCTGGTTCGAATCCGGCACGACCCACCATTGTAAAAAAGCGCCCTAACGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGa                                                                                                                        <  1:2278510/147‑1 (MQ=11)
                             gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAAAGGAGcc                                                                                                                     >  2:1694196/1‑150 (MQ=17)
                               tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCCGGTTTGAGTTGACCGCAGCGAAACAACGGAGCCGc                                                                                                                   >  2:1466201/1‑150 (MQ=17)
                                 aCTTTTATTCAATGGGTCGCTGTTTCTAATCCTGCACGACCCACAATGGTAAAAACGCGCCCTAAGGGCGATTTTTTACTATCTGCGATACTAAAATTTTCGAACCTGCAGCAGTTTTGTGTTGTGCGCTGCTAACCAACGGAGCCGCTc                                                                                                                 <  2:3846103/150‑1 (MQ=255)
                                   ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGACGGTTTGAGGTGAGCGCAGCGACCCACCGGAGCCGCTcgg                                                                                                               >  2:3599466/1‑149 (MQ=11)
                                    ttaaaTCATTTGGTCGCTGGTTCTAATCCTGAACCACCCACCAATTTAAAAACGCGCCCTAAAGGCGATTTTTTAATATATGCGATAATAAATTATTCGAACCTGAAGAAGTTTTTTTTTGAGCGAAGCTAAACAACGGAGCCGCTcgcg                                                                                                              <  2:4520046/147‑1 (MQ=2)
                                     ttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                                                                                                             >  1:2560056/1‑150 (MQ=32)
                                      tAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAATCAAAGATTCGAAACTGCAGAATGTTTGAGTTGAGCTCAGAGAAACAACTCATCCGCACtcggc                                                                                                            >  2:1599161/1‑150 (MQ=11)
                                           aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                                                                       >  2:1342112/1‑150 (MQ=32)
                                              tGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTTCGATACTCAAAGATTCGAAACTGCATCATGTTTGAGGTGAGCGCAGCGAAACAACGGATCCGCTCGCGGCGACGGccc                                                                                                    >  1:343699/1‑150 (MQ=11)
                                               ggTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCg                                                                                                   >  1:3557770/1‑150 (MQ=32)
                                                  cGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAg                                                                                                >  2:1396725/1‑150 (MQ=32)
                                                    cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGTTTCCAACCTGCAGACGCTTTGAGGTGAGCTCAGAGAAACAACTGAGCCGCACGCGGCGACGGCCCGAAggg                                                                                              >  2:2396477/1‑150 (MQ=11)
                                                    cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                                                                                              >  1:1546060/1‑150 (MQ=32)
                                                    cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                                                                                              >  1:3955825/1‑150 (MQ=32)
                                                       gTTCGAATCCGCCACCACCCACCAATTTAAAAAAGCGCCATAAGGGCGCTTTTTTACTATCTGCGATACTAAACGATTCGAACCTGCAGCTGGTTTGTGTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcga                                                                                           <  2:2560056/150‑1 (MQ=11)
                                                        ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGACCCTGCAGCAGGTTTGAGTTGAGCGCAGCGACACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                                                                          >  2:1820783/1‑150 (MQ=17)
                                                         tCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagc                                                                                         >  2:4553351/1‑150 (MQ=32)
                                                          cGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGAGAAACAACTGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcg                                                                                        >  2:811382/1‑150 (MQ=17)
                                                          cGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAATCCAATATTCGACCCTGCCGGAGGTTTGAGTTGAGCGCATCGAAACAACGGAGCCGCACGCGGCGACGGCCCGAAGGgcgcgcg                                                                                        >  2:811388/1‑150 (MQ=11)
                                                           gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                                                                       >  2:1817591/1‑150 (MQ=32)
                                                           gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                                                                       >  1:2993333/1‑150 (MQ=32)
                                                            aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGCACCTGCATCACGGTTGAGTTGAGCCCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGCGCGagcgaa                                                                                      >  1:2492433/1‑150 (MQ=11)
                                                            aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaa                                                                                      >  1:1960832/1‑150 (MQ=32)
                                                             aTCCGGCACGACCCACCACTGTAAAAAACCGCCCCAAGGGCGCTTTTTTACTATCTGCGATACTCAACTATTCGAACCTGCCGCCGGTGTGTGTTGCGCCCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaag                                                                                     <  2:2607905/150‑1 (MQ=11)
                                                               ccTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCCAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCCAAGGGCGagcgaagcg                                                                                   >  2:677157/1‑150 (MQ=17)
                                                                ctccaccacCCCCAAATTAAAAAACGCCCCCAAACGGCGTTTTTTTAATATTTGCTAAACACAAATTTTCCAACCTGCAGAAGGTTTTTGTTGAGCGCCGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                                                                                  <  1:4330332/147‑1 (MQ=2)
                                                                   cACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTc                                                                               >  1:2304504/1‑150 (MQ=32)
                                                                     cGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCAt                                                                             >  2:1455326/1‑150 (MQ=32)
                                                                     cGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAACTTTCCACCCTGCCGCAGCTTTGAGTTGTGCGAAGAGAAACAACAGATCCGATCGCGGCGACGGCCCCAAGGCCGAGCCAATCGAGTcca                                                                             >  1:2635824/1‑148 (MQ=11)
                                                                       aCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATcc                                                                           >  1:278159/1‑150 (MQ=32)
                                                                                    aaaaaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTTATTTGAGCGCAGCGAAACAACGGATCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCAGCACGAcccccc                                                              >  1:3154592/1‑150 (MQ=11)
                                                                                     aaaaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca                                                             >  2:2006414/1‑150 (MQ=32)
                                                                                     aaaaGCGCCCGAACGGCGATTTTTTAATATCTGCTATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca                                                             <  2:4209783/150‑1 (MQ=11)
                                                                                      aaaGCGCCCTAACGGCGATTTTTTACTATCTGCGATACTCAAAGTTTCTAACCTGCAGCAGTTTTGTGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa                                                            <  2:4243607/150‑1 (MQ=11)
                                                                                        aGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATg                                                          >  1:3812848/1‑150 (MQ=32)
                                                                                           gcCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGTTTCGAACCTGCAGCAGGTTTGCGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa                                                       <  2:1934027/150‑1 (MQ=17)
                                                                                           gcCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa                                                       >  2:4334881/1‑150 (MQ=32)
                                                                                           gcCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCAGCACGACCCACCAAAGTaa                                                       >  1:3710518/1‑150 (MQ=17)
                                                                                           ccccTAACGGCGCTTATTAACTATCTGCTTTACTAAATTCTTCCACCCGCCTGATGTTTTGAGTTGTGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa                                                       <  2:4435884/149‑1 (MQ=2)
                                                                                             ccTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa                                                     >  2:2407109/1‑150 (MQ=32)
                                                                                              cTAAAGGCGCTTTTTTAATATATGCGATAATCAAAGATTCGAACCTGCAGCAGGTTTTAGTTTAGCGCAGCGAAACAACGGAGCCGCTACCGGCGACGGCCCCAAGGGCGAGCGAACCCAGTAAAACAGCACCACCCACCAAAGTaaaaa                                                    >  1:2810877/1‑150 (MQ=11)
                                                                                               tAATGGGGCTTTTTTAATTTGTGCTATACTCAATGATTCGAACCTGCTGCTGTTTTGAGTTGAGCGCCGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa                                                   <  2:2168665/150‑1 (MQ=11)
                                                                                               tAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaaa                                                   >  1:4452836/1‑150 (MQ=32)
                                                                                                  tggCGTTTTTTTACTATATGCTATACTAAATTTTTCTAACCTGCTGCTGGTTTTTTTTGAGCGCAGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAgcg                                                <  1:4389389/149‑1 (MQ=11)
                                                                                                  aGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAgcg                                                >  1:1008345/1‑150 (MQ=32)
                                                                                                  aGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAgcg                                                >  1:1008334/1‑150 (MQ=32)
                                                                                                   ggCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGCGTCATCCTGCACGACCCACCAATGTAAAAAAgcgc                                               >  2:4633022/1‑150 (MQ=17)
                                                                                                   ggCGCTTTTTTACTATCTGCGATACTCAAAGATTAGAACCTGCAGCAGGGTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATACTGCACGACCCACCAATGTAAAAAAgcgc                                               >  1:3707001/1‑150 (MQ=12)
                                                                                                    tgtatttttaAATATCTGCCATACTAAAAGTTCCCACCCTGCTGCGGTTGTGTGTGTCGCCCCGCCAAACCACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGcc                                              <  1:2217118/146‑1 (MQ=2)
                                                                                                        ttttttACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                          >  2:884055/1‑150 (MQ=32)
                                                                                                        ttttttACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                          >  2:4192390/1‑150 (MQ=32)
                                                                                                            ttACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                      >  1:1487249/1‑150 (MQ=32)
                                                                                                            ttACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgc                                      >  1:2277215/1‑150 (MQ=32)
                                                                                                             tACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcg                                     >  2:2573161/1‑150 (MQ=31)
                                                                                                             tACTATCTGCGATACGCAAAGATTCGAACCTGCGGCTGGTTTGAGTTGAGCGCAGCGAAACCACGGAGCCGCTCGCGGCGACGGCCCGTCGGGGGAGCGAAGCGAGTCATCCTGCACGACCCTCCAATGTACAGAAGCGCCCTATAGgcg                                     >  2:2573074/1‑150 (MQ=11)
                                                                                                              aCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcgc                                    >  2:4025043/1‑150 (MQ=11)
                                                                                                                           cTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCt                       >  1:185959/1‑150 (MQ=31)
                                                                                                                           cTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCt                       >  1:2578471/1‑150 (MQ=31)
                                                                                                                            tCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTg                      >  1:708737/1‑150 (MQ=32)
                                                                                                                             cAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGc                     >  2:1719952/1‑150 (MQ=32)
                                                                                                                              aatgATTCCACCCGGCAGCTGGTTGGCGTGGCGCCCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCg                    <  1:2091343/147‑1 (MQ=11)
                                                                                                                              aaaGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCg                    >  2:4749268/1‑150 (MQ=33)
                                                                                                                                 gATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAtt                 >  2:4337802/1‑150 (MQ=35)
                                                                                                                                  aTTCGAACCTGCCGCAGGTTTGCGTTGAGCGCCGCGAAACCACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAttt                <  2:3512209/150‑1 (MQ=255)
                                                                                                                                  aTTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAttt                >  1:2942069/1‑150 (MQ=37)
                                                                                                                                  aTTCGAACCTGACGAGGGTTGTAGTTGCGCGCCGCAAACAAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCATTTTTGCTATCTGCGAttt                <  1:3079891/150‑1 (MQ=11)
                                                                                                                                   ttCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTg               >  1:3274476/1‑150 (MQ=38)
                                                                                                                                   ttCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTg               >  2:3002113/1‑150 (MQ=38)
                                                                                                                                     cGACCCTGCAGCAGGTTTGATTTGCGCCCAGCCAACCAACGGCGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCACCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCg             <  1:4707199/150‑1 (MQ=14)
                                                                                                                                      gcaCCTGCAGCAGTTTCGAGTTGACCGCTCCTAAACAACGGAGCCTCTCGCGGCGACGGCCCGAAGGGCGATCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGa            <  2:3004606/148‑1 (MQ=14)
                                                                                                                                       aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa           >  2:1753271/1‑150 (MQ=40)
                                                                                                                                       aaCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaa           >  1:1961439/1‑150 (MQ=40)
                                                                                                                                        aCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGaaa          >  1:2898362/1‑150 (MQ=255)
                                                                                                                                         ccTGCTGCTGGTTGGAGTTTAGCGCCGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAt         <  1:4388753/150‑1 (MQ=14)
                                                                                                                                          cTGCAGCTGTATTGATTTGTGCGCAGCTAAACCACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt        <  1:2323961/150‑1 (MQ=14)
                                                                                                                                          cTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAAtt        >  2:2207149/1‑150 (MQ=255)
                                                                                                                                              agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCt    >  2:1996198/1‑150 (MQ=255)
                                                                                                                                              agcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCt    >  2:935416/1‑150 (MQ=255)
                                                                                                                                                cagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGa  >  2:3494387/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGA  >  CP009273/2514609‑2514902

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: