Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,463,550 G→C 7.7% R270G (CGT→GGT)  PP_3078 ← ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,463,5500GC 7.7% 200.1 / 5.9 65R270G (CGT→GGT) PP_3078ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base C (3/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

GCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTACCGG  >  NC_002947/3463404‑3463697
                                                                                                                                                   |                                                                                                                                                   
gCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGt                                                                                                                                                    <  5:216244/148‑1 (MQ=255)
 cccAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGc                                                                                                                                                  >  4:384907/1‑149 (MQ=255)
  ccAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCATGGTGGGCACGTGACCTGc                                                                                                                                                  >  1:322095/1‑149 (MQ=255)
   cAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGcc                                                                                                                                                 >  7:9525/1‑148 (MQ=255)
    aaGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCcca                                                                                                                                               >  6:75198/1‑149 (MQ=255)
    aaGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCcca                                                                                                                                               <  3:384907/149‑1 (MQ=255)
    aaGGTGTAGGCCGTGGCGGAGAACTCGCGGGTGTAGTCGTAGCGACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCcca                                                                                                                                               <  4:262506/149‑1 (MQ=255)
     aGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCccac                                                                                                                                              <  5:209328/149‑1 (MQ=255)
      ggTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCccacc                                                                                                                                             <  1:266294/149‑1 (MQ=255)
         gtAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATg                                                                                                                                          >  7:80738/1‑149 (MQ=255)
          tAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGCATGGTGGGCACGTGACCTGCccgccag                                                                                                                                           <  2:322095/148‑2 (MQ=255)
          tAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGc                                                                                                                                         <  8:9525/149‑1 (MQ=255)
           aGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCTCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGt                                                                                                                                                    >  5:397063/1‑137 (MQ=255)
           aGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGt                                                                                                                                                    <  6:397063/137‑1 (MQ=255)
            ggCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGcc                                                                                                                                        <  5:75198/148‑1 (MQ=255)
            ggCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGcc                                                                                                                                        <  6:286263/148‑1 (MQ=255)
            ggCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCg                                                                                                                                       <  8:126987/149‑1 (MQ=255)
                gTGGCGGTGTACTCGCGGTTGTAGTCGTATGCACCTCCTTTGAGTACCTGGCGAGCCATTTATGGTTCTCAGAAGTAATCGATGCGCTGTTTGTCTAATTTGTAGCGGCCACGTC‑TTGCGTGCAGTTCACTTGCCCAACATTCAGGGTc                                                                                                                                   <  4:390562/149‑1 (MQ=255)
                gTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTc                                                                                                                                   <  7:167346/149‑1 (MQ=255)
                    cGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGt                                                                                                                               <  6:341873/149‑1 (MQ=255)
                              gcgGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTg                                                                                                                      <  5:108896/148‑1 (MQ=255)
                                        tCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCggg                                                                                                                                     >  6:209215/1‑123 (MQ=255)
                                        tCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCggg                                                                                                                                     <  5:209215/123‑1 (MQ=255)
                                         cGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGcc                                                                                                           <  6:219551/148‑1 (MQ=255)
                                              ccACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTc                                                                                                                                   <  7:97108/119‑1 (MQ=255)
                                              ccACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTc                                                                                                                                   >  8:97108/1‑119 (MQ=255)
                                                               cTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGt                                                                                    <  7:146857/149‑1 (MQ=255)
                                                                  gcgcggCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGGGGGCAGGTCACGTGCCCAGCATGCCGGGTCGCGTTCGTAGGTGATGCTGCGGCCGTTGTCGATGGGCCCAATGCGGTATg                                                                                 <  5:195442/149‑1 (MQ=255)
                                                                  gcgcgcCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGCTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATg                                                                                 <  5:189553/149‑1 (MQ=255)
                                                                   tcgagtaCTTCGGCGACACCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATgg                                                                                >  3:8928/7‑149 (MQ=255)
                                                                    gcgcCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATg                                                                                                 <  7:9917/131‑1 (MQ=255)
                                                                    gcgcCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATg                                                                                                 >  8:9917/1‑131 (MQ=255)
                                                                       ccACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCg                                                                            >  7:35682/1‑149 (MQ=255)
                                                                        cACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGc                                                                           >  8:374088/1‑149 (MQ=255)
                                                                          cTTCGGCGACACCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCGCCCGGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTg                                                                         >  1:131681/1‑149 (MQ=255)
                                                                           ttCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGc                                                                        >  4:344290/1‑149 (MQ=255)
                                                                                tgtCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATgg                                                                                >  8:135398/1‑136 (MQ=255)
                                                                                tgtCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATgg                                                                                <  7:135398/136‑1 (MQ=255)
                                                                                tgtCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACt                                                                   <  7:374088/149‑1 (MQ=255)
                                                                                  tCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTgg                                                                 >  5:212251/1‑149 (MQ=255)
                                                                                     ccGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACt                                                                   <  1:150436/144‑1 (MQ=255)
                                                                                     ccGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACt                                                                   >  2:150436/1‑144 (MQ=255)
                                                                                       gAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGt                                                            >  7:83416/1‑149 (MQ=255)
                                                                                             cTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGc                                                      >  7:149030/1‑149 (MQ=255)
                                                                                                gATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCgg                                                   >  6:16847/1‑149 (MQ=255)
                                                                                                gATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCgg                                                   >  6:1634/1‑149 (MQ=255)
                                                                                                     gcAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCg                                              >  7:161695/1‑149 (MQ=255)
                                                                                                      cAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGt                                             >  2:315491/1‑149 (MQ=255)
                                                                                                      cAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGt                                             >  6:183851/1‑149 (MQ=255)
                                                                                                        gTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTg                                           >  2:127923/1‑149 (MQ=255)
                                                                                                               gAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGt                                    >  6:11315/1‑149 (MQ=255)
                                                                                                               gAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGt                                    >  2:261666/1‑149 (MQ=255)
                                                                                                                        gCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCgg                                                   >  7:239714/1‑125 (MQ=255)
                                                                                                                        gCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCgg                                                   <  8:239714/125‑1 (MQ=255)
                                                                                                                               cGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTc                    <  6:212251/149‑1 (MQ=255)
                                                                                                                                   atggtggGCACGTGACCTGCCCGCCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTg                 >  4:305314/2‑149 (MQ=255)
                                                                                                                                   atggtggGCACGTGACCTGCCCGCCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTg                 >  2:56868/2‑149 (MQ=255)
                                                                                                                                       cGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCg             <  5:16847/149‑1 (MQ=255)
                                                                                                                                         ggCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCGGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCggg           <  8:80738/149‑1 (MQ=255)
                                                                                                                                          gCACGTGACCTGCCCGCCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGc          <  1:56868/149‑1 (MQ=255)
                                                                                                                                           cAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCTGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCa         >  6:367257/1‑149 (MQ=255)
                                                                                                                                              gTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGa                                                                     >  2:108169/1‑86 (MQ=255)
                                                                                                                                              gTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGa                                                                     <  1:108169/86‑1 (MQ=255)
                                                                                                                                                 aCGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTACCg   >  5:255077/1‑149 (MQ=255)
                                                                                                                                                  cGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTACCgg  <  4:8928/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCCGTTGTCGATGCGCCCAATGCGGTATGGCCCGCTGCCGACTGGCTTGTCGAAGCCGGCACCGTTGGCAAAGTCGCGCTGCTGCAGGTCGTGCTCGGGCATTACCGG  >  NC_002947/3463404‑3463697

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑