Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,508,226 C→G 6.4% intergenic (+80/‑170) cycA → / → urtA D‑alanine/beta‑alanine/D‑serine/glycine permease/urea ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,508,2260CG 6.4% 238.6 / 4.3 75intergenic (+80/‑170)cycA/urtAD‑alanine/beta‑alanine/D‑serine/glycine permease/urea ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (34/36);  new base G (2/3);  total (36/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 = 9.53e-01

GCGGTGACCTACCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGCCCAACACCA  >  NC_002947/5508084‑5508371
                                                                                                                                                |                                                                                                                                                 
gCGGTGACCTACCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCgg                                                                                                                                              <  2:270756/148‑1 (MQ=255)
gCGGTGACCTACCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGc                                                                                                                                             <  5:174126/149‑1 (MQ=255)
 cGGTGACCTACCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCTCACCGCGGcc                                                                                                                                            <  2:281187/149‑1 (MQ=255)
  ggTGACCTACCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGGGCCCTTGTGTCGCGATGGGCCGCCGCGCGGccc                                                                                                                                           >  6:327069/1‑149 (MQ=255)
      aCCTACCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCCCGAAGGCACGGTGTACGCCCTCCGATGTGGGAGCGCCCTTGGGTCGCGATGGGCCGCGGTGCGGcccatc                                                                                                                                        <  7:152299/148‑4 (MQ=255)
      aCCTACCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCGGTGCGGcccatc                                                                                                                                        <  2:74591/148‑4 (MQ=255)
       ccTACCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATc                                                                                                                                      <  3:192437/149‑1 (MQ=255)
        cTACCAGCTCGTACGCAGCCGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCCATcc                                                                                                                                     >  4:287972/1‑149 (MQ=255)
        cTACCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGCATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCGGTGCGGcccatcgc                                                                                                                                      <  1:92096/149‑6 (MQ=255)
        cTACCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATcc                                                                                                                                     >  3:99532/1‑149 (MQ=255)
         tACCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTATGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCt                                                                                                                                    >  7:271185/1‑149 (MQ=255)
           ccAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCtgt                                                                                                                                  >  1:113449/1‑149 (MQ=255)
            cAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCtgtg                                                                                                                                 <  3:25650/149‑1 (MQ=255)
            cAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCtgtg                                                                                                                                 <  6:328936/149‑1 (MQ=255)
            cAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCtgtg                                                                                                                                 <  2:260682/149‑1 (MQ=255)
            cAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCtgtg                                                                                                                                 <  8:160449/149‑1 (MQ=255)
                  gTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGcacac                                                                                                                            >  3:326870/1‑148 (MQ=255)
                   tACGCAGCAGGCGCCATCCCCGCACGGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGcacac                                                                                                                            >  1:204290/1‑148 (MQ=255)
                       cagcagGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATc                                                                                                                       >  7:213369/1‑148 (MQ=255)
                         gcagGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTACCACGAAGGCACGGTGTACGCCGTCCGAAGTGGGAGCGCCCTTGTGTCGCGATGGGCCGAACCGCGGACCCAATCATGTGCACACAAATCCa                                                                                                                     >  7:312240/1‑148 (MQ=255)
                           agGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCa                                                                                                                  <  2:113449/149‑1 (MQ=255)
                             gcgcCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAgg                                                                                                                >  2:127010/1‑149 (MQ=255)
                              cgcCATCCCCGCACGGCCGTGCGTAACGGCTGATGTCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGc                                                                                                               <  8:26763/149‑1 (MQ=255)
                                  aTCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCAcc                                                                                                            >  4:366062/1‑148 (MQ=255)
                                     cccGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGccc                                                                                                                                           >  6:341199/1‑114 (MQ=255)
                                     cccGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGccc                                                                                                                                           <  5:341199/114‑1 (MQ=255)
                                      ccGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAgg                                                                                                                >  7:139117/1‑140 (MQ=255)
                                      ccGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAgg                                                                                                                <  8:139117/140‑1 (MQ=255)
                                      ccGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCt                                                                                                       >  7:212431/1‑149 (MQ=255)
                                        gCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTgg                                                                                                     <  4:99532/149‑1 (MQ=255)
                                         cACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGa                                                                                                    <  8:212431/149‑1 (MQ=255)
                                          aCGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGAt                                                                                                   <  3:87273/149‑1 (MQ=255)
                                           cGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATc                                                                                                  <  5:266309/149‑1 (MQ=255)
                                             ggCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATcc                                                                                                 <  2:204290/149‑1 (MQ=255)
                                              gCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAAt                                                                                                                                       <  7:120284/110‑1 (MQ=255)
                                              gCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAAt                                                                                                                                       >  8:120284/1‑110 (MQ=255)
                                                     gTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACg                                                                                         >  7:72399/1‑149 (MQ=255)
                                                     gTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACg                                                                                         >  2:368749/1‑149 (MQ=255)
                                                                   ccTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCaa                                                                                   >  6:80203/1‑141 (MQ=255)
                                                                   ccTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCaa                                                                                   <  5:80203/141‑1 (MQ=255)
                                                                   ccTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCg                                                                            >  7:15434/1‑148 (MQ=255)
                                                                                 cccACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTaaa                                                             >  4:192324/1‑149 (MQ=255)
                                                                                  ccACGAAGGCACGGTGTATGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAAc                                                            <  8:271185/149‑1 (MQ=255)
                                                                                  ccACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAAc                                                            >  7:260469/1‑149 (MQ=255)
                                                                                  ccACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAAc                                                            >  8:331068/1‑149 (MQ=255)
                                                                                  ccACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAAc                                                            <  8:72399/149‑1 (MQ=255)
                                                                                    aCGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAAcc                                                           >  8:39316/1‑148 (MQ=255)
                                                                                     cGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGt                                                         <  6:99879/149‑1 (MQ=255)
                                                                                      gAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTGGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTa                                                        <  6:91974/149‑1 (MQ=255)
                                                                                           cACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCTAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCg                                                                            <  4:342650/124‑1 (MQ=255)
                                                                                           cACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCAAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCg                                                                            >  3:342650/1‑124 (MQ=255)
                                                                                                   aCGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTaaa                                                             <  2:212160/131‑1 (MQ=255)
                                                                                                   aCGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTaaa                                                             >  1:212160/1‑131 (MQ=255)
                                                                                                    cGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTg                                                                                                      >  4:86929/1‑89 (MQ=255)
                                                                                                       cGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGca                                       <  8:260469/149‑1 (MQ=255)
                                                                                                            gATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGcc                                  >  3:329847/1‑149 (MQ=255)
                                                                                                                 tgGGAGCGCCCTTGTGTCGCGATTGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCGCTCTAAACCGTACCTTTCGCCAGTAGGCACAGcccat                               <  4:326870/148‑3 (MQ=255)
                                                                                                                 tgGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCtt                               <  8:213369/148‑1 (MQ=255)
                                                                                                                 tgGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTg                              <  6:325746/149‑1 (MQ=255)
                                                                                                                     aGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAAc                          >  2:107492/1‑149 (MQ=255)
                                                                                                                          ccTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTAc                                                       <  1:252327/115‑1 (MQ=255)
                                                                                                                          ccTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTAc                                                       >  2:252327/1‑115 (MQ=255)
                                                                                                                          ccTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAgg                                         <  3:103903/129‑1 (MQ=255)
                                                                                                                          ccTTGTGTCGAGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAgg                                         >  4:103903/1‑129 (MQ=255)
                                                                                                                                 tCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCt              >  1:214031/1‑149 (MQ=255)
                                                                                                                                    cgATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTa                                           >  7:351283/1‑117 (MQ=255)
                                                                                                                                    cgATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTa                                           <  8:351283/117‑1 (MQ=255)
                                                                                                                                    cgATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGc           <  7:39316/149‑1 (MQ=255)
                                                                                                                                     ttggggCCGCGGTGCGGCCCCAATCCTGTGCACCCAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGcc          >  7:210684/4‑149 (MQ=255)
                                                                                                                                      aTGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGccc         <  7:331068/149‑1 (MQ=255)
                                                                                                                                         ggCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGCCCaa       >  1:364743/1‑148 (MQ=255)
                                                                                                                                           ccGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGCCCAacac    <  3:366062/149‑1 (MQ=255)
                                                                                                                                           ccGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGCCCAacac    >  8:169251/1‑149 (MQ=255)
                                                                                                                                             gCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGCCCAACACCa  >  4:326250/1‑149 (MQ=255)
                                                                                                                                             gCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGCCCAACACCa  <  2:226081/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GCGGTGACCTACCAGCTCGTACGCAGCAGGCGCCATCCCCGCACGGCCGTGCGTAACGGCTGATGGCCTCGGGGGCTGTTCCCACGAAGGCACGGTGTACGCCGTCCGATGTGGGAGCGCCCTTGTGTCGCGATGGGCCGCACCGCGGCCCCAATCCTGTGCACACAAATCCATCAGGCACCAACCTGGATCCAACCAACGCCCCAAGACCTCGCACAGCCCCTCTAAACCGTACCTTTCGCCACTAGGCACAGCCCTTGCAACGCCCTCCCCGCTTGCCCAACACCA  >  NC_002947/5508084‑5508371

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