Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,545,022 C→T 5.4% G138D (GGC→GAC)  PP_3130 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,545,0220CT 5.4% 229.9 / 2.2 74G138D (GGC→GAC) PP_3130hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/36);  new base T (2/2);  total (36/38)
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.94e-01

GAATGCACGCAGCTCGGACGGGTTGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACGGATGTCAGCCTGCAGGTCGAAGTC  >  NC_002947/3544877‑3545164
                                                                                                                                                 |                                                                                                                                              
gAATGCACGCAGCTCGGACGGGTTGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATgccgc                                                                                                                                             >  5:336882/1‑149 (MQ=255)
  aTGCACGCAGCTCGGACGGGTTGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGccgccg                                                                                                                                           >  5:209809/1‑149 (MQ=255)
  aTGCACGCAGCTCGGACGGGTTGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGccgccg                                                                                                                                           >  2:59614/1‑149 (MQ=255)
      aCGCAGCTCGGACGGGTTGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGa                                                                                                                                       <  4:23933/149‑1 (MQ=255)
           gCTCGGACGGGTTGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTgc                                                                                                                                  <  3:119036/149‑1 (MQ=255)
               ggACGGGTTGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCg                                                                                                                              >  5:177853/1‑149 (MQ=255)
                gACGGGTTGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCgg                                                                                                                             >  7:78002/1‑149 (MQ=255)
                 aCGGGTTGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCggg                                                                                                                            <  8:78002/149‑1 (MQ=255)
                  cGGGTTGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGt                                                                                                                           <  8:94688/149‑1 (MQ=255)
                     gTTGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGa                                                                                                                        <  6:177853/149‑1 (MQ=255)
                      ttGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGCCGGGCAACAATTGCGGGGCGCCCATGTCGCCGGCGATTTGCGCCGGGTCGaa                                                                                                                       >  6:218811/1‑149 (MQ=255)
                      ttGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGACGGGCAACAATTGCGGGGCGCCCATGTCGCCGGCGATTTGCGCCGGGTCGaa                                                                                                                       >  6:213557/1‑149 (MQ=255)
                        gCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAAc                                                                                                                      <  8:256355/148‑1 (MQ=255)
                        gCCGCGCTCGAAATAGACCTTGGCCAGT‑GCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACt                                                                                                                     <  2:52488/148‑1 (MQ=255)
                          cgcgCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTgg                                                                                                                   <  8:170133/149‑1 (MQ=255)
                          cgcgCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTgg                                                                                                                   <  6:381381/149‑1 (MQ=255)
                            cgcTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGaa                                                                                                                 <  7:71978/149‑1 (MQ=255)
                             gcTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAac                                                                                                                <  6:336882/149‑1 (MQ=255)
                               tCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAacac                                                                                                              >  5:152222/1‑149 (MQ=255)
                                 gAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAacac                                                                                                              <  3:245079/147‑1 (MQ=255)
                                 gAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCa                                                                                                            <  6:152222/149‑1 (MQ=255)
                                    aTAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGgc                                                                                                         <  4:14176/149‑1 (MQ=255)
                                        aCCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGgccgcc                                                                                                     >  3:326952/1‑149 (MQ=255)
                                         ccTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCa                                                                                                    >  5:86930/1‑149 (MQ=255)
                                         ccTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCa                                                                                                    >  5:92394/1‑149 (MQ=255)
                                          cTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCAc                                                                                                   >  6:366010/1‑149 (MQ=255)
                                          cTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCAc                                                                                                   >  4:360080/1‑149 (MQ=255)
                                              gCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCAt                                                                                               <  5:366010/149‑1 (MQ=255)
                                                cAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGa                                                                                             >  7:112645/1‑149 (MQ=255)
                                                     gCGGGTTGAAACGCATTTCGATGCGTAGCCCCTGGTCAGCCTTGCCGCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGCTTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGacacag                                                                                        >  6:104041/1‑147 (MQ=255)
                                                     gCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACCCCAGGCCGCCACCCATGACACCg                                                                                        >  4:369705/1‑149 (MQ=255)
                                                     gCGGGTCGAAACGCATTTCGATGCGTAGCCCCTGGTCAGCCTTGCCCCCCGTCCAGCCGCGCTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGCCGATTTGCGCCGGGTCGAACTGCAACACCAGGCCGCCACCCATGACACCg                                                                                        >  7:322420/1‑149 (MQ=255)
                                                           cGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTgg                                                                                  >  6:42013/1‑149 (MQ=255)
                                                            gAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTggc                                                                                 >  3:216917/1‑149 (MQ=255)
                                                              aaCGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTggcgg                                                                               >  5:129679/1‑149 (MQ=255)
                                                                cGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCt                                                                             <  8:49076/149‑1 (MQ=255)
                                                                    tttCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTgg                                                                         <  5:224109/149‑1 (MQ=255)
                                                                     ttCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCAc                                                                                                   >  1:318473/1‑122 (MQ=255)
                                                                     ttCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCAc                                                                                                   <  2:318473/122‑1 (MQ=255)
                                                                        gATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGt                                                                     <  2:312029/149‑1 (MQ=255)
                                                                         aTGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTc                                                                    <  6:92394/149‑1 (MQ=255)
                                                                         aTGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTc                                                                    <  6:86930/149‑1 (MQ=255)
                                                                             gTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTc                                                                    >  4:255917/1‑145 (MQ=255)
                                                                                gCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATc                                                              <  6:209809/148‑1 (MQ=255)
                                                                                   ccTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAg                                                          <  1:59614/149‑1 (MQ=255)
                                                                                            ccTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGccc                                                      >  3:328281/1‑144 (MQ=255)
                                                                                            ccTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGccc                                                      <  4:328281/144‑1 (MQ=255)
                                                                                            ccTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCt                                                  >  2:223228/1‑148 (MQ=255)
                                                                                            ccTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCt                                                  >  6:240757/1‑148 (MQ=255)
                                                                                             cTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCa                                                <  1:223228/149‑1 (MQ=255)
                                                                                                  cccACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGtt                                           >  3:334616/1‑149 (MQ=255)
                                                                                                  cccACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGtt                                           >  7:54192/1‑149 (MQ=255)
                                                                                                   ccACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGttt                                          >  6:200987/1‑149 (MQ=255)
                                                                                                     aCGTCCAGCCTCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTgc                                        <  7:96441/149‑1 (MQ=255)
                                                                                                     aCGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTgc                                        <  8:112645/149‑1 (MQ=255)
                                                                                                          cAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCctc                                   <  5:42013/149‑1 (MQ=255)
                                                                                                             ccGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTc                                >  2:154914/1‑149 (MQ=255)
                                                                                                             ccGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTc                                >  3:79410/1‑149 (MQ=255)
                                                                                                                  gTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTAc                           <  6:129679/149‑1 (MQ=255)
                                                                                                                   ttGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACg                          <  2:34202/149‑1 (MQ=255)
                                                                                                                      acGGGCAACAATTGCGGGGCGCCCATGTCGCCGGCGATTTTCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACGGAt                       <  5:218811/148‑1 (MQ=255)
                                                                                                                      acGGGCAACAATTGCGGGGCGCCCATGTCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACGGAt                       <  5:213557/148‑1 (MQ=255)
                                                                                                                        gggCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACGGATGt                     >  3:69574/1‑149 (MQ=255)
                                                                                                                          gCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACGGATGt                     >  8:185489/1‑147 (MQ=255)
                                                                                                                          gCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACGGATGTCa                   <  4:334616/149‑1 (MQ=255)
                                                                                                                            aacaaTTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACGGATGTCAGc                 <  4:326952/149‑1 (MQ=255)
                                                                                                                              caaTTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACGGATGTCAGCCt               >  5:270323/1‑149 (MQ=255)
                                                                                                                               aaTTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCGCAGGGTTTGCGCCTCTTCCTTACGGATAGCAGCCTg              >  2:236514/1‑149 (MQ=255)
                                                                                                                                aTTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACGGATGTCAGCCTg              >  1:106782/1‑148 (MQ=255)
                                                                                                                                 ttGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACGGATGTCAGCCTg              >  3:111739/1‑147 (MQ=255)
                                                                                                                                   gCGGGGCGCCCATGCCGCCGGCGATTTTCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTgg                                                                         >  1:15859/1‑86 (MQ=255)
                                                                                                                                   gCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTgg                                                                         <  2:15859/86‑1 (MQ=255)
                                                                                                                                   gCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACGGATGt                     <  7:247387/138‑1 (MQ=255)
                                                                                                                                   gCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACGGATGt                     >  8:247387/1‑138 (MQ=255)
                                                                                                                                      gggCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTgcg                                       >  2:172340/1‑117 (MQ=255)
                                                                                                                                      gggCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTgcg                                       <  1:172340/117‑1 (MQ=255)
                                                                                                                                           cccATGCCGCCGGAGATATGCGCCTGGTCGAACAGGAACACCATGCCGCCTCCCATGACACCGCTCTGGCGGCTGTGTGCGTCGAGATCGAACGCCCATCTCAGGCTTTGCGCCTCTTCCTTACGGATGTCAGCCTGCAGGTCGAAGTc  <  5:104041/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GAATGCACGCAGCTCGGACGGGTTGCCGCGCTCGAAATAGACCTTGGCCAGTGGCGGGTCGAAACGCATTTCGATGCGTCGCCCCTGGTCAGCCTTGCCCCACGTCCAGCCGCGGTTGTCGGGCAACAATTGCGGGGCGCCCATGCCGCCGGCGATTTGCGCCGGGTCGAACTGGAACACCAGGCCGCCACCCATGACACCGCTCTGGCGGCTGTGGGCGTCGAGATCGAAGGCCCAGCTCAGGGTTTGCGCCTCTTCCTTACGGATGTCAGCCTGCAGGTCGAAGTC  >  NC_002947/3544877‑3545164

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