Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 1,441,861 (C)6→4 intergenic (+31/+11) mexK → / ← PA14_RS06750 multidrug efflux RND transporter permease subunit MexK/DUF4197 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084631,441,8570C.98.9% 1775.3 / 10.0 366intergenic (+27/+16)mexK/PA14_RS06750multidrug efflux RND transporter permease subunit MexK/DUF4197 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base . (173/189);  total (174/192)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0084631,441,8580C.99.2% 1778.8 / 6.5 366intergenic (+28/+15)mexK/PA14_RS06750multidrug efflux RND transporter permease subunit MexK/DUF4197 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base . (173/189);  minor base G (1/2);  total (174/192)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 420 total aligned reads displayed.

CTTCGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGCCCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGCCCTC  >  NC_008463/1441707‑1442004
                                                                                                                                                      ||                                                                                                                                                  
  tcGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccc                                                                                                                                                   >  1:1370878/1‑151 (MQ=255)
  tcGGGCCGAGGGCGGTGGCGATCATGGGCGGGCGGATCGGCGGCACGGTGCTCTCCCTGCTGTTCCTGCCGGCGCTCTGTGCCGCCTGGTTGGGTGTGATGATGGACGAGGGGCGGTCCTGGCGGGCGGTGGGGACGGGGCGATCGCGccc                                                                                                                                                   >  1:2407003/1‑151 (MQ=255)
      gCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccccgtc                                                                                                                                               >  2:9624847/1‑148 (MQ=255)
       ccGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGGCGGGAGAAGGGGACGGGGCGATCGCGccccgtcc                                                                                                                                              >  1:7687104/1‑147 (MQ=255)
       ccGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccccgtcc                                                                                                                                              <  2:3976857/151‑5 (MQ=255)
           tGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCg                                                                                                                                        >  2:1594828/1‑151 (MQ=255)
               ggTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTc                                                                                                                                    >  2:10627906/1‑151 (MQ=255)
               ggTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTc                                                                                                                                    <  1:3705059/151‑1 (MQ=255)
                gTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCa                                                                                                                                   >  1:3455619/1‑151 (MQ=255)
                 tGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGATGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCAc                                                                                                                                  >  1:828635/1‑151 (MQ=255)
                    cGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTg                                                                                                                               >  2:6386056/1‑151 (MQ=255)
                    cGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTg                                                                                                                               >  2:2572604/1‑151 (MQ=255)
                     gATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTgg                                                                                                                              >  2:773534/1‑151 (MQ=255)
                      aTCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGTACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGc                                                                                                                             <  1:5560730/151‑1 (MQ=255)
                      aTCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGc                                                                                                                             >  1:9630216/1‑151 (MQ=255)
                        cATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCcg                                                                                                                           <  2:8053029/151‑1 (MQ=255)
                        cATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCcg                                                                                                                           >  2:7221424/1‑151 (MQ=255)
                          tGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCcgcg                                                                                                                         <  1:1872879/151‑1 (MQ=255)
                            ggcgggcTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGcc                                                                                                                       <  1:3202817/151‑1 (MQ=255)
                            ggcgggcTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGcc                                                                                                                       <  1:6088441/151‑1 (MQ=255)
                              cgggcTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGa                                                                                                                     >  1:4874463/1‑151 (MQ=255)
                                  cTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGAc                                                                                                                 >  1:9783961/1‑151 (MQ=255)
                                   tGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGAcc                                                                                                                >  2:7920069/1‑151 (MQ=255)
                                        gtcgCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACcttctt                                                                                                           >  1:4560949/1‑151 (MQ=255)
                                          cgCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTgg                                                                                                         <  2:9078732/151‑1 (MQ=255)
                                           gCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGc                                                                                                        <  2:1452292/151‑1 (MQ=255)
                                             cACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCa                                                                                                      >  2:8975645/1‑151 (MQ=255)
                                              aCCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAg                                                                                                     <  1:7063120/151‑1 (MQ=255)
                                              aCCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTGCTTGGCCAg                                                                                                     >  1:10276264/1‑151 (MQ=255)
                                               ccGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAgg                                                                                                    >  1:9582402/1‑151 (MQ=255)
                                                 gTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGcc                                                                                                  >  2:5191233/1‑151 (MQ=255)
                                                 gTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGcc                                                                                                  >  2:7258989/1‑151 (MQ=255)
                                                   gCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCgg                                                                                                >  1:9873316/1‑151 (MQ=255)
                                                   gCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTACTTGGCCAGGCCgg                                                                                                >  1:5059667/1‑151 (MQ=255)
                                                    cTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGCGACGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCCGGCCGGt                                                                                               >  2:4550580/1‑151 (MQ=255)
                                                    cTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGGGTGACGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGACGACCTTCTTGGCCAGGCCGGt                                                                                               >  1:8234003/1‑151 (MQ=255)
                                                      cACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGACGACCTTCTTGGCCAGGCCGGTgg                                                                                             >  2:10126214/1‑151 (MQ=255)
                                                              tgTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGgccgcc                                                                                     >  2:5689548/1‑151 (MQ=255)
                                                               gTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGgccgccg                                                                                    >  2:3085368/1‑151 (MQ=255)
                                                                 tCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGgccgccgcc                                                                                  <  1:629307/151‑1 (MQ=255)
                                                                 tCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGgccgccgcc                                                                                  <  2:187076/151‑1 (MQ=255)
                                                                    tGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCggg                                                                               <  1:10713478/151‑1 (MQ=255)
                                                                    tGCCGGCGCTATATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCggg                                                                               <  1:2059068/151‑1 (MQ=255)
                                                                       cGGCGCTCTAGGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTc                                                                            <  1:9278199/151‑1 (MQ=255)
                                                                        ggCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCt                                                                           <  1:6911387/151‑1 (MQ=255)
                                                                         gcgcTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTg                                                                          >  2:3701245/1‑151 (MQ=255)
                                                                          cgcTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTgg                                                                         >  2:7627054/1‑151 (MQ=255)
                                                                           gcTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGc                                                                        >  2:10375732/1‑151 (MQ=255)
                                                                            ctctATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCg                                                                       <  2:8371823/151‑1 (MQ=255)
                                                                             tctATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCgg                                                                      >  2:9396228/1‑151 (MQ=255)
                                                                              ctATGCCGCCTGGTTCCTTGTGAAGAAGTACGAGTCGCGGGCCTGACGGGAGAAGGGGACGGGGCGAACGCG‑‑CCCCGACCGTCGACTCACCTGGACGCGCCGAAGACCTTCTTGGCCAGGCCGGTTGCGGCCGCCGCCGGGTTCTGGCGGa                                                                     >  1:1093171/1‑151 (MQ=255)
                                                                               tATGCCGCCTGGTTCCGTGTGACGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGAt                                                                    <  1:4643023/151‑1 (MQ=255)
                                                                                      ccTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCtcttc                                                             <  1:2944398/151‑1 (MQ=255)
                                                                                      ccTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCtcttc                                                             >  1:539683/1‑151 (MQ=255)
                                                                                       cTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCtcttct                                                            >  2:9181521/1‑151 (MQ=255)
                                                                                       cTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCtcttct                                                            >  2:560176/1‑151 (MQ=255)
                                                                                         ggTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGACGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGCCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTcc                                                          <  1:10031209/151‑1 (MQ=255)
                                                                                           ttCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTg                                                        <  2:2764790/151‑1 (MQ=255)
                                                                                             ccGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCt                                                      <  1:8896515/151‑1 (MQ=255)
                                                                                              cGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCtt                                                     <  2:380186/151‑1 (MQ=255)
                                                                                               gtgtGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTg                                                    <  1:9363516/151‑1 (MQ=255)
                                                                                               gtgtGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTg                                                    >  1:3869082/1‑151 (MQ=255)
                                                                                                 gtgAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCCGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGc                                                  >  2:5914493/1‑151 (MQ=255)
                                                                                                 gtgAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGc                                                  >  1:8067476/1‑151 (MQ=255)
                                                                                                     agaagGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGatc                                              <  1:5223980/151‑1 (MQ=255)
                                                                                                     agaagGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGatc                                              <  2:8488454/151‑1 (MQ=255)
                                                                                                       aagGACGCGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGatcat                                            <  1:10429122/151‑1 (MQ=255)
                                                                                                        agGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGatcatc                                           >  1:7486670/1‑151 (MQ=255)
                                                                                                         ggACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCAtct                                          <  2:6610016/151‑1 (MQ=255)
                                                                                                         ggACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCAtct                                          >  1:1708361/1‑151 (MQ=255)
                                                                                                             gAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCCGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGaa                                      >  1:8169059/1‑151 (MQ=255)
                                                                                                             gAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGaa                                      <  1:8749638/151‑1 (MQ=255)
                                                                                                             gAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGaa                                      >  1:9007496/1‑151 (MQ=255)
                                                                                                              cggCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGGGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGaac                                     <  1:1997392/150‑2 (MQ=255)
                                                                                                              aGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAag                                     <  2:2951660/151‑1 (MQ=255)
                                                                                                                  gcgGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGcc                                 >  2:1257961/1‑151 (MQ=255)
                                                                                                                    gggCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCg                               >  2:10841229/1‑151 (MQ=255)
                                                                                                                    gggCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCg                               >  1:1801579/1‑151 (MQ=255)
                                                                                                                         tGACGGGAGAAGGGGACGGGGCGCTCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAg                          <  2:3464154/151‑1 (MQ=255)
                                                                                                                         tGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGTCCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAg                          >  2:1056711/1‑151 (MQ=255)
                                                                                                                            cGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTACTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGc                       >  1:423740/1‑151 (MQ=255)
                                                                                                                                agaAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGt                   >  1:6693662/1‑151 (MQ=255)
                                                                                                                                agaAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGt                   >  2:2229088/1‑151 (MQ=255)
                                                                                                                                agaAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCGGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGt                   <  1:7185074/151‑1 (MQ=255)
                                                                                                                                 gaAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGtt                  <  1:8425055/151‑1 (MQ=255)
                                                                                                                                 gaAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGtt                  >  1:6192603/1‑151 (MQ=255)
                                                                                                                                   aGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCg                <  1:4335652/151‑1 (MQ=255)
                                                                                                                                    ggggAGGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCgg               <  1:15052/151‑1 (MQ=255)
                                                                                                                                     gggACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCAGTCGAGGGCCTGTTCGGt              <  2:4954332/151‑1 (MQ=255)
                                                                                                                                          ggggCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTGCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACAt         >  2:7411027/1‑151 (MQ=255)
                                                                                                                                           gggCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATa        >  2:6946438/1‑151 (MQ=255)
                                                                                                                                           gggCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATa        <  2:3543199/151‑1 (MQ=255)
                                                                                                                                            ggCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAg       <  2:5798075/151‑1 (MQ=255)
                                                                                                                                               gcgatcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGc      <  1:2572604/146‑1 (MQ=255)
                                                                                                                                                cgatcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGcc     <  2:8278636/147‑1 (MQ=255)
                                                                                                                                                 gatcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGccc    >  1:5619133/4‑151 (MQ=255)
                                                                                                                                                  atcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGCCCt   <  2:2377285/149‑1 (MQ=255)
                                                                                                                                                  atcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGAGCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGCCCt   >  2:10034818/3‑151 (MQ=255)
                                                                                                                                                   tcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGCCCTc  <  2:6130472/150‑1 (MQ=255)
                                                                                                                                                   tcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGCCCTc  <  1:2190425/150‑1 (MQ=255)
                                                                                                                                                      ||                                                                                                                                                  
CTTCGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGCCCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGCCCTC  >  NC_008463/1441707‑1442004

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: