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.99.7% 1529.1 / ‑2.6 311intergenic (+27/+16)mexK/PA14_RS06750multidrug efflux RND transporter permease subunit MexK/DUF4197 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (127/183);  total (128/183)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
*NC_0084631,441,8580C.99.7% 1528.4 / ‑2.7 311intergenic (+28/+15)mexK/PA14_RS06750multidrug efflux RND transporter permease subunit MexK/DUF4197 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (127/183);  minor base G (1/0);  total (128/183)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 357 total aligned reads displayed.

CTTCGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGCCCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGCCCTC  >  NC_008463/1441707‑1442004
                                                                                                                                                      ||                                                                                                                                                  
  tcGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccc                                                                                                                                                   >  2:5388743/1‑151 (MQ=255)
   cGGGCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGcccc                                                                                                                                                  <  1:7008580/151‑1 (MQ=255)
    gggCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccccg                                                                                                                                                 >  2:7690829/1‑150 (MQ=255)
    gggCCGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccccg                                                                                                                                                 <  1:6546301/151‑2 (MQ=255)
       ccGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccccgtcc                                                                                                                                              >  1:9573087/1‑147 (MQ=255)
       ccGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccccgtcc                                                                                                                                              >  1:8536334/1‑147 (MQ=255)
        cGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTCTGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccccgtccg                                                                                                                                             <  2:8640404/151‑6 (MQ=255)
        cGATGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCGccccgtccg                                                                                                                                             <  2:9854228/151‑6 (MQ=255)
           tGGCGGTGGCGATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGATCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCg                                                                                                                                        >  2:356791/1‑151 (MQ=255)
            ggCGGTGGCGATCATGGGCGGGCTGTTCGTCGCCACCGTGCTCACCCTGCGGTTACTGCCGGCGCGCTATGCCGCCTGGTTCCGTGGGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGA‑‑GCGCGCCCCGTCCGTCGa                                                                                                                                       >  2:10798268/1‑151 (MQ=255)
                    cGATCTTGGGCGGGCGGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTg                                                                                                                               <  1:3730137/151‑1 (MQ=255)
                     gATCATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTgg                                                                                                                              <  1:8705699/151‑1 (MQ=255)
                        cATGGGCGGGCTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCcg                                                                                                                           <  2:320813/151‑1 (MQ=255)
                           gggcgggcTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGTCGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCgcgc                                                                                                                        >  1:9763946/1‑151 (MQ=255)
                            ggcgggcTGATCGTCGCCACCGTGCACACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGcc                                                                                                                       <  1:2707942/151‑1 (MQ=255)
                             gcgggcTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCg                                                                                                                      <  1:9785488/151‑1 (MQ=255)
                              cgggcTGATCGTCGCCACCGTGCTCACCATGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGa                                                                                                                     >  2:2756535/1‑151 (MQ=255)
                               gggcTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGaa                                                                                                                    >  1:4167445/1‑151 (MQ=255)
                                ggcTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAg                                                                                                                   <  2:1921248/151‑1 (MQ=255)
                                 gcTGATCGTCGCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGa                                                                                                                  <  1:10693832/151‑1 (MQ=255)
                                        gtcgCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACcttctt                                                                                                           <  1:5158939/151‑1 (MQ=255)
                                        cccgCCACCGCGCTCACCCTGCCGTTCCTGCCGGCGCGCCCCGCCGCCTGGTTCCGTGTGAAGGCGGACGAGGCGCGCGCCTGACGGGTGAAGGGTACGGGGCGATCGCG‑‑CGCCGTCCGTCGACGCACCTGGCCGCGCCGAAGACcttctt                                                                                                           <  2:2368630/149‑1 (MQ=255)
                                         tcgCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGGGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTg                                                                                                          <  1:3451782/151‑1 (MQ=255)
                                          cgCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTgg                                                                                                         <  1:4262809/151‑1 (MQ=255)
                                          cgCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTgg                                                                                                         <  1:10494681/151‑1 (MQ=255)
                                          cgCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTgg                                                                                                         <  1:7063835/151‑1 (MQ=255)
                                           gCCACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGc                                                                                                        >  2:9589046/1‑151 (MQ=255)
                                            ccACCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGcc                                                                                                       <  2:7760534/151‑1 (MQ=255)
                                              aCCGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAg                                                                                                     <  2:171165/151‑1 (MQ=255)
                                               ccGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAgg                                                                                                    <  1:4492057/151‑1 (MQ=255)
                                                cGTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGc                                                                                                   >  1:4964613/1‑151 (MQ=255)
                                                 gTGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGcc                                                                                                  >  2:9206830/1‑151 (MQ=255)
                                                  tGCTCACCCTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCg                                                                                                 >  1:201472/1‑151 (MQ=255)
                                                        cccTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTggcg                                                                                           <  2:4988202/151‑1 (MQ=255)
                                                        cccTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTggcg                                                                                           >  2:3966070/1‑151 (MQ=255)
                                                        cccTGCTGTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTggcg                                                                                           >  2:8986586/1‑151 (MQ=255)
                                                          ctgctgTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTggcggc                                                                                         >  1:9274005/1‑151 (MQ=255)
                                                              tgTTCCTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGgccgcc                                                                                     <  1:4739396/151‑1 (MQ=255)
                                                                  ccTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGccgccgccg                                                                                 <  2:10684457/151‑1 (MQ=255)
                                                                  ccTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGccgccgccg                                                                                 <  2:3688753/151‑1 (MQ=255)
                                                                  ccTGCCGGCGCTCTATGCCGCCGGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGccgccgccg                                                                                 >  2:5819134/1‑151 (MQ=255)
                                                                   cTGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCgg                                                                                >  2:4674415/1‑151 (MQ=255)
                                                                    tGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGACGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCggg                                                                               <  1:5439002/151‑1 (MQ=255)
                                                                    tGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCggg                                                                               <  1:2105643/151‑1 (MQ=255)
                                                                    tGCCGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCggg                                                                               <  1:845296/151‑1 (MQ=255)
                                                                     gCCGGCGCTCTAGGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGt                                                                              <  2:1549377/151‑1 (MQ=255)
                                                                       cGGCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTc                                                                            >  2:7206337/1‑151 (MQ=255)
                                                                        ggCGCTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCt                                                                           >  1:1575753/1‑151 (MQ=255)
                                                                           gcTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGc                                                                        >  2:2998274/1‑151 (MQ=255)
                                                                           gcTCTATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGc                                                                        <  1:1783749/151‑1 (MQ=255)
                                                                            ctctATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCg                                                                       >  2:1201926/1‑151 (MQ=255)
                                                                             tctATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCgg                                                                      <  1:1237002/151‑1 (MQ=255)
                                                                              ctATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGa                                                                     >  1:1322878/1‑151 (MQ=255)
                                                                              ctATGCCGCCTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGa                                                                     <  1:5648686/151‑1 (MQ=255)
                                                                               tATGCCGCCTGGTGCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGAt                                                                    <  2:6622455/151‑1 (MQ=255)
                                                                                  gccgccTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGct                                                                 <  2:2144670/151‑1 (MQ=255)
                                                                                    cgccTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGACGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGctct                                                               <  1:8149162/151‑1 (MQ=255)
                                                                                    cgccTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGctct                                                               <  1:288525/151‑1 (MQ=255)
                                                                                     gccTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCtctt                                                              <  1:1005123/151‑1 (MQ=255)
                                                                                       cTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCTTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCtcttct                                                            <  1:9143527/151‑1 (MQ=255)
                                                                                       cTGGTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCCGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCtcttct                                                            >  2:4244452/1‑151 (MQ=255)
                                                                                         ggTTCCGTGTGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTcc                                                          >  1:6124816/1‑151 (MQ=255)
                                                                                               gtgtGAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTg                                                    <  2:2624702/151‑1 (MQ=255)
                                                                                                tgtgAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTgg                                                   <  1:3988694/151‑1 (MQ=255)
                                                                                                  tgAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCg                                                 <  1:2593119/151‑1 (MQ=255)
                                                                                                  tgAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCg                                                 <  1:3162108/151‑1 (MQ=255)
                                                                                                  tgAAGAAGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTGCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCGGGCGGATGCTCTTCTCCTGCTTGGCg                                                 >  2:8052636/1‑151 (MQ=255)
                                                                                                   gaagaaGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGCTCGCG‑‑CCCCGTCCGTCGACTCACCGGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCGGGCGGATGCTCTTCTCCTGCTTGGCGa                                                <  2:1319136/151‑1 (MQ=255)
                                                                                                   gaagaaGGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGa                                                >  2:6359165/1‑151 (MQ=255)
                                                                                                      gaagGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTTGCGGATGCTCTTCTCCTGCTTGGCGatca                                             >  1:9002741/1‑151 (MQ=255)
                                                                                                      gaagGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGatca                                             <  2:8890603/151‑1 (MQ=255)
                                                                                                       aagGACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGatcat                                            >  2:9508610/1‑151 (MQ=255)
                                                                                                          gACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCAtctc                                         >  1:4154521/1‑151 (MQ=255)
                                                                                                          gACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCAtctc                                         <  1:5792165/151‑1 (MQ=255)
                                                                                                          gACGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCAtctc                                         <  1:4875514/151‑1 (MQ=255)
                                                                                                           aCGAGGCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCg                                        >  1:368462/1‑151 (MQ=255)
                                                                                                               ggCGCGGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAaga                                    >  1:7584413/1‑151 (MQ=255)
                                                                                                                 cgcgGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGc                                  <  2:10603305/151‑1 (MQ=255)
                                                                                                                  gcgGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGTTGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGcc                                 <  2:7081217/151‑1 (MQ=255)
                                                                                                                   cgcgCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCCTCTCGAAGAGccc                                <  2:10503065/148‑1 (MQ=255)
                                                                                                                   cgGGCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCCGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGccc                                >  1:4566381/1‑151 (MQ=255)
                                                                                                                    gggCCTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCCCCTGGCCGCGCCGACGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCg                               <  2:2171443/151‑1 (MQ=255)
                                                                                                                       ccTGACGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTGGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCg                            <  2:5695973/151‑1 (MQ=255)
                                                                                                                           aCGGGAGACGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAggg                        <  2:8356940/151‑1 (MQ=255)
                                                                                                                           aCGGGAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAggg                        <  2:8503638/151‑1 (MQ=255)
                                                                                                                             gggAGAAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGcc                      >  2:2521669/1‑151 (MQ=255)
                                                                                                                               gagaAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTg                    <  2:6087436/151‑1 (MQ=255)
                                                                                                                                 gaAGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGtt                  <  1:7066597/151‑1 (MQ=255)
                                                                                                                                  aaGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCGGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTc                 <  2:6937506/151‑1 (MQ=255)
                                                                                                                                   aGGGGACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGACGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCg                <  2:4007760/151‑1 (MQ=255)
                                                                                                                                       gACGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGa            <  1:6246399/151‑1 (MQ=255)
                                                                                                                                         cGGGGCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACa          <  1:10176207/151‑1 (MQ=255)
                                                                                                                                            ggCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAg       <  2:5241381/151‑1 (MQ=255)
                                                                                                                                            ggCGATCGCG‑‑CCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAg       <  2:6543710/151‑1 (MQ=255)
                                                                                                                                               gcgatcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGc      <  1:6775725/146‑1 (MQ=255)
                                                                                                                                                cgatcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGcc     <  1:2981452/147‑1 (MQ=255)
                                                                                                                                                 gatcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAATACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGccc    <  1:772108/148‑1 (MQ=255)
                                                                                                                                                 gatcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGccc    >  1:8139712/4‑151 (MQ=255)
                                                                                                                                                 gatcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGccc    <  1:8871571/148‑1 (MQ=255)
                                                                                                                                                   tcgcgCCCCGTCCGTCGACTCACCTGGCCGCGCCGAAGACCTTCTTGGCCAGGCCGGTGGCGGCCGCCGCCGGGTTCTGGCGGATGCTCTTCTCCTGCTTGGCGATCATCTCGAAGAGCCCGTCGAGGGCCTGTTCGGTGACATAGCCCTc  <  1:2320421/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: