Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,602,056 T→G 5.5% intergenic (+47/+55) zntA → / ← tusA zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,602,0560TG 5.5% 462.8 / 4.7 188intergenic (+47/+55)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (121/52);  new base G (3/11);  total (125/63)
Fisher's exact test for biased strand distribution p-value = 4.77e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

ACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATA  >  CP009273/3601908‑3602204
                                                                                                                                                    |                                                                                                                                                    
acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTg                                                                                                                                                     <  2:4453079/150‑1 (MQ=255)
  gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   >  2:389368/1‑150 (MQ=255)
   cTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                  <  2:4424491/150‑1 (MQ=255)
      ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               >  1:1845193/1‑150 (MQ=255)
             tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTGAg                                                                                                                                        >  2:1027992/1‑150 (MQ=255)
              gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCTACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                       <  2:4275095/150‑1 (MQ=255)
               aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                      >  2:4041975/1‑150 (MQ=255)
                   ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                  >  2:1805868/1‑150 (MQ=255)
                   ggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                  <  1:4038914/150‑1 (MQ=255)
                     ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                <  1:314258/150‑1 (MQ=255)
                      tgtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGTGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                               >  1:3655447/1‑150 (MQ=255)
                       gtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGa                                                                                                                              >  1:692936/1‑150 (MQ=255)
                       gtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGa                                                                                                                              >  1:2924885/1‑150 (MQ=255)
                          gcgggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                           <  2:55687/147‑1 (MQ=255)
                           ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                                          >  1:1932383/1‑150 (MQ=255)
                           ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                                          >  1:2629876/1‑150 (MQ=255)
                            tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                         >  2:4278695/1‑150 (MQ=255)
                            tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCACTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                         >  1:11972/1‑150 (MQ=255)
                             ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCACCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggc                                                                                                                        >  2:809781/1‑150 (MQ=255)
                             ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggc                                                                                                                        >  2:3575119/1‑150 (MQ=255)
                                   gCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATc                                                                                                                  >  1:24661/1‑150 (MQ=255)
                                    cTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCACCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCg                                                                                                                 >  2:3416832/1‑150 (MQ=255)
                                     tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                >  2:3548081/1‑150 (MQ=255)
                                     tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATGGTTGCGCAGGAGAGAAGGCAAACCGATCGCGACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGCGAGGGTTAGGGTGAGGGGAAGAGGCGGCATCGa                                                                                                                >  1:276093/1‑150 (MQ=39)
                                      ggCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                               >  2:611809/1‑150 (MQ=255)
                                           aTACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCaa                                                                                                          >  1:3562479/1‑150 (MQ=255)
                                               gggggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACAGTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGGGGGGAAAAGGCGGCATCGAAGCCAGTCAg                                                                                                      >  1:3834781/1‑150 (MQ=255)
                                                    cgccGGTGGTGGTGACCGCCATTGCTTTAATTTTTTTGCCCATGACATACGCCAAACCGATCGAAACATTGCGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGGGAGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                 <  2:4537955/147‑1 (MQ=19)
                                                       cgttGCGGGTGACAGCGATTGCGTTAAGATTGTGGGCCTGTCGAAAAGTCAAACCGACCGCAACATTGAGCGCTATCGGCCCCCTCGCCCCTGTGGGGTGCGGGTTGGGTTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATc                                                                                              <  2:1397080/147‑1 (MQ=255)
                                                        ggTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCa                                                                                             >  2:268991/1‑150 (MQ=255)
                                                        ggTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCa                                                                                             >  2:4215859/1‑150 (MQ=255)
                                                         gTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCaa                                                                                            >  2:1462674/1‑150 (MQ=255)
                                                           gCTGGTGACTGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCTTTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                          <  2:3375933/150‑1 (MQ=255)
                                                           gCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                          >  1:2520216/1‑150 (MQ=255)
                                                             tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                        >  2:4175931/1‑150 (MQ=255)
                                                              ggTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGcc                                                                                       >  2:2840901/1‑150 (MQ=255)
                                                               gTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                      >  1:4055190/1‑150 (MQ=255)
                                                               gTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                      >  1:1576866/1‑150 (MQ=255)
                                                                tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGTTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                     >  2:3843617/1‑150 (MQ=255)
                                                                  aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                   >  2:220477/1‑150 (MQ=255)
                                                                   cAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                  >  2:1250056/1‑150 (MQ=255)
                                                                      cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                               >  1:2690325/1‑150 (MQ=255)
                                                                      cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                               >  1:176101/1‑150 (MQ=255)
                                                                      cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                               >  1:176183/1‑150 (MQ=255)
                                                                       tatggtgTTATGTTTGTTGCGCAGGAGATACGGCAACCCGCTCGCAACATTGCGCGCGAGGGGTCCCCTCGCCCCTCTGGGGAGCGGGTTAGGGTGTGGGGAAAAGGCGGCATAGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                              <  2:1439992/144‑1 (MQ=37)
                                                                       gAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                              >  1:1153224/1‑150 (MQ=255)
                                                                       gAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                              >  1:1593823/1‑150 (MQ=255)
                                                                         aTGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCa                                                                            >  2:2785542/1‑150 (MQ=255)
                                                                             tttatttttgtcgCCCAGGATATATGCCAAACCGATCGCAACTTTGCGCCCTTTCGGCCCCCACGCCCCTCGGGGTATTGGTTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        <  1:4212959/143‑1 (MQ=255)
                                                                             gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  2:2608610/1‑150 (MQ=255)
                                                                             gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  2:2003166/1‑150 (MQ=255)
                                                                                  gATTGTTGCGCTGGATATATGGCAAACCGATGGAAACTTTGCGCGCTATGGGTCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGAt                                                                   <  2:2353248/150‑1 (MQ=37)
                                                                                  gATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGAt                                                                   >  2:1739686/1‑150 (MQ=255)
                                                                                   agtgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCACAGGGGAGAGGGTTAGGGGGAGGGGAAAAGGCGGCAGCGAATCCCATCAGCCCCTCTCAACCGCCTGTACGAATCACATAACGAtc                                                                  >  2:2802128/3‑149 (MQ=255)
                                                                                     tgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                                >  1:3399196/1‑150 (MQ=255)
                                                                                     tgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                                >  1:4212988/1‑150 (MQ=255)
                                                                                     tgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                                >  1:3324342/1‑150 (MQ=255)
                                                                                      gttgCGCAGTCTATACGCCAAACGCACCGCAACTTGGCGCGCTATCGGTCCCCTCCCCCCTGTGGGTAGTGGGTTAGGGTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                               <  1:4537743/150‑1 (MQ=255)
                                                                                            cAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTcc                                                         >  2:1483416/1‑150 (MQ=255)
                                                                                            atgtataTACGGCAAACCGATCGAAACATTGTGCGCTATCGGTCCCCTCGCCCCTCTGGGGTTTGGGTTAGGGTGTGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCATATAACGATAAGGCAGTcc                                                         <  2:657613/144‑1 (MQ=255)
                                                                                             tggAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                        <  2:4107306/149‑1 (MQ=255)
                                                                                             cgtcgaTACGCCAACCCGATCGCAAAATGGCGCGCGATGGGTCCCCTCGCCCCTCGGGGGAGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                        <  2:2359553/146‑1 (MQ=255)
                                                                                             aGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                        >  1:4614148/1‑150 (MQ=255)
                                                                                                ataTAAGGCAACCCGATCGCAAAATTGCGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGGGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                     <  1:4169252/148‑1 (MQ=38)
                                                                                                    acgcaaaaCCCATCCAAACATTGCGCGCGCGGGGGCCCCCCGCCCCTGGGGGGTGTGGGTTGGGGTGCGGGACCAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                 <  1:2307714/145‑1 (MQ=255)
                                                                                                    aaGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                 >  1:1417967/1‑150 (MQ=255)
                                                                                                    aaGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                 >  2:157266/1‑150 (MQ=255)
                                                                                                      ggCAAACCGATCGCAACCTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGAGAAGGCAGTCCATCCGTctct                                               >  2:1771963/1‑150 (MQ=255)
                                                                                                       gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              >  1:1505792/1‑150 (MQ=255)
                                                                                                        cAACCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                             <  2:2483763/150‑1 (MQ=255)
                                                                                                        cAACCCGATCGCAACATTGAGCGCGATCGGGCCCCGCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                             <  2:597892/150‑1 (MQ=255)
                                                                                                        cAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCACTCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                             >  2:3363626/1‑150 (MQ=255)
                                                                                                         accccGATCGCAACATTGACCGCGATCGGCCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            <  2:2714255/147‑1 (MQ=255)
                                                                                                         aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            >  1:581484/1‑150 (MQ=255)
                                                                                                             cGATCGCAACATTGAGCGCGATCGGTCCCCNCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCaa                                        >  1:4283978/1‑150 (MQ=255)
                                                                                                              gaccGCACCATGGCGCGCGCTCGGCCCCCGCGCCCCTCGGGGGAGCGGGTTCGGTTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                       <  1:3710861/147‑1 (MQ=255)
                                                                                                              gATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCacc                                       >  2:2374584/1‑148 (MQ=255)
                                                                                                              gATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                       >  2:3120645/1‑150 (MQ=255)
                                                                                                               aTCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcc                                      <  2:3768308/150‑1 (MQ=255)
                                                                                                                 cGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAg                                    >  1:165100/1‑150 (MQ=255)
                                                                                                                     aCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCg                                <  2:502106/150‑1 (MQ=255)
                                                                                                                     aCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCg                                >  1:1263906/1‑150 (MQ=255)
                                                                                                                       aTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                              >  2:4609053/1‑150 (MQ=255)
                                                                                                                         tGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                            <  2:1143394/150‑1 (MQ=255)
                                                                                                                           aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                          >  2:3762253/1‑150 (MQ=255)
                                                                                                                                gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                     >  2:2134190/1‑150 (MQ=255)
                                                                                                                                  tCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                   >  2:3775439/1‑150 (MQ=255)
                                                                                                                                    ggccccctcgcccctcgGGGGAGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa                 <  1:2557435/147‑1 (MQ=255)
                                                                                                                                    ggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa                 >  1:4220001/1‑150 (MQ=255)
                                                                                                                                         ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa            >  2:682036/1‑150 (MQ=255)
                                                                                                                                         ccTCGCCCATCTGGGGAGTGGTTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa            <  2:1663124/150‑1 (MQ=255)
                                                                                                                                          cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAGAAAGGTACAAAAc           >  2:3401018/1‑150 (MQ=255)
                                                                                                                                            cGCCCCTGTGGGGAGTGGGTTCGGGTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAccc         <  2:3721983/150‑1 (MQ=255)
                                                                                                                                            cGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAccc         >  1:928584/1‑150 (MQ=255)
                                                                                                                                             gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        >  1:3059323/1‑150 (MQ=255)
                                                                                                                                             gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        >  1:743636/1‑150 (MQ=255)
                                                                                                                                               cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAgg      >  2:3174138/1‑150 (MQ=255)
                                                                                                                                                  tgtGGGGTGTGGGTTTGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat   <  1:105296/148‑1 (MQ=255)
                                                                                                                                                   ctGGGGAGAGGGTTCGGGTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAata  <  1:2054291/150‑1 (MQ=255)
                                                                                                                                                   ctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAata  >  1:4404661/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATA  >  CP009273/3601908‑3602204

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