Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,778,106 T→G 5.2% N2052T (AAC→ACC)  pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,778,1060TG 5.2% 1153.4 / 25.3 358N2052T (AAC→ACC) pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (144/190);  new base G (14/9);  total (159/199)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 363 total aligned reads displayed.

GTTGTCCTGCATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACC  >  NC_002947/4777958‑4778255
                                                                                                                                                    |                                                                                                                                                     
gTTGTCCTGCATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGtt                                                                                                                                                      <  2:5942888/150‑1 (MQ=255)
   gTCCTGCATACCGACATGGGGCAGGGTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGc                                                                                                                                                   <  2:2127352/150‑1 (MQ=255)
    tCCTGCATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCg                                                                                                                                                  <  2:3990909/150‑1 (MQ=255)
    tCCTGCATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCg                                                                                                                                                  >  2:7624745/1‑150 (MQ=255)
      cggCATACCGACTTGGGGCAGGTTCAGCACCGTCCCCCCGATCGACGCCAGTTGCTGCGCCAGCGCCGTGCGCGGAGCGACATACACCGGCCTGGCGGCGGTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACGTGCCGTTGGCGGt                                                                                                                                                <  2:4920428/148‑1 (MQ=255)
         cATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGAc                                                                                                                                             >  1:422317/1‑150 (MQ=255)
         cATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGGTGGCGGGGAc                                                                                                                                             >  1:5221099/1‑150 (MQ=255)
             ccGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGc                                                                                                                                         <  1:8532953/150‑1 (MQ=255)
             ccGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGc                                                                                                                                         >  2:2048510/1‑150 (MQ=255)
              cGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTGGGCGGTGACCGGCa                                                                                                                                        >  1:4786952/1‑150 (MQ=255)
               gACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCaa                                                                                                                                       <  1:7577813/150‑1 (MQ=255)
                  ttGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGtt                                                                                                                                    <  1:7227036/150‑1 (MQ=255)
                  ttGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGtt                                                                                                                                    <  2:5439767/150‑1 (MQ=255)
                  ttGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGtt                                                                                                                                    <  2:5928034/150‑1 (MQ=255)
                  ttGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGCTGCTCGCCAACTTGCCGTTGGCGGTGACCGGCAAGtt                                                                                                                                    <  1:4361735/150‑1 (MQ=255)
                   tGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTg                                                                                                                                   <  2:1551048/150‑1 (MQ=255)
                    ggggCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCTCGCCAACTTGCCGGTCACCGCCAACGGCAAGTTgg                                                                                                                                  <  2:5340280/150‑1 (MQ=38)
                       gCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCTCGCCAACTTGCCGGTCACCGCCAACGGCAAGTTGGCgc                                                                                                                               <  1:3582887/150‑2 (MQ=38)
                       gCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGa                                                                                                                               <  1:2397130/150‑1 (MQ=255)
                         aGGGTCAGCACCGTCTCCCAGATCGACGCCCGTTGCTGCTCCCGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGAAGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGagc                                                                                                                             <  1:6647036/150‑1 (MQ=255)
                          ggTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGagca                                                                                                                            >  1:3358073/1‑150 (MQ=255)
                           gTTCAGCACCGTCTCCCAGATCGACGCCTGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGagcag                                                                                                                           <  2:1880934/150‑1 (MQ=255)
                               aGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCCAGGTGGCGAGCAGCATc                                                                                                                       >  1:8130054/1‑150 (MQ=255)
                               aGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATc                                                                                                                       <  2:6096304/150‑1 (MQ=255)
                                gCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCa                                                                                                                      >  1:2549465/1‑150 (MQ=255)
                                gCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCa                                                                                                                      <  2:7972433/150‑1 (MQ=255)
                                   ccGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGt                                                                                                                   <  2:1052548/150‑1 (MQ=255)
                                     gTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTgc                                                                                                                 <  1:5884465/150‑1 (MQ=255)
                                     gTCACCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGAGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTAGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTgc                                                                                                                 >  2:7916190/1‑150 (MQ=255)
                                         cccAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCg                                                                                                             <  2:4302133/150‑1 (MQ=255)
                                         cccAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCg                                                                                                             >  2:4723915/1‑150 (MQ=255)
                                         cccAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCg                                                                                                             <  1:7530301/150‑1 (MQ=255)
                                          ccAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCgg                                                                                                            <  1:6296747/150‑1 (MQ=255)
                                          ccAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCgg                                                                                                            >  2:2951870/1‑150 (MQ=255)
                                          ccAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCgg                                                                                                            >  2:7869204/1‑150 (MQ=255)
                                          ccAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCgg                                                                                                            <  1:3221715/150‑1 (MQ=255)
                                            aGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCa                                                                                                          >  2:3429061/1‑150 (MQ=255)
                                            aGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGACGATGGCGGTGACCGGCACGTTGGCGAGCAGCATCAGGTGCGCCGGCa                                                                                                          >  2:4979172/1‑150 (MQ=255)
                                              aTCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCAcc                                                                                                        <  2:5662549/150‑1 (MQ=255)
                                              aTCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCAcc                                                                                                        <  2:8074341/150‑1 (MQ=255)
                                                    gCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTaa                                                                                                  >  1:1250705/1‑150 (MQ=255)
                                                    gCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACATGCCGTTGGCGGTGACCGGCCAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCAGGTaa                                                                                                  >  1:3673933/1‑150 (MQ=255)
                                                     ccAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAAt                                                                                                 <  1:8260628/150‑1 (MQ=255)
                                                        gTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCgg                                                                                              <  2:5602382/150‑1 (MQ=255)
                                                         ttGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCggg                                                                                             <  1:1652819/150‑1 (MQ=255)
                                                          tgctgcTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGc                                                                                            <  2:3066764/150‑1 (MQ=255)
                                                               ctCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCCCCATGTAATCGGGCAAGcc                                                                                       >  1:5377309/1‑150 (MQ=255)
                                                                tCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGTGCAGCATCAGGTGCGCCGGCACCATGTCATCGGGCAAGccc                                                                                      >  1:8054243/1‑150 (MQ=255)
                                                                tCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGccc                                                                                      <  1:2211143/150‑1 (MQ=255)
                                                                tCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGccc                                                                                      <  1:5844322/150‑1 (MQ=255)
                                                                 ccAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCt                                                                                     <  2:2437033/150‑1 (MQ=255)
                                                                  cAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTg                                                                                    >  1:2066954/1‑150 (MQ=255)
                                                                    gcgcCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGAc                                                                                  <  1:2636292/150‑1 (MQ=255)
                                                                    gcgcCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGAc                                                                                  <  2:865117/150‑1 (MQ=255)
                                                                       ccGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGTCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGa                                                                               <  1:6767926/150‑1 (MQ=255)
                                                                        cGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAg                                                                              <  1:5105735/150‑1 (MQ=255)
                                                                         gTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAgg                                                                             <  1:703516/150‑1 (MQ=255)
                                                                          tGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATAGGGCAAGCCCTGACCGAGGa                                                                            >  1:1149617/1‑150 (MQ=255)
                                                                             gcgGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAAt                                                                         >  1:648336/1‑150 (MQ=255)
                                                                              cgGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGTCAGGGCTTGCCCGattacatg                                                                        <  2:5104660/150‑8 (MQ=255)
                                                                                   cGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCa                                                                   <  2:8369067/150‑1 (MQ=255)
                                                                                       aTACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTg                                                               <  2:5256560/150‑1 (MQ=255)
                                                                                       aTACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTg                                                               >  2:8696465/1‑150 (MQ=255)
                                                                                       aTACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTg                                                               >  1:6124493/1‑150 (MQ=255)
                                                                                         acacCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCt                                                             <  2:6716802/150‑1 (MQ=255)
                                                                                           acCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCa                                                           >  1:2710097/1‑150 (MQ=255)
                                                                                           acCTGCCTGGCGGCGTTGAAATACGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCTGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCAAGTAATCGGGCAAGCCCTGACCGAGGAAATGAGCCAGTTGCTCa                                                           >  2:4766310/1‑150 (MQ=255)
                                                                                            cctgcctgGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCAc                                                          <  2:3931693/150‑1 (MQ=255)
                                                                                             ctgcctgGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACg                                                         <  2:6309154/150‑1 (MQ=255)
                                                                                             ctgcctgGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACATGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCAACAGGTGCGCCGGCACCATGTAATAGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCAcc                                                         >  1:5161158/1‑149 (MQ=255)
                                                                                                    gcggcgTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGc                                                  <  1:6891580/150‑1 (MQ=255)
                                                                                                     cggcgTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCt                                                 >  2:149706/1‑150 (MQ=255)
                                                                                                       gcgTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTgg                                               <  1:1386864/150‑1 (MQ=255)
                                                                                                        cgTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGc                                              <  1:3853886/150‑1 (MQ=255)
                                                                                                         gTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCCTGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGcc                                             >  2:7487300/1‑150 (MQ=255)
                                                                                                         gTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGcc                                             >  2:8183371/1‑150 (MQ=255)
                                                                                                          ttGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCCGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGccc                                            >  2:3554592/1‑150 (MQ=255)
                                                                                                          ttGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACATGCCGTTTGCGGTTAACGGCAACTTTGCGAGCAGCCTCAGGTGCGCCCGCAACAAGTAATAGGGCAAGCCCAGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTTGccc                                            >  1:1020256/1‑150 (MQ=255)
                                                                                                               aTCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTc                                       <  1:2986021/150‑1 (MQ=255)
                                                                                                                tCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCg                                      <  2:2549465/150‑1 (MQ=255)
                                                                                                                  cGGTGCCGGCAAGAGGTTGCGCGCCAACATGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCt                                    >  2:2231121/1‑150 (MQ=255)
                                                                                                                      gCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGccc                                <  2:7294449/150‑1 (MQ=255)
                                                                                                                          gCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACTCAGTCGCTGGCCCGCCTCGCTGCCCGGCt                            >  2:7036456/1‑150 (MQ=255)
                                                                                                                           cAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTc                           >  2:2318235/1‑150 (MQ=255)
                                                                                                                            aaGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCg                          <  2:895400/150‑1 (MQ=255)
                                                                                                                            aaGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCg                          >  2:319858/1‑150 (MQ=255)
                                                                                                                            aaGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCg                          >  2:3021041/1‑150 (MQ=255)
                                                                                                                              gagGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGc                        <  1:1840717/150‑1 (MQ=255)
                                                                                                                                 gTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTACTCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGcggg                     <  1:4020864/150‑4 (MQ=255)
                                                                                                                                 gTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGcggg                     <  1:3137656/150‑4 (MQ=255)
                                                                                                                                        gcCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCa              <  1:5477899/150‑1 (MQ=255)
                                                                                                                                        gcCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCa              <  2:5221099/150‑1 (MQ=255)
                                                                                                                                        gcCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCa              <  1:5615328/150‑1 (MQ=255)
                                                                                                                                        gcCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCa              <  2:422317/150‑1 (MQ=255)
                                                                                                                                        gcCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCa              <  1:732206/150‑1 (MQ=255)
                                                                                                                                        gcCAACTTGCCGGTCACCGCCAACGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCCGCTCGGCa              >  2:626124/1‑150 (MQ=38)
                                                                                                                                         ccAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAg             <  1:4623255/150‑1 (MQ=255)
                                                                                                                                         ccAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAg             >  1:876844/1‑150 (MQ=255)
                                                                                                                                                 ccGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGc     <  1:5081841/150‑1 (MQ=255)
                                                                                                                                                    ttGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGcacc  >  1:378937/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTGTCCTGCATACCGACTTGGGGCAGGTTCAGCACCGTCTCCCAGATCGACGCCAGTTGCTGCTCCAGCGCCGTGCGCGGAGCGACATACACCTGCCTGGCGGCGTTGAAATCCGGTGCCGGCAAGAGGTTGCGCGCCAACTTGCCGTTGGCGGTGACCGGCAAGTTGGCGAGCAGCATCAGGTGCGCCGGCACCATGTAATCGGGCAAGCCCTGACCGAGGAAATGCGCCAGTTGCTCACGCAGTCGCTGGCCCGCCTCGCTGCCCGGCTCGGCCTGCTCGGCAGCCATCGGCACC  >  NC_002947/4777958‑4778255

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: