Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,787,820 +A pseudogene (898/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,8191.A93.2% 1473.8 / 92.8 396pseudogene (897/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base . (8/19);  new base A (170/199);  total (178/218)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 404 total aligned reads displayed.

CGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCA  >  NC_008463/2787669‑2787967
                                                                                                                                                       |                                                                                                                                                     
 ggCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCa                                                                                                                                                      <  1:9427059/151‑1 (MQ=255)
    cTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAagg                                                                                                                                                   >  1:8684579/1‑151 (MQ=255)
     tGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAaggg                                                                                                                                                  <  2:2733683/151‑1 (MQ=255)
      gggCAAGGAGCTGCGCTCCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAagggg                                                                                                                                                 <  2:9669513/151‑1 (MQ=255)
      gggCAAGGAGCGGCGCACCGGCGCGGCCACCGAGAACGGCCGCGCAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCTagggg                                                                                                                                                 <  2:1553117/151‑1 (MQ=255)
       ggCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAaggggg                                                                                                                                                <  2:11003727/151‑2 (MQ=255)
        gCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAagggggt                                                                                                                                               <  2:4070347/151‑3 (MQ=255)
         cAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTc                                                                                                                                               >  1:1010296/1‑151 (MQ=255)
          aaGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGGTCATGCTGATCGGCGAGACCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCa                                                                                                                                              >  2:596849/1‑151 (MQ=255)
           aGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCaa                                                                                                                                             >  1:431972/1‑150 (MQ=255)
             gAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCaagg                                                                                                                                           <  1:5193386/151‑4 (MQ=255)
             gAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCaagg                                                                                                                                           >  2:2475049/1‑148 (MQ=255)
              aGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCaaggc                                                                                                                                          <  2:10538374/151‑5 (MQ=255)
              aGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCaaggc                                                                                                                                          <  2:9193955/151‑5 (MQ=255)
                cTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTccaggcta                                                                                                                                        >  2:949289/1‑144 (MQ=255)
                     cACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGACGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCAcc                                                                                                                                   >  1:4612724/1‑151 (MQ=255)
                     cACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCAcc                                                                                                                                   >  1:1656658/1‑151 (MQ=255)
                     cACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCAcc                                                                                                                                   <  2:4333682/151‑1 (MQ=255)
                        cGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTc                                                                                                                                >  2:1909003/1‑151 (MQ=255)
                         ggCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCt                                                                                                                               <  1:8092683/151‑1 (MQ=255)
                           cgcgGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTAc                                                                                                                             <  1:2165611/151‑1 (MQ=255)
                              ggCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGAc                                                                                                                          <  1:6421998/151‑1 (MQ=255)
                               gCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGAcc                                                                                                                         <  2:2918416/151‑1 (MQ=255)
                                 cACCGAGAACGGCCGCGAATTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCAGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGAACACCGTCTACGACCGc                                                                                                                       >  2:1614387/1‑151 (MQ=255)
                                   ccGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGcac                                                                                                                     >  2:5081577/1‑151 (MQ=255)
                                    cGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGcacc                                                                                                                    >  2:10693808/1‑151 (MQ=255)
                                       gaACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCaccacc                                                                                                                 <  2:8797238/151‑1 (MQ=255)
                                         aCGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCt                                                                                                               <  1:3333694/151‑1 (MQ=255)
                                          cGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTg                                                                                                              >  1:903744/1‑151 (MQ=255)
                                            gCCGCGACGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCACCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggt                                                                                                            <  1:2833976/151‑3 (MQ=255)
                                            gCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggt                                                                                                            <  2:2286639/151‑3 (MQ=255)
                                             ccGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCGGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtg                                                                                                           >  2:6240328/1‑148 (MQ=255)
                                              cgcgAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtgg                                                                                                          >  2:11367490/1‑147 (MQ=255)
                                                 gAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtggaca                                                                                                       >  2:6458745/1‑144 (MQ=255)
                                                                       tCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaa                                                                                >  1:1254142/1‑151 (MQ=255)
                                                                        cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaag                                                                               >  2:5636114/1‑151 (MQ=255)
                                                                          tGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaagaa                                                                             <  1:7722774/151‑1 (MQ=255)
                                                                           gCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAc                                                                            >  2:8318948/1‑151 (MQ=255)
                                                                           gCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAc                                                                            <  1:10326630/151‑1 (MQ=255)
                                                                             tGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcct                                                                          <  2:3571213/151‑2 (MQ=255)
                                                                             tGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcct                                                                          <  1:8905687/151‑2 (MQ=255)
                                                                              gATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctg                                                                         <  1:6646749/151‑3 (MQ=255)
                                                                              gATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctg                                                                         >  2:4183747/1‑149 (MQ=255)
                                                                              gATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctg                                                                         <  1:4204351/151‑3 (MQ=255)
                                                                               aTCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctgg                                                                        <  2:11246936/151‑4 (MQ=255)
                                                                                  ggCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCtggtgg                                                                    <  2:7617991/151‑1 (MQ=255)
                                                                                     gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAg                                                                 >  2:4454582/1‑151 (MQ=255)
                                                                                     gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCGGGTGGAAg                                                                 >  2:2209169/1‑151 (MQ=255)
                                                                                      agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCCCGGTGAAGAAGAACCTGGTGGAAgg                                                                >  2:3290497/1‑151 (MQ=255)
                                                                                      agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAgg                                                                <  1:3027922/151‑1 (MQ=255)
                                                                                      agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAgg                                                                <  1:4092827/151‑1 (MQ=255)
                                                                                       gaACAGCCGGGACGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggc                                                               <  2:9015614/151‑1 (MQ=255)
                                                                                          cAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcg                                                            >  1:397075/1‑151 (MQ=255)
                                                                                           aGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcgg                                                           >  1:6243114/1‑151 (MQ=255)
                                                                                            gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggc                                                          <  1:4558199/151‑1 (MQ=255)
                                                                                             ccGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGCCAAGGCCGTGGCCCCGGTGACGAAGAACCTGGTGGAAggcgcggcg                                                         >  1:10508264/1‑151 (MQ=255)
                                                                                              cGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGGGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGACggcgcggcgc                                                        <  1:2911860/151‑1 (MQ=255)
                                                                                               gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggcgcg                                                       >  1:1058109/1‑150 (MQ=255)
                                                                                               gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCt                                                       <  2:2359068/151‑1 (MQ=255)
                                                                                                    gTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGa                                                  <  2:5566876/151‑1 (MQ=255)
                                                                                                     tGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTATACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGAt                                                 <  1:606373/151‑1 (MQ=255)
                                                                                                      ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATc                                                <  1:5770943/151‑1 (MQ=255)
                                                                                                         ccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcgcg                                             <  1:4886662/151‑1 (MQ=255)
                                                                                                           aggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGt                                           <  1:8958556/149‑1 (MQ=255)
                                                                                                            ggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGACGAAGACCCTGGTGGAAGGCGCGGCGCTGGGGATCGCGGTg                                          >  1:4384650/2‑151 (MQ=255)
                                                                                                            ggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTg                                          <  1:1452031/150‑1 (MQ=255)
                                                                                                            ggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTg                                          <  1:1733868/150‑1 (MQ=255)
                                                                                                              cgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGct                                        >  2:4704572/1‑151 (MQ=255)
                                                                                                              cgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGct                                        <  1:3286997/151‑1 (MQ=255)
                                                                                                              cgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGct                                        <  1:2403033/151‑1 (MQ=255)
                                                                                                              cgGTCAGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAGGGCGCGGCGCTGGTGATCGCGGTGct                                        <  1:5907935/151‑1 (MQ=255)
                                                                                                                   ggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                   >  2:5798692/2‑151 (MQ=255)
                                                                                                                   ggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                   <  1:357241/150‑1 (MQ=255)
                                                                                                                     ccAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTg                                 >  2:10085863/1‑151 (MQ=255)
                                                                                                                      cAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGt                                >  2:3045885/1‑151 (MQ=255)
                                                                                                                      cAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGt                                >  1:8463091/1‑151 (MQ=255)
                                                                                                                        gcgcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTc                              <  2:1456007/151‑1 (MQ=255)
                                                                                                                         cgcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTcc                             <  2:7243054/151‑1 (MQ=255)
                                                                                                                          gcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCt                            <  2:4215564/151‑1 (MQ=255)
                                                                                                                          gcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCt                            >  1:54927/1‑151 (MQ=255)
                                                                                                                          gcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCt                            <  1:6899605/151‑1 (MQ=255)
                                                                                                                            cTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                          <  1:10934744/151‑1 (MQ=255)
                                                                                                                            cTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                          <  1:7085576/151‑1 (MQ=255)
                                                                                                                               gaggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCa                       <  1:9792844/151‑1 (MQ=255)
                                                                                                                               gaggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCa                       <  2:1254422/151‑1 (MQ=255)
                                                                                                                                 ggagATCTACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAAc                     <  2:494295/151‑1 (MQ=255)
                                                                                                                                 ggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAAc                     <  2:3755677/151‑1 (MQ=255)
                                                                                                                                  gagaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACa                    <  2:8018091/151‑1 (MQ=255)
                                                                                                                                  gagaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACa                    <  1:10643040/151‑1 (MQ=255)
                                                                                                                                        aaCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcg              >  1:9334923/1‑151 (MQ=255)
                                                                                                                                        aaCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcg              >  1:7286742/1‑151 (MQ=255)
                                                                                                                                         aCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACCAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgc             >  2:6476041/1‑151 (MQ=255)
                                                                                                                                         aCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgc             <  2:7939856/151‑1 (MQ=255)
                                                                                                                                          ccGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgcg            <  1:8045096/151‑1 (MQ=255)
                                                                                                                                          ccGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgcg            <  1:11096835/151‑1 (MQ=255)
                                                                                                                                           cGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcgg           >  1:82126/1‑151 (MQ=255)
                                                                                                                                            ggACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggc          >  2:5311752/1‑151 (MQ=255)
                                                                                                                                            ggACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggc          >  1:1985780/1‑151 (MQ=255)
                                                                                                                                                ccTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTg      <  2:8593605/151‑1 (MQ=255)
                                                                                                                                                   ccccAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCa  >  2:9716185/2‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCA  >  NC_008463/2787669‑2787967

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: