Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 169,937 C→A 26.0% intergenic (+6/+152) fhuB → / ← hemL fused iron‑hydroxamate transporter subunits of ABC superfamily: membrane components/glutamate‑1‑semialdehyde aminotransferase (aminomutase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273169,9370CA26.0% 586.8 / 33.6 252intergenic (+6/+152)fhuB/hemLfused iron‑hydroxamate transporter subunits of ABC superfamily: membrane components/glutamate‑1‑semialdehyde aminotransferase (aminomutase)
Reads supporting (aligned to +/- strand):  ref base C (26/138);  new base A (61/2);  total (112/140)
Fisher's exact test for biased strand distribution p-value = 1.02e-31
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

GCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTG  >  CP009273/169788‑170086
                                                                                                                                                     |                                                                                                                                                      
gcgcTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTc                                                                                                                                                        <  2:868458/150‑1 (MQ=255)
gcgcTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTc                                                                                                                                                        <  2:46914/150‑1 (MQ=255)
gcgcTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTc                                                                                                                                                        <  2:1773381/150‑1 (MQ=255)
  gcTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCa                                                                                                                                                      <  1:4631347/150‑1 (MQ=255)
  gcTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCa                                                                                                                                                      >  2:2287539/1‑150 (MQ=255)
    tggtggGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTACaaa                                                                                                                                                    >  1:3787660/1‑148 (MQ=255)
    tggtggGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTGCata                                                                                                                                                    >  1:4677108/1‑150 (MQ=255)
     ggtgggtgGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCatat                                                                                                                                                   <  1:4876010/150‑1 (MQ=255)
      gtgggtggTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTACAaata                                                                                                                                                  >  1:2381155/1‑150 (MQ=255)
       tgggtggTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATaa                                                                                                                                                 <  1:2528894/150‑1 (MQ=255)
          gtggtTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTGTTTTTTTAAATATaaaag                                                                                                                                              >  1:3348718/1‑147 (MQ=255)
                 agtggtgGTGTTGGCTGATTGGTGGGGGGGTTTGGTGTTGTTTCCTTTCCAGATCCGGGGGGGGCTGATGTCAACATTTATGGGCGCGCATTATTTTATCTATTTGTTTAGAAAGCAGAGCCTTTAATTTTTCCATATAACGGCCCTGTa                                                                                                                                      <  1:87566/145‑1 (MQ=255)
                  ctgctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAAAAGCCCAGTaa                                                                                                                                     >  1:3187471/1‑149 (MQ=255)
                   tgctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTGTATTTTTAAAAATAAAAGCCCAGTAca                                                                                                                                    >  1:2077705/1‑149 (MQ=255)
                    gctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTTTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGc                                                                                                                                   <  1:1933968/150‑1 (MQ=255)
                    gctgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGc                                                                                                                                   <  2:715439/150‑1 (MQ=255)
                      tgGTGTTCGCTGAGTGGTGTGGGGGGTTGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACATTTATCGGCGCGCAATATTTTATCTATTTGTTGAGAAAGCAGAGCCTTTAATTTTTCCATATCACGGCCCTGTACGCCt                                                                                                                                 <  1:2939997/150‑1 (MQ=255)
                      tgGTGTTCGCTGAGTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCt                                                                                                                                 <  1:3861324/150‑1 (MQ=255)
                      tgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTACATAAAAAGGCCCGGTACGcca                                                                                                                                 >  1:1467884/1‑149 (MQ=255)
                      tgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCt                                                                                                                                 <  1:2230465/150‑1 (MQ=255)
                      tgGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTGAATTTTTACATATAAAAGCCCGGTAAGcca                                                                                                                                 >  2:656609/1‑149 (MQ=255)
                       ggTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGTGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAAATTTTAAATAAAAAGGACCTGTAAGCCTg                                                                                                                                >  2:4333631/1‑150 (MQ=255)
                       ggTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTGTATTTTTTCTTATTAACGCCCCGTTCGccag                                                                                                                                >  1:80626/1‑148 (MQ=255)
                        gtgtTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTgg                                                                                                                               <  1:1514621/150‑1 (MQ=255)
                        gtgtTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAAACGCCCTGTAAGCCTgg                                                                                                                               >  2:2977637/1‑150 (MQ=255)
                         tgtTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTACATATAAAAGCCCAGTTAGCCTggt                                                                                                                              >  1:4319347/1‑149 (MQ=255)
                          gtTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTACATATAACCGCCCTGTACGCCTGGAc                                                                                                                             >  2:1910941/1‑150 (MQ=255)
                             cGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTAAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGa                                                                                                                          <  1:214323/150‑1 (MQ=255)
                              gCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGAt                                                                                                                         <  2:1831644/150‑1 (MQ=255)
                               cTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATa                                                                                                                        <  2:965379/150‑1 (MQ=255)
                                tGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTGAATATTTTCATAAAAACACCCCGTAAGCccggccggcaaa                                                                                                                       >  1:4092706/1‑146 (MQ=255)
                                        gtgGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTa                                                                                                               <  2:792509/150‑1 (MQ=255)
                                         tgGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTTAATATAAAGGCCCAGTAAGCCAGGACCGATAAAGAGtttta                                                                                                              >  1:3092913/1‑148 (MQ=255)
                                         tgGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTATGCCTGGACGGATAAAGCGTTTAc                                                                                                              <  1:3593139/150‑1 (MQ=255)
                                             cGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTTAATATAAAGGCCCGGTACGCCAGGCCGGCTAAAGCGGTTAAGCGa                                                                                                          >  1:3502500/1‑150 (MQ=255)
                                               gATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACa                                                                                                        <  1:295295/150‑1 (MQ=255)
                                               gATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACa                                                                                                        <  1:1188683/150‑1 (MQ=255)
                                                  ggTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTCTCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCAGATCACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACg                                                                                                     <  1:2177974/150‑1 (MQ=255)
                                                  ggTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACg                                                                                                     <  2:1178608/150‑1 (MQ=255)
                                                       tGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTTCAGATAAAGGCCCAGTAAGCCAGGACGGCTAAAGCGTTTTAGCGACAACCGATGc                                                                                                >  1:4716661/1‑150 (MQ=255)
                                                         tttCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTACATATAAAGGCCCGGTAAGCCGGGACGGATAAAGCGTTTAAGCGACAACGGATGccc                                                                                              >  2:3428234/1‑150 (MQ=255)
                                                         tttCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGccc                                                                                              <  2:3389868/150‑1 (MQ=255)
                                                             cattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTGCAGATAAACGCCCGGGACCCCGGGGCGGGGAAACCCGTTACGCCCCAACCGAGGcccccag                                                                                          >  2:4166012/1‑146 (MQ=255)
                                                             cattccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACAAAAAAAAGCCCATGACGCCAGGACGGATAAATAGGTTTAGCGACAACCGCTGCCCCATg                                                                                          >  1:1095354/1‑150 (MQ=255)
                                                                tccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTGAATTTTTGCATAAAAACGCCCCGGAAGCCCGGACGGGGAAAAAGGTTACGCGCCCACCGCTGcccccggcga                                                                                       >  1:1395337/1‑143 (MQ=255)
                                                                accaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTGGAGAACGCAGCGCCTTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGa                                                                                       <  1:1612277/149‑1 (MQ=255)
                                                                 ccaGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTGAATTTTTGCATAAAAACGCCCTGTAAGCCAGGGAGGATAAAACGGTTAAGCGACAACCGATGCCCCCTGCGAc                                                                                      >  1:1522587/1‑150 (MQ=255)
                                                                  caGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTTCATAAAAAAGCCCGGTAAGCCTGGACGGATAAAGCGTTTAAGCGACAACCGAATCCCCCTGCGaag                                                                                     >  1:1451512/1‑148 (MQ=255)
                                                                  caGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACg                                                                                     <  2:3286919/150‑1 (MQ=255)
                                                                  caGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACg                                                                                     <  1:75474/150‑1 (MQ=255)
                                                                  caGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTGAATTTTTTAAAAAAACTGGCCAGTAAGCCTGGTCGGTGAAATCGTTTAATCGACAAACGATGCCCCAAGCGACg                                                                                     >  1:978443/1‑150 (MQ=255)
                                                                    gATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTTCATATAAAGGCCCTGTAAGCCAGGGCGGGTAAATCGTTTTAGCGACAACGGGTGCCCGATGCGAccca                                                                                   >  1:2730159/1‑147 (MQ=255)
                                                                       cccGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGgc                                                                                <  1:4577332/150‑1 (MQ=255)
                                                                       cccGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAAAGGCCCAGTACGCCTGGACGGCTAAAGCGTTTACGCGACAACGGATGCCCGCTGCGACGCTGgc                                                                                >  2:2796476/1‑150 (MQ=255)
                                                                          ggcggGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTCCCAAAAAAGGCCCGGTACGCCAGGCCGGATAAAGCGGTTACGCGACCACCGATGCCCCAGGCGACGCTGGcgcg                                                                             >  2:2130348/1‑150 (MQ=255)
                                                                             ggggCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACAAATAACGGCCCTGTACGCCTGGACGGATAAATCGGTTAAGCGACAACCGATGCCCGATGCGACGCTGGCGCGTCt                                                                          >  2:3077377/1‑150 (MQ=255)
                                                                               ggCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTTCATAAAAAGGCCCAGTACGCCAGGACGGATAAATCGGTTTAGCGACAAACGATTACCCAAGCGACGCTGGCGCGTCTTa                                                                        >  2:164694/1‑150 (MQ=255)
                                                                               ggCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTa                                                                        <  2:4471341/150‑1 (MQ=255)
                                                                               ggCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTa                                                                        <  1:2778370/150‑1 (MQ=255)
                                                                                 ctgctgTCACCCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATc                                                                      <  2:3513163/150‑1 (MQ=255)
                                                                                  tgctgTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCg                                                                     <  1:656566/150‑1 (MQ=255)
                                                                                   gctgTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCgg                                                                    <  2:2603487/150‑1 (MQ=255)
                                                                                      gTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATATAACGGCCCTGTACGCCTGGACGGATAAAGCGGTTACGCGCCAACGGATGCCCCATGCGACGCTGGCGCGTCGTATAGGGcc                                                                 >  2:4784196/1‑150 (MQ=255)
                                                                                       tCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCt                                                                <  1:2394607/150‑1 (MQ=255)
                                                                                       tCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACAAAAAAAAGCCCAGTACACCAGGTCCGATAAAACCTTTAAGCGACAAACGATGACCCAAGCGACGCAGGCGCGGCATAAAGGGccg                                                                >  1:341549/1‑149 (MQ=255)
                                                                                            cTTTATCGGCGCTCAATATTTTATCTATTTTTTGAGAAAGCAGAGCCTTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACaaa                                                           <  1:1335108/150‑1 (MQ=255)
                                                                                            cTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAAATTTTCCATATAAAGGCCCAGTACGCCTGGACGGATAAAGCGTTTAAGCGACAACGGATTCCCGATGCGACGCTGGCGCGTCTTATAGGGCCAACaaa                                                           >  2:2393863/1‑150 (MQ=255)
                                                                                               tATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATAAAAACGCCCTGTAAGCCTGGACGGATAAATCGTTTACGCGACAAACGATTCCCGAAGCGACGCAGGCGCGTCTTATAGGGCCTACaaaacg                                                        >  2:291631/1‑147 (MQ=255)
                                                                                                       cgcCATATTTTATCTATTTGTTGATAAAGCTGCGCTGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTCCGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAAt                                                <  1:1735430/150‑1 (MQ=39)
                                                                                                       cgcCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTAATTTTACAAATAAACGCCCAGTACGCCGGGACGGGTAAAAAGTTTTACAGCCAAACCATTCCCCAATCGAAGCTGGCGCGGCGTATAGGGCCAAAAAAAAGccccccaa                                                >  1:4421218/1‑146 (MQ=255)
                                                                                                            tatTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACATAAAACGGCCCGGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATTCCCGATGCGACGCTGGCGCGTCTTATCGGGCCGAAAAACCGCCCCCAAACGTAg                                           >  2:3194210/1‑150 (MQ=255)
                                                                                                            tatTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTACAAATAAAGGCCCTCTAACCCTGGACGGATAAAGCGTTTACGCGACAAAGGATGCCCCAAGCGACCCGTGCGCGTCATATAGGGCCTACAAAAGGCCCCCAAACGTAg                                           >  2:143499/1‑150 (MQ=255)
                                                                                                            tatTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAAATTTTCCAAATAACGGCCCAGTACGCCTGGACGGATAAAGCGTTTACGCGACAAAGGATTCCCGATTCGACGCTGGCGCGTCTTATCGGGCCAACAAAAGGCCCCGaaaaggag                                           >  2:1527092/1‑144 (MQ=255)
                                                                                                                  aTCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGa                                     <  2:4825834/150‑1 (MQ=255)
                                                                                                                   tCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGAt                                    <  2:2817757/150‑1 (MQ=255)
                                                                                                                    cTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATa                                   <  1:249550/150‑1 (MQ=255)
                                                                                                                         tgttgAGAAAGCAGAGCCGTTAATTTTTACATATAAAGGCCCAGTACGCCTGGACGGATAAAGCGTTTACGCGACAAAGGATTCCCCATGCGACGCTGGCGCGTCTTATAGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCg                              >  2:4374506/1‑150 (MQ=255)
                                                                                                                          gttgAGAAAGCAGAGCCGTTAAATTTTACAAAAAAAGGCCCAGTACGCCGGGACGGAAAAAGCGTTTAAGCCACAAACGAAGCCCGAATCGACGCTTGCGCGTCATATCGGGCCTAAAAACCGCCCCCCATAGGAGGTCGGATAAGTCGt                             >  2:3082612/1‑150 (MQ=255)
                                                                                                                           ttgAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGtt                            <  2:275079/150‑1 (MQ=255)
                                                                                                                           ttgAGAAAGCAGAGCCGTTAATTTTTACAAATAAAAGCCCAGGAAGCCGGGACGGGTAAAGAGTTTAAGCGACAAACGATTCCCGAAGCGACGCAGGCGCGTCTTATACGGCCAACAAACCGCCCCCAAACGTAGGGCGGATAAGGCGtt                            >  2:453271/1‑150 (MQ=255)
                                                                                                                             gcgaACGCAGAGCCTTTAATTTTTCCATAAACCGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCa                          <  1:2455000/148‑1 (MQ=255)
                                                                                                                             gagaAAGCAGAGCCGTTAATTTTTACATATAAACGCCCGGTAAGCCGGGGCCGATAAAACGTTTAAGCGACAACGGATGGCCCCTGCGGCGCTGGCGCGGCGTATAGGGCCAACAAAAGGCCCCCCAACGGAGGTCGGGGAAGGCggtaa                          >  2:543305/1‑146 (MQ=255)
                                                                                                                              agaAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCAc                         <  2:851933/150‑1 (MQ=255)
                                                                                                                                 aaGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCAcgcc                      <  2:2165598/150‑1 (MQ=255)
                                                                                                                                     agagCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCAt                  <  2:3360116/150‑1 (MQ=255)
                                                                                                                                       cgCCGTTATTTTTTCCATATAACGGCCCGGTACGCCTGGCCGTATAAAGGGTTTACGCGACAAGGGATGCCCGATGCGACGCTGGCGCGTCTTAGCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATcc                <  1:3276156/149‑1 (MQ=38)
                                                                                                                                       cgCCGTTAATTTTTCCATAGAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATcc                <  1:1712406/149‑1 (MQ=255)
                                                                                                                                       agCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGGTAAGGCGTTCACGCCGCATcc                >  1:292624/1‑150 (MQ=255)
                                                                                                                                       agCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGGCTTATCGGGCCGACAAACGGCCCCGAATCGTAGGTCGGATAATGCGTTCACGCCGCATcc                >  1:3275613/1‑150 (MQ=255)
                                                                                                                                         ccGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa              >  1:1858777/1‑150 (MQ=255)
                                                                                                                                         ccGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa              >  1:4310819/1‑150 (MQ=255)
                                                                                                                                          cGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGAc             >  2:4486423/1‑150 (MQ=255)
                                                                                                                                          cGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGAc             <  2:4848950/150‑1 (MQ=255)
                                                                                                                                            ttAATTTTTCCATAGAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaa           <  2:197340/150‑1 (MQ=255)
                                                                                                                                            tgaaTTGTTCCCTATCCCGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaa           <  2:3977342/148‑1 (MQ=255)
                                                                                                                                              taTTTTTCCAAACACCCGGCCCTGTACGCCTGCACGGATAACGGGTTTACGCTACAACGGATGCCCGAGGCGACGCTGGCGCGTCTTATCGGGCCAACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACaaa          <  1:416485/149‑1 (MQ=37)
                                                                                                                                               aTTTTGACATATAACGGCCCTGTACGCCTGGACGTATAATGCTTTTACGCGACAACGGTTGCCCGATGCGACGCTGGCGCGTCTTATCGGGCATACAAACGGCCCCGAATCGTAGGTCGGATAAGGCTTTCACGCCGCATCCGACAAAcc        <  1:1037364/150‑1 (MQ=38)
                                                                                                                                                  tttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATg     <  2:3961974/150‑1 (MQ=255)
                                                                                                                                                  tttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAAGCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATg     >  2:3561989/1‑150 (MQ=255)
                                                                                                                                                   ttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGACTCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGc    >  2:525795/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTG  >  CP009273/169788‑170086

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: