Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,670,537 Δ1 bp coding (537/1020 nt) PA14_RS16675 → ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,670,5340T.99.7% 1854.9 / ‑2.3 377coding (534/1020 nt)PA14_RS16675ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (184/192);  total (184/193)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 392 total aligned reads displayed.

GCCCGCGACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTTCCTGGAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTT  >  NC_008463/3670385‑3670677
                                                                                                                                                     |                                                                                                                                                 
gCCCGCGACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGt                                                                                                                                                   <  1:8297470/150‑1 (MQ=255)
gCCCGCGACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGt                                                                                                                                                   >  2:4951823/1‑150 (MQ=255)
  ccGCGACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGttt                                                                                                                                                 <  1:1801975/150‑1 (MQ=255)
  ccGCGACGTGCTCGCCCGGGTGATGTACGGCGTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGttt                                                                                                                                                 <  2:7647957/150‑1 (MQ=255)
   cgcgACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGtttc                                                                                                                                                >  2:6661524/1‑149 (MQ=255)
   cgcgACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGtttc                                                                                                                                                >  1:7994681/1‑149 (MQ=255)
       aCGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTgg                                                                                                                                           <  1:61426/150‑1 (MQ=255)
              cGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTg                                                                                                                                   <  1:9704120/150‑1 (MQ=255)
              cGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTg                                                                                                                                   >  2:2342063/1‑150 (MQ=255)
               gCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTgg                                                                                                                                  <  1:7449338/150‑1 (MQ=255)
                cccGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGt                                                                                                                                 >  1:8765133/1‑150 (MQ=255)
                 ccGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGAGCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTc                                                                                                                                >  1:10110162/1‑150 (MQ=255)
                 ccGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTc                                                                                                                                >  1:4530408/1‑150 (MQ=255)
                 ccGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTc                                                                                                                                <  2:9247626/150‑1 (MQ=255)
                   gggTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCgg                                                                                                                              >  2:2850417/1‑150 (MQ=255)
                     gTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCgggg                                                                                                                            >  2:7871274/1‑150 (MQ=255)
                       gATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCt                                                                                                                          <  2:1777635/150‑1 (MQ=255)
                       gATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCt                                                                                                                          <  1:963477/150‑1 (MQ=255)
                        aTTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTg                                                                                                                         <  2:3806804/150‑1 (MQ=255)
                          ttACGGCTTCCGCGTGTCGATCCTCGTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGcc                                                                                                                       <  2:4373299/150‑1 (MQ=255)
                           tACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCg                                                                                                                      >  1:4391014/1‑150 (MQ=255)
                            aCGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCGGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCACGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCgg                                                                                                                     <  1:2322356/150‑1 (MQ=255)
                            aCGGCGTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCgg                                                                                                                     <  2:5158110/150‑1 (MQ=255)
                                 gtCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCTGCGCCCTGCAAGGCTAATACGGCTGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGATGTCTGGTCGGGGCTGCCGGTGctc                                                                                                                <  2:5385582/149‑1 (MQ=38)
                                 gtCCGCGTGGCGATCCTCTACGCCCTGACCCCGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTCCGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGctc                                                                                                                <  2:9926107/149‑1 (MQ=255)
                                  tCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGctct                                                                                                               <  1:3340810/150‑1 (MQ=255)
                                    cgcgTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTAc                                                                                                             >  2:7155685/1‑150 (MQ=255)
                                     gcgTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTAcc                                                                                                            <  2:9472328/150‑1 (MQ=255)
                                       gtgtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCCAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctg                                                                                                          >  2:5266653/1‑150 (MQ=255)
                                       gtgtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctg                                                                                                          >  2:566917/1‑150 (MQ=255)
                                        tgtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgc                                                                                                         >  1:3953049/1‑150 (MQ=255)
                                        tgtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgc                                                                                                         <  2:9743517/150‑1 (MQ=255)
                                         gtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgct                                                                                                        <  1:2757925/150‑1 (MQ=255)
                                         gtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCGGGTCGGGGCTGCCGGTGCTCTACctgct                                                                                                        >  2:5959212/1‑150 (MQ=255)
                                          tCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgctg                                                                                                       <  1:3427530/150‑1 (MQ=255)
                                          tCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgctg                                                                                                       >  1:3080622/1‑150 (MQ=255)
                                           cGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGa                                                                                                      >  2:8674291/1‑150 (MQ=255)
                                                 tctTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATcc                                                                                                <  1:9928342/150‑1 (MQ=255)
                                                  ctTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCt                                                                                               >  1:5244596/1‑150 (MQ=255)
                                                   ttCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTg                                                                                              <  1:1351163/150‑1 (MQ=255)
                                                    tCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTgg                                                                                             <  1:1301961/150‑1 (MQ=255)
                                                     cGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGc                                                                                            <  1:4244392/150‑1 (MQ=255)
                                                     cGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGc                                                                                            <  1:6735965/150‑1 (MQ=255)
                                                       ccctgaccctgaCCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCa                                                                                          <  1:1744294/150‑1 (MQ=255)
                                                       ccctgaccctgaCCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCa                                                                                          <  1:7143024/150‑1 (MQ=255)
                                                          tgaccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCt                                                                                       <  1:2888545/150‑1 (MQ=255)
                                                           gaccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCtt                                                                                      <  2:10239268/150‑1 (MQ=255)
                                                             ccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCg                                                                                    <  1:2125076/150‑1 (MQ=255)
                                                             ccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCg                                                                                    >  1:2815952/1‑150 (MQ=255)
                                                             ccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCg                                                                                    >  1:5468315/1‑150 (MQ=255)
                                                               ctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTc                                                                                  >  2:8937235/1‑150 (MQ=255)
                                                               ctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTc                                                                                  <  1:10031120/150‑1 (MQ=255)
                                                                tgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTcc                                                                                 <  2:5730575/150‑1 (MQ=255)
                                                                 gaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCa                                                                                >  2:4997313/1‑150 (MQ=255)
                                                                  accGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAg                                                                               >  1:2581044/1‑150 (MQ=255)
                                                                  accGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAg                                                                               >  2:1203428/1‑150 (MQ=255)
                                                                   ccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGc                                                                              >  2:9302418/1‑150 (MQ=255)
                                                                    cGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGcc                                                                             >  2:643805/1‑150 (MQ=255)
                                                                     gTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCg                                                                            >  2:2479441/1‑150 (MQ=255)
                                                                      tGCTCAGGTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGa                                                                           <  2:2722101/150‑1 (MQ=255)
                                                                       gCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGaa                                                                          <  2:4172400/150‑1 (MQ=255)
                                                                        cTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAAc                                                                         <  1:8279020/150‑1 (MQ=255)
                                                                          cAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACtt                                                                       <  1:1035566/150‑1 (MQ=255)
                                                                           aGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTc                                                                      >  2:3593059/1‑150 (MQ=255)
                                                                            gTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCt                                                                     <  2:6402332/150‑1 (MQ=255)
                                                                             ttCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCtg                                                                    <  1:6944354/150‑1 (MQ=255)
                                                                                 ctGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCtggtgg                                                                <  1:6504496/150‑1 (MQ=255)
                                                                                    aTCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgctg                                                             >  1:5115154/1‑150 (MQ=255)
                                                                                       ggCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGctgctg                                                          <  1:8194247/150‑1 (MQ=255)
                                                                                         cGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTggg                                                        <  1:7722776/150‑1 (MQ=255)
                                                                                           tCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTAGCTGCTGGGca                                                      >  1:4844609/1‑150 (MQ=255)
                                                                                              ccGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcatca                                                   <  1:84721/150‑1 (MQ=255)
                                                                                              ccGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcatca                                                   <  1:9041303/150‑1 (MQ=255)
                                                                                               cGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcatcat                                                  <  1:3539829/150‑1 (MQ=255)
                                                                                               cGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcatcat                                                  <  2:8740888/150‑1 (MQ=255)
                                                                                                   ggCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgctg                                              >  2:7300257/1‑150 (MQ=255)
                                                                                                      gcCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGctgctg                                           >  2:4718299/1‑150 (MQ=255)
                                                                                                       cccTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGt                                          >  1:4700523/1‑150 (MQ=255)
                                                                                                         cTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGttt                                        >  2:3324630/1‑150 (MQ=255)
                                                                                                         cTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGttt                                        <  2:3021733/150‑1 (MQ=255)
                                                                                                         cTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGttt                                        >  1:3200403/1‑150 (MQ=255)
                                                                                                            cAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTcc                                     <  1:9323178/150‑1 (MQ=255)
                                                                                                             aaGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCt                                    >  1:1530344/1‑150 (MQ=255)
                                                                                                               ggCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTgg                                  <  2:140805/150‑1 (MQ=255)
                                                                                                                 ctactaCGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGAt                                <  2:1688777/150‑1 (MQ=255)
                                                                                                                   actccGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGATGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGTTCATCCTGGCCAGCTTCGTCCTGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGa                              <  2:3514834/150‑1 (MQ=255)
                                                                                                                         ccggcTGGCTCGACCTGCTCGGCCAGCGTTTCCAGGAAGTCAGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCtgg                        <  2:9124828/149‑1 (MQ=255)
                                                                                                                            gcTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCtggtgg                     >  1:3834226/1‑150 (MQ=255)
                                                                                                                             cTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGa                    <  2:8399185/150‑1 (MQ=255)
                                                                                                                             cTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCGTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGa                    <  1:3628748/150‑1 (MQ=255)
                                                                                                                              tGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGAc                   >  2:6732097/1‑150 (MQ=255)
                                                                                                                                 gTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTg                <  1:5624486/150‑1 (MQ=255)
                                                                                                                                    gACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTc             >  1:5884103/1‑150 (MQ=255)
                                                                                                                                     aCCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTcc            >  1:5769161/1‑150 (MQ=255)
                                                                                                                                     aCCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTcc            >  2:4783823/1‑150 (MQ=255)
                                                                                                                                     aCCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTcc            >  1:7622778/1‑150 (MQ=255)
                                                                                                                                     aCCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTcc            >  2:5868941/1‑150 (MQ=255)
                                                                                                                                       cTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCcgc          <  1:10562074/150‑1 (MQ=255)
                                                                                                                                            cGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGACCTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGa     <  2:8021242/150‑1 (MQ=255)
                                                                                                                                              gCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGt   >  1:7351732/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCCCGCGACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTTCCTGGAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTT  >  NC_008463/3670385‑3670677

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: