Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,463,450 C→A 6.0% G303V (GGC→GTC)  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,4500CA 6.0% 214.2 / 3.5 67G303V (GGC→GTC) PP_3078ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (26/36);  new base A (2/2);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

GCTGGTCGAGGTTGAACACGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCC  >  NC_002947/3463304‑3463592
                                                                                                                                                  |                                                                                                                                               
gCTGGTCGAGGTTGAACACGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCa                                                                                                                                               <  5:24521/149‑1 (MQ=255)
    gTCGAGGTTGAACACGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGc                                                                                                                                           >  6:216244/1‑149 (MQ=255)
    gTCGAGGTTGAACACGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGc                                                                                                                                           <  3:319728/149‑1 (MQ=255)
      cGAGGTTGAACACGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCt                                                                                                                                         >  4:129738/1‑149 (MQ=255)
        aGGTTGAACACGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCttt                                                                                                                                       <  5:257941/149‑1 (MQ=255)
        aGGTTGAACACGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCttt                                                                                                                                       <  6:84903/149‑1 (MQ=255)
        aGGTTGAACACGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCttt                                                                                                                                       <  4:153375/149‑1 (MQ=255)
          gTTGAACACGAACCCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTTTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGa                                                                                                                                     <  8:270817/149‑1 (MQ=255)
           ttGAACACGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAg                                                                                                                                    <  1:108455/149‑1 (MQ=255)
           ttGAACACGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAg                                                                                                                                    >  6:362759/1‑149 (MQ=255)
               acacGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTAcc                                                                                                                                >  6:104040/1‑149 (MQ=255)
               acacGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTAcc                                                                                                                                >  6:108896/1‑149 (MQ=255)
               acacGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTAc                                                                                                                                 <  6:69328/148‑1 (MQ=255)
                        ccTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGgcgcgc                                                                                                                        >  4:302346/1‑148 (MQ=255)
                             gCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTc                                                                                                                  >  5:188280/1‑149 (MQ=255)
                              ccACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCg                                                                                                                 <  3:129738/149‑1 (MQ=255)
                                aCCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGgt                                                                                                               <  6:188280/149‑1 (MQ=255)
                                aCCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGgt                                                                                                               <  1:93991/149‑1 (MQ=255)
                                 ccTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGgtg                                                                                                              >  3:390562/1‑149 (MQ=255)
                                      ggTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGcc                                                                                                         <  7:176350/149‑1 (MQ=255)
                                           ggCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCTCCGAAg                                                                                                     <  6:307410/148‑1 (MQ=255)
                                           ggCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGGACCTCGCGCGCCACTTCGGTGTCGCCGAAg                                                                                                     >  2:310602/1‑148 (MQ=255)
                                             ccGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTAc                                                                                                  >  7:27436/1‑149 (MQ=255)
                                                ggCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTc                                                                                                >  6:209328/1‑148 (MQ=255)
                                                 gCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGa                                                                                              >  7:126987/1‑149 (MQ=255)
                                                                   cAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGaaa                                                                             <  5:104040/148‑1 (MQ=255)
                                                                      gcggcCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGGTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTg                                                                          >  5:286263/1‑148 (MQ=255)
                                                                         gcCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGc                                                                      >  6:189553/1‑149 (MQ=255)
                                                                         gcCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGAAGGCCGTGGCGGAGAACTCGCCGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCCCTTCGGTGTCGCCCAACTACTCGATGCCCACTTTCTCCAAATACTAGc                                                                      >  6:195442/1‑149 (MQ=255)
                                                                               gtcCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCAc                                                                >  5:272674/1‑149 (MQ=255)
                                                                                tcCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACg                                                               <  6:151368/149‑1 (MQ=255)
                                                                                  cAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGTCACgg                                                              <  8:27436/148‑1 (MQ=255)
                                                                                        cgcACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCg                                                       <  5:362759/148‑1 (MQ=255)
                                                                                           aCCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGATTGTCGAAATTGTAGCGGCCACGGC‑TGGCggg                                                     >  7:285817/1‑147 (MQ=255)
                                                                                             cGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGGGGGCAgg                                                 <  6:272674/149‑1 (MQ=255)
                                                                                               ttGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGt                                                >  5:219551/1‑148 (MQ=255)
                                                                                                    gCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGt                                           <  5:216244/148‑1 (MQ=255)
                                                                                                     cccAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGc                                         >  4:384907/1‑149 (MQ=255)
                                                                                                      ccAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGgc                                                                                                                            <  4:53480/66‑1 (MQ=255)
                                                                                                      ccAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGgc                                                                                                                            >  3:53480/1‑66 (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)
                                                                                                                  ccGTGGCGGAGAACTCGCGGTTGTAGTCGTAGACACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTg                                                                          >  7:37399/1‑104 (MQ=255)
                                                                                                                  ccGTGGCGGAGAACTCGCGGTTGTAGTCGTAGACACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTg                                                                          <  4:195858/104‑1 (MQ=255)
                                                                                                                  ccGTGGCGGAGAACTCGCGGTTGTAGTCGTAGACACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTg                                                                          <  8:37399/104‑1 (MQ=255)
                                                                                                                  ccGTGGCGGAGAACTCGCGGTTGTAGTCGTAGACACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTg                                                                          >  3:195858/1‑104 (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)
                                                                                                                                   cgGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGc                                                    <  4:133411/108‑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)
                                                                                                                                                  |                                                                                                                                               
GCTGGTCGAGGTTGAACACGAAACCCTGGGCCACCTGCGGTTTGGCCGGGCCCAGGTGGGCGCGCTGCAGGCGGCCGTCGTCCAGTTGCGCACCGTTGTAGCCCAAGGTGTAGGCCGTGGCGGAGAACTCGCGGTTGTAGTCGTAGCCACCGCCTTTGAGTACCTGGCGCGCCACTTCGGTGTCGCCGAAGTACTCGATGCGCAGTTTGTCGAAATTGTAGCGGCCACGGC‑TGGCGGGCAGGTCACGTGCCCACCATGCCGGGTCGCGTTCGAAGGTGATGCTGCGGCC  >  NC_002947/3463304‑3463592

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: ‑