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.7% 1481.5 / 9.3 308intergenic (+27/+16)mexK/PA14_RS06750multidrug efflux RND transporter permease subunit MexK/DUF4197 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base . (155/149);  total (159/149)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0084631,441,8580C.99.0% 1483.5 / 7.5 308intergenic (+28/+15)mexK/PA14_RS06750multidrug efflux RND transporter permease subunit MexK/DUF4197 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base . (155/149);  minor base G (3/0);  total (159/149)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 348 total aligned reads displayed.

TTCGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGCCCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGCCCTC  >  NC_008463/1441708‑1442004
                                                                                                                                                     ||                                                                                                                                                  
 tcGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGcc                                                                                                                                                    >  1:4975865/1‑150 (MQ=255)
   gggCCGTTGGCGTTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCGGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGACGGGGACGGGGCGATCGCGcccc                                                                                                                                                  <  2:8408249/150‑1 (MQ=255)
   gggCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGcccc                                                                                                                                                  >  1:8318275/1‑150 (MQ=255)
   gggCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGcccc                                                                                                                                                  <  1:10794338/150‑1 (MQ=255)
   gggCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGcccc                                                                                                                                                  <  2:7908463/150‑1 (MQ=255)
    ggCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccccg                                                                                                                                                 >  2:4513743/1‑149 (MQ=255)
       cGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccccgtcc                                                                                                                                              <  2:5790442/150‑5 (MQ=255)
        gATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccccgtccg                                                                                                                                             <  2:8808681/150‑6 (MQ=255)
         aTGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccccgtccgt                                                                                                                                            >  2:5150820/1‑144 (MQ=255)
             cGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGAc                                                                                                                                      <  1:6580873/150‑1 (MQ=255)
               gTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTc                                                                                                                                    <  1:5660390/150‑1 (MQ=255)
               gTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTc                                                                                                                                    >  2:3411746/1‑150 (MQ=255)
               gTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTc                                                                                                                                    <  2:4856720/150‑1 (MQ=255)
               gTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTACCGTGTGAAGAAGGACGTGGCGCGGGCCTGCCGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTc                                                                                                                                    <  1:9452212/150‑1 (MQ=255)
                 ggCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCAc                                                                                                                                  <  1:3535223/150‑1 (MQ=255)
                    gATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTg                                                                                                                               >  2:8150319/1‑150 (MQ=255)
                     agcctGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGGTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTgg                                                                                                                              <  1:5156352/146‑1 (MQ=255)
                         tGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCcgc                                                                                                                          >  2:6148596/1‑150 (MQ=255)
                            gcgggcTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGcc                                                                                                                       >  2:1696603/1‑150 (MQ=255)
                               ggcTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGaa                                                                                                                    >  2:5939252/1‑150 (MQ=255)
                                 cTGATCGTCGCCACCGGGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGa                                                                                                                  <  1:4208270/150‑1 (MQ=255)
                                    aTCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACct                                                                                                               <  2:2544572/150‑1 (MQ=255)
                                        tcgCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACcttctt                                                                                                           <  1:10244024/150‑1 (MQ=255)
                                        tcgCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACcttctt                                                                                                           >  1:10869701/1‑150 (MQ=255)
                                        tcgCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACcttctt                                                                                                           >  1:1676061/1‑150 (MQ=255)
                                        tcgCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACcttctt                                                                                                           >  1:2507413/1‑150 (MQ=255)
                                        tcgCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACcttctt                                                                                                           >  1:3377712/1‑150 (MQ=255)
                                          gCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTgg                                                                                                         >  1:9147847/1‑150 (MQ=255)
                                           ccACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGc                                                                                                        >  2:9225489/1‑150 (MQ=255)
                                           ccACCGCGCACACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGCGTGAAGAAGGACGCGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGc                                                                                                        <  1:1790284/150‑1 (MQ=255)
                                             aCCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCa                                                                                                      >  2:9779125/1‑150 (MQ=255)
                                              ccGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAg                                                                                                     >  2:10001476/1‑150 (MQ=255)
                                              ccGTGCTCACACGGCTGTTCCTGCCGGCGCTCAATGCCGACTGGTTCCGTGAGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCTAACGCG‑‑CACCGTCCGTCGACTCACCTGGCCGCGCCGAATCCCATCTTGGCCAg                                                                                                     >  2:6322323/1‑150 (MQ=255)
                                               cGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAgg                                                                                                    <  1:1448759/150‑1 (MQ=255)
                                               attGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGAGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAgg                                                                                                    >  2:1464368/3‑150 (MQ=255)
                                                 tGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTGCTTGGCCAGGcc                                                                                                  >  2:98542/1‑150 (MQ=255)
                                                  gCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCg                                                                                                 >  2:9307884/1‑150 (MQ=255)
                                                  gCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCg                                                                                                 <  1:1281827/150‑1 (MQ=255)
                                                    tCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGt                                                                                               <  1:800118/150‑1 (MQ=255)
                                                    tCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGt                                                                                               >  2:4159091/1‑150 (MQ=255)
                                                    tCACCCTGCTGGTCCTGCCGGCGCTCTATGCCGCGTGGTTCCGTGTGAAGAAGGACTAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGTTCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGt                                                                                               <  1:785860/150‑1 (MQ=255)
                                                          tgctgTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGACGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTggcggc                                                                                         <  2:2662501/150‑1 (MQ=255)
                                                          tgctgTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTggcggc                                                                                         >  2:10401661/1‑150 (MQ=255)
                                                           gctgTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGgcc                                                                                        <  2:2849968/150‑1 (MQ=255)
                                                            ctgTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGgccg                                                                                       <  2:9518559/150‑1 (MQ=255)
                                                             tgTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGgccgc                                                                                      <  2:10517299/150‑1 (MQ=255)
                                                              gTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGgccgcc                                                                                     <  1:10461818/150‑1 (MQ=255)
                                                               ttCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCGTCTTGGCCAGGCCGGTGGCGgccgccg                                                                                    >  1:9249478/1‑150 (MQ=255)
                                                                  cTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGccgccgccg                                                                                 <  1:4865917/150‑1 (MQ=255)
                                                                          gcTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTgg                                                                         <  2:9083060/150‑1 (MQ=255)
                                                                          gcTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTgg                                                                         >  1:6741813/1‑150 (MQ=255)
                                                                          gcTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTgg                                                                         >  1:2938235/1‑150 (MQ=255)
                                                                          gcTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCGTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTgg                                                                         >  2:8554876/1‑150 (MQ=255)
                                                                            tctATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCTATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCg                                                                       <  1:4055007/150‑1 (MQ=255)
                                                                            tctATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCg                                                                       <  2:1203307/150‑1 (MQ=255)
                                                                             ctATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCgg                                                                      >  2:11149274/1‑150 (MQ=255)
                                                                              tATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGa                                                                     <  1:5811758/150‑1 (MQ=255)
                                                                              tATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGa                                                                     >  2:10544791/1‑150 (MQ=255)
                                                                               aTGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGAt                                                                    <  2:2119441/150‑1 (MQ=255)
                                                                                 gccgccTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGc                                                                  >  1:6782187/1‑150 (MQ=255)
                                                                                 gccgccTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGc                                                                  >  1:6379057/1‑150 (MQ=255)
                                                                                   cgccTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACGCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGctc                                                                <  1:720054/150‑1 (MQ=255)
                                                                                    gccTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGctct                                                               >  1:5949799/1‑150 (MQ=255)
                                                                                      cTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCtcttc                                                             >  2:5044906/1‑150 (MQ=255)
                                                                                       tGGTTCCGTGTGAAGAAGGACGAGGCGCGGGACTGACGGGAGAAGGGTACGGGGCGATCGCG‑‑CCCAGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCAAGGCCGGTGGCGGCCGCCGCAGGTTTCTGGCGGATGCtcttct                                                            >  2:134095/1‑150 (MQ=255)
                                                                                         gTTCCTTGTTAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGACCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTcc                                                          <  1:156883/150‑1 (MQ=255)
                                                                                          ttCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGGGAAGGGGACGGGGCGATCGCG‑‑CGCCGTCCGTCGACTCTGCTGGCCGCGCCGATGTCGGTGTTGGCCAGGCCGGTGGCTGACGGCGCGGGGTTCTCTCGGATGATGTtctcag                                                         >  2:6783943/1‑148 (MQ=255)
                                                                                            ccGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGc                                                       <  1:5461478/150‑1 (MQ=255)
                                                                                              gtgtGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCtt                                                     >  1:8499948/1‑150 (MQ=255)
                                                                                                gtgAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTgg                                                   <  1:6443550/150‑1 (MQ=255)
                                                                                                gtgAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTgg                                                   <  1:6175594/150‑1 (MQ=255)
                                                                                                 tgAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGc                                                  <  1:2736405/150‑1 (MQ=255)
                                                                                                        ggACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCCCCGGGCCGCGCCGACGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGatcatc                                           <  2:568909/150‑1 (MQ=255)
                                                                                                         gACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCAtct                                          <  2:2939964/150‑1 (MQ=255)
                                                                                                            gAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGACGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGa                                       <  2:7274976/150‑1 (MQ=255)
                                                                                                             aGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGaa                                      >  2:7754274/1‑150 (MQ=255)
                                                                                                             aGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGaa                                      >  2:4307498/1‑150 (MQ=255)
                                                                                                                 gcgGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGc                                  >  2:8390784/1‑150 (MQ=255)
                                                                                                                  cgGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGcc                                 >  2:830142/1‑150 (MQ=255)
                                                                                                                  cgGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGcc                                 <  1:4402852/150‑1 (MQ=255)
                                                                                                                     gCCTGACGGTAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGt                              <  2:8165770/150‑1 (MQ=255)
                                                                                                                       cTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCg                            <  2:9475698/150‑1 (MQ=255)
                                                                                                                          aCGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAgg                         <  1:3165612/150‑1 (MQ=255)
                                                                                                                            gggAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGc                       >  2:7092407/1‑150 (MQ=255)
                                                                                                                               agaAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTg                    >  1:10041014/1‑150 (MQ=255)
                                                                                                                                   ggggACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCg                >  2:6413077/1‑150 (MQ=255)
                                                                                                                                    gggACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGTATTCTCTTCTCCTGCTTTGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCgg               <  1:10324529/150‑1 (MQ=255)
                                                                                                                                    gggACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCgg               >  1:301704/1‑150 (MQ=255)
                                                                                                                                      gACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGACGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTg             <  2:5386841/150‑1 (MQ=255)
                                                                                                                                       aCGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGGTCTGGCGGATGCTCTTCTCCTGCTTGGCGCTCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGa            >  1:7529794/1‑150 (MQ=255)
                                                                                                                                           ggCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATa        >  2:5564096/1‑150 (MQ=255)
                                                                                                                                           ggCGATCGCG‑‑CCCCGTCCGTCGACACACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATa        <  1:9786241/150‑1 (MQ=255)
                                                                                                                                              gcgatcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAg       >  1:8777462/6‑150 (MQ=255)
                                                                                                                                              gcgatcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAg       <  1:7623665/145‑1 (MQ=255)
                                                                                                                                               cgatcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGc      <  1:5335224/146‑1 (MQ=255)
                                                                                                                                                 atcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGccc    >  1:8974384/3‑150 (MQ=255)
                                                                                                                                                 atcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGccc    <  1:8105135/148‑1 (MQ=255)
                                                                                                                                                 agcgcgcCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGccc    <  1:2248084/149‑1 (MQ=255)
                                                                                                                                                   cgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGCCCTc  <  2:319546/150‑1 (MQ=255)
                                                                                                                                                   cgcgCCCCGACCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGCCCTc  <  1:7156/150‑1 (MQ=255)
                                                                                                                                                     ||                                                                                                                                                  
TTCGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGCCCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGCCCTC  >  NC_008463/1441708‑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: