Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,787,777 C→G pseudogene (855/3090 nt) PA14_RS13085 → CusA/CzcA family heavy metal efflux RND transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084632,787,7770CG99.7% 1527.2 / NA 390pseudogene (855/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (187/202);  minor base T (1/0);  total (188/202)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 399 total aligned reads displayed.

TCGACGGCGTGCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCG  >  NC_008463/2787629‑2787924
                                                                                                                                                    |                                                                                                                                                      
tCGACGGCGTGCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAggc                                                                                                                                                      >  2:10022893/1‑151 (MQ=255)
 cGACGGCGTGCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAggcg                                                                                                                                                     >  1:2945028/1‑151 (MQ=255)
 cGACGGCGTGCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAggcg                                                                                                                                                     <  1:8864275/151‑1 (MQ=255)
  gACGGCGTGCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCgg                                                                                                                                                    >  1:4755933/1‑151 (MQ=255)
        gTGCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGcc                                                                                                                                              >  2:7308152/1‑151 (MQ=255)
         tGCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCa                                                                                                                                             <  1:11222598/151‑1 (MQ=255)
         tGCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCa                                                                                                                                             <  2:8000845/151‑1 (MQ=255)
          gCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGACGTGGCCCAGGCGGTCGGCCAg                                                                                                                                            >  2:729112/1‑151 (MQ=255)
          gCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAg                                                                                                                                            >  2:1337233/1‑151 (MQ=255)
           ccGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGACCGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAgc                                                                                                                                           <  2:2116314/151‑1 (MQ=255)
            cGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAgcg                                                                                                                                          >  1:2054788/1‑151 (MQ=255)
            cGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAgcg                                                                                                                                          <  2:63146/151‑1 (MQ=255)
                ccGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTg                                                                                                                                      <  1:1580454/151‑1 (MQ=255)
                  gTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTgga                                                                                                                                    <  1:5375382/151‑1 (MQ=255)
                  gTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTgga                                                                                                                                    >  1:3556512/1‑151 (MQ=255)
                   tATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTggag                                                                                                                                   <  1:7201471/151‑1 (MQ=255)
                    aTCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTggagg                                                                                                                                  <  1:3788825/151‑1 (MQ=255)
                         cgACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATc                                                                                                                             <  2:11393449/151‑1 (MQ=255)
                          gACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCa                                                                                                                            <  2:951622/151‑1 (MQ=255)
                            cGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAAc                                                                                                                          <  1:8582622/151‑1 (MQ=255)
                             gTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAAcc                                                                                                                         <  2:1447309/151‑1 (MQ=255)
                              tCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCg                                                                                                                        <  2:1169837/151‑1 (MQ=255)
                               cGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCgg                                                                                                                       <  2:3167417/151‑1 (MQ=255)
                                  cGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGAcc                                                                                                                    <  2:4132228/151‑1 (MQ=255)
                                    aGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCt                                                                                                                  <  2:11037675/151‑1 (MQ=255)
                                    aGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCt                                                                                                                  <  2:3402272/151‑1 (MQ=255)
                                     ggTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCGCTGGCGGTCGGCGAGCGCCTGGAGGTGATCAACCGGTCCCTg                                                                                                                 <  2:722098/151‑1 (MQ=255)
                                        cGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGccc                                                                                                             <  2:11492076/151‑1 (MQ=255)
                                        cGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCAGCGAAGTTGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCAAGCGCCTGGAGGAGATCAACCGGACCCTGccc                                                                                                             <  1:10869122/151‑1 (MQ=255)
                                         ggCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCa                                                                                                           >  2:2897681/1‑151 (MQ=255)
                                           ccTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCaag                                                                                                         >  2:5582575/1‑149 (MQ=255)
                                            cTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAagg                                                                                                        >  1:5601488/1‑151 (MQ=255)
                                               ggCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAaggggg                                                                                                     <  2:3279074/151‑2 (MQ=255)
                                                gCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAagggggt                                                                                                    <  1:4947643/151‑3 (MQ=255)
                                                  aaGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCa                                                                                                   >  1:7528017/1‑151 (MQ=255)
                                                  aaGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCa                                                                                                   >  1:5052733/1‑151 (MQ=255)
                                                  aaGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCa                                                                                                   >  1:269754/1‑151 (MQ=255)
                                                  aaGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGAGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCAGCACAAGGGGGTCa                                                                                                   >  2:7475060/1‑151 (MQ=255)
                                                     gAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCaagg                                                                                                >  1:6276639/1‑148 (MQ=255)
                                                      aGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGAGCAACCGGACCCTGCCCAAGGGGGTCaaggc                                                                                               >  2:4526035/1‑147 (MQ=255)
                                                         tGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGAt                                                                                            <  2:10436308/151‑1 (MQ=255)
                                                          gcgcACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATc                                                                                           <  1:1676545/151‑1 (MQ=255)
                                                              aCCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCg                                                                                       <  1:6737041/151‑1 (MQ=255)
                                                              aCCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCg                                                                                       <  1:9620124/151‑1 (MQ=255)
                                                               ccGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGt                                                                                      <  2:7688837/151‑1 (MQ=255)
                                                                   cgcgGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTAc                                                                                  >  1:9112838/1‑151 (MQ=255)
                                                                   cgcgGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTAc                                                                                  <  2:9121151/151‑1 (MQ=255)
                                                                     cgGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCACCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGa                                                                                <  1:9837698/151‑1 (MQ=255)
                                                                          aCCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGca                                                                           >  2:9685985/1‑151 (MQ=255)
                                                                            cGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGcacc                                                                         >  2:9010109/1‑151 (MQ=255)
                                                                            cGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGcacc                                                                         <  1:2413645/151‑1 (MQ=255)
                                                                             gagaACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGcacca                                                                        >  2:2680152/1‑151 (MQ=255)
                                                                               gaACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCaccacc                                                                      >  2:8862717/1‑151 (MQ=255)
                                                                                aaCGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCAccc                                                                     <  1:3341470/151‑1 (MQ=255)
                                                                                  cGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTg                                                                   <  2:9615177/151‑1 (MQ=255)
                                                                                  cGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTg                                                                   <  2:4432445/151‑1 (MQ=255)
                                                                                   ggCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTgg                                                                  <  1:5865546/151‑2 (MQ=255)
                                                                                   ggCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTgg                                                                  <  2:6583509/151‑2 (MQ=255)
                                                                                     ccGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtg                                                                <  2:8585962/151‑4 (MQ=255)
                                                                                      cgcgAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtgg                                                               <  1:6231617/151‑5 (MQ=255)
                                                                                        cgAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCTGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtggac                                                             <  2:611285/151‑7 (MQ=255)
                                                                                        cgAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtggac                                                             <  2:5135877/151‑7 (MQ=255)
                                                                                        cgAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtggac                                                             <  1:7607787/151‑7 (MQ=255)
                                                                                        cgAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtggac                                                             >  2:8453568/1‑145 (MQ=255)
                                                                                        cgAAGTGGTGCACGGCACGGTGTTCATGCTGCTCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGAACAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACCACCGCACCACCCTggtggac                                                             <  2:6921778/151‑7 (MQ=255)
                                                                                                                 aTGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaaga                                   <  1:6495191/151‑1 (MQ=255)
                                                                                                                   gCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAc                                 <  1:7732356/151‑1 (MQ=255)
                                                                                                                     tGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcct                               >  1:8444812/1‑150 (MQ=255)
                                                                                                                       aTCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctgg                             >  2:10237874/1‑148 (MQ=255)
                                                                                                                       aTCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctgg                             >  1:1276348/1‑148 (MQ=255)
                                                                                                                       aTCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGGCAAGGCGATCACCGTCTACGACCGCACCCCCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctgg                             <  2:8158551/151‑4 (MQ=255)
                                                                                                                        tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctggt                            <  2:7267939/151‑5 (MQ=255)
                                                                                                                        tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctggt                            >  2:3844419/1‑147 (MQ=255)
                                                                                                                          ggCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCtggtgg                         >  1:1047154/1‑151 (MQ=255)
                                                                                                                          ggCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCtggtgg                         >  2:4970934/1‑151 (MQ=255)
                                                                                                                          ggCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGAACGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCtggtgg                         >  2:3579120/1‑151 (MQ=255)
                                                                                                                            cGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGGGGCCACGGTGAAGAAGAACCGGGTGGaa                       >  1:11363815/1‑151 (MQ=255)
                                                                                                                            cGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGaa                       <  1:10021710/151‑1 (MQ=255)
                                                                                                                              agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAgg                     <  2:488980/151‑1 (MQ=255)
                                                                                                                              agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAgg                     >  2:6673776/1‑151 (MQ=255)
                                                                                                                               gaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggc                    >  2:6715355/1‑151 (MQ=255)
                                                                                                                                  cAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcg                 <  1:2568470/151‑1 (MQ=255)
                                                                                                                                  cAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcg                 <  1:3202179/151‑1 (MQ=255)
                                                                                                                                   aGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcgg                <  1:4417183/151‑1 (MQ=255)
                                                                                                                                    gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGCAGAAGAACCTGGTGGAAggcgcggc               >  1:5270435/1‑151 (MQ=255)
                                                                                                                                    gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggc               <  2:10738038/151‑1 (MQ=255)
                                                                                                                                    gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggc               <  2:8255732/151‑1 (MQ=255)
                                                                                                                                     ccGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggcg              <  1:7968473/151‑1 (MQ=255)
                                                                                                                                     ccGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggcg              >  1:10646213/1‑151 (MQ=255)
                                                                                                                                      cGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggcgc             <  1:7512225/151‑1 (MQ=255)
                                                                                                                                        ggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTg           <  1:5442387/151‑1 (MQ=255)
                                                                                                                                         gAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCCCCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTgg          >  1:6340160/1‑151 (MQ=255)
                                                                                                                                           aGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTg        <  1:11115715/151‑1 (MQ=255)
                                                                                                                                            gTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCTGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGa       <  2:9376371/151‑1 (MQ=255)
                                                                                                                                              ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATc     >  1:9291483/1‑151 (MQ=255)
                                                                                                                                              ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATc     <  2:2374468/151‑1 (MQ=255)
                                                                                                                                              ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATc     <  2:3073676/151‑1 (MQ=255)
                                                                                                                                               gCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcg    >  2:10874877/1‑151 (MQ=255)
                                                                                                                                                cccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcgc   <  2:5857542/151‑1 (MQ=255)
                                                                                                                                                 ccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcgcg  <  1:2917471/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TCGACGGCGTGCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCG  >  NC_008463/2787629‑2787924

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