Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,019,083 G→C 7.1% H39Q (CAC→CAG)  dadX ← alanine racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,019,0830GC 7.1% 265.4 / 8.9 84H39Q (CAC→CAG) dadXalanine racemase
Reads supporting (aligned to +/- strand):  ref base G (38/40);  new base C (3/3);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

GCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAAT  >  NC_002947/6018935‑6019231
                                                                                                                                                    |                                                                                                                                                    
gcgcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCg                                                                                                                                                      >  4:63933/1‑149 (MQ=255)
  gcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTg                                                                                                                                                    <  4:307211/149‑1 (MQ=255)
     aGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGcc                                                                                                                                                 <  8:115507/149‑1 (MQ=255)
          gCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAgg                                                                                                                                            <  2:221942/149‑1 (MQ=255)
          gCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAgg                                                                                                                                            >  8:181233/1‑149 (MQ=255)
          gCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAg                                                                                                                                             <  6:210628/148‑1 (MQ=255)
             ggCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCAt                                                                                                                                         >  7:127319/1‑149 (MQ=255)
                ttCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCgg                                                                                                                                      >  5:91194/1‑149 (MQ=255)
                ttCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCgg                                                                                                                                      >  6:129160/1‑149 (MQ=255)
                  cAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGc                                                                                                                                     <  6:117520/148‑1 (MQ=255)
                   aaaaaaaCCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCt                                                                                                                                   >  8:104599/1‑149 (MQ=255)
                    aaaaaaCCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGCCATCGGCCtt                                                                                                                                  >  7:128982/1‑149 (MQ=255)
                     aaaaaCCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATACAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCtt                                                                                                                                  >  8:106501/1‑148 (MQ=255)
                       aaaCCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGAt                                                                                                                               <  3:41114/149‑1 (MQ=255)
                         aCCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCa                                                                                                                             >  2:316556/1‑149 (MQ=255)
                              ccAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGc                                                                                                                         <  5:16469/148‑1 (MQ=255)
                                agcagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGcc                                                                                                                                    <  5:95509/135‑1 (MQ=255)
                                agcagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGcc                                                                                                                                    >  6:95509/1‑135 (MQ=255)
                                 gcagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGcc                                                                                                                                    <  2:177312/134‑1 (MQ=255)
                                 gcagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGcc                                                                                                                                    >  1:177312/1‑134 (MQ=255)
                                 gcagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAgg                                                                                                                     >  5:364433/1‑149 (MQ=255)
                                  cagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCACGGCGCCGTGCGCTGCGCCCTGGccctcgag                                                                                                                                           <  3:152649/126‑7 (MQ=37)
                                  cagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCACGGCGCCGTGCGCTGCGCCCTGGccctcgag                                                                                                                                           >  4:152649/1‑120 (MQ=37)
                                  cagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCACGGCGCCGTGCGCTGCGCCCTGGcc                                                                                                                                                 <  7:266885/120‑1 (MQ=38)
                                  cagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCACGGCGCCGTGCGCTGCGCCCTGGcc                                                                                                                                                 >  8:266885/1‑120 (MQ=38)
                                       gccccGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTg                                                                                                               <  1:316556/146‑1 (MQ=255)
                                       gTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTg                                                                                                               <  3:233257/149‑1 (MQ=255)
                                       gTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTg                                                                                                               <  3:209808/149‑1 (MQ=255)
                                          ccGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGcc                                                                                                                                    >  1:185102/1‑125 (MQ=255)
                                          ccGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGcc                                                                                                                                    <  2:185102/125‑1 (MQ=255)
                                           cGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCAc                                                                                                           >  2:17211/1‑149 (MQ=255)
                                             gggCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCtt                                                                                                                 <  5:242255/141‑1 (MQ=255)
                                             gggCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCtt                                                                                                                 >  6:242255/1‑141 (MQ=255)
                                                 tttGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATACAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCa                                                                                                     <  7:106501/149‑1 (MQ=255)
                                                    gATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCtt                                                                                                                 >  7:75262/1‑134 (MQ=255)
                                                    gATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCtt                                                                                                                 <  8:75262/134‑1 (MQ=255)
                                                     aTACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGtt                                                                                                  >  4:247831/1‑148 (MQ=255)
                                                         cggcggCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACgg                                                                                             >  6:275695/1‑149 (MQ=255)
                                                         cggcggCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACgg                                                                                             <  8:286011/149‑1 (MQ=255)
                                                            cggcgCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCa                                                                                                                             <  2:153003/114‑1 (MQ=255)
                                                            cggcgCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCa                                                                                                                             >  1:153003/1‑114 (MQ=255)
                                                                    aGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTa                                                                                  <  5:129160/149‑1 (MQ=255)
                                                                     gTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAg                                                                                 <  3:63933/149‑1 (MQ=255)
                                                                        ccAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTtg                                                                              >  5:302417/1‑149 (MQ=255)
                                                                              ccTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACg                                                                                              <  5:6604/127‑1 (MQ=255)
                                                                              ccTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACg                                                                                              >  6:6604/1‑127 (MQ=255)
                                                                              ccTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGgcg                                                                         >  6:284789/1‑148 (MQ=255)
                                                                                tctTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTgg                                                                                                              <  3:336448/109‑1 (MQ=255)
                                                                                tctTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTgg                                                                                                              >  4:336448/1‑109 (MQ=255)
                                                                                  ttCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAg                                                                      <  6:364433/147‑1 (MQ=255)
                                                                                         cAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGa                                                              <  1:17211/148‑1 (MQ=255)
                                                                                         cAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAg                                                             <  8:127319/149‑1 (MQ=255)
                                                                                          aGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAgg                                                            >  1:301929/1‑149 (MQ=255)
                                                                                            gCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGt                                                           <  6:302417/148‑1 (MQ=255)
                                                                                              cACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGa                                                        >  4:233952/1‑149 (MQ=255)
                                                                                              cACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGa                                                        >  7:76013/1‑149 (MQ=255)
                                                                                               aCGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTc                                                          >  8:232156/1‑146 (MQ=255)
                                                                                               aCGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGAt                                                       >  4:289872/1‑149 (MQ=255)
                                                                                                   cAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAgg                                                   >  3:176696/1‑149 (MQ=255)
                                                                                                    aaaGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATc                                                      <  6:326058/145‑1 (MQ=255)
                                                                                                    aaaGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATc                                                      >  5:326058/1‑145 (MQ=255)
                                                                                                         cGTCGGCTTCGGCCTCGAGGGCGAGGGCGCAGCGCACGGCGCCGTGGCGGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCgcgg                                             <  5:127963/149‑1 (MQ=255)
                                                                                                             ggCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCgcgggcgg                                         >  3:37303/1‑149 (MQ=255)
                                                                                                              gCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCgcgggcgg                                         >  7:309399/1‑148 (MQ=255)
                                                                                                              gCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGcgggcggg                                        <  3:233952/149‑1 (MQ=255)
                                                                                                               cTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGa                                       <  7:104599/149‑1 (MQ=255)
                                                                                                                ttCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGAc                                      <  7:181233/149‑1 (MQ=255)
                                                                                                                     ccTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATa                                 <  4:37303/149‑1 (MQ=255)
                                                                                                                      cTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAg                                                                                                                      >  8:2664/1‑63 (MQ=255)
                                                                                                                      cTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAg                                                                                                                      <  7:2664/63‑1 (MQ=255)
                                                                                                                      cTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAg                                >  6:5398/1‑149 (MQ=255)
                                                                                                                      cTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCATTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAg                                >  1:85904/1‑149 (MQ=255)
                                                                                                                        ctacggCCACGGCGCCGTGCGCTGCGCCCTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCa                                                     <  2:133568/122‑1 (MQ=38)
                                                                                                                        ctacggCCACGGCGCCGTGCGCTGCGCCCTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCa                                                     >  1:133568/5‑126 (MQ=38)
                                                                                                                            ggCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCt                           <  7:232156/148‑1 (MQ=255)
                                                                                                                             gCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGc                         <  3:247831/149‑1 (MQ=255)
                                                                                                                             gCCAGGGCGCAGCGAACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGACAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGCACGAATAGCGGCtca                         >  4:82467/1‑147 (MQ=255)
                                                                                                                              ccAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGcc                        >  8:355888/1‑149 (MQ=255)
                                                                                                                                  ggCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTAt                    <  5:275695/149‑1 (MQ=255)
                                                                                                                                    cgcagcgcaCGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTAtgg                  <  4:176696/149‑1 (MQ=255)
                                                                                                                                      cagcgcaCGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTAtggtg                <  5:5398/149‑1 (MQ=255)
                                                                                                                                            aCGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGc                                                                   <  2:212332/92‑1 (MQ=255)
                                                                                                                                            aCGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGc                                                                   >  1:212332/1‑92 (MQ=255)
                                                                                                                                              ggCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGt        >  2:83682/1‑149 (MQ=255)
                                                                                                                                                cgcCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTg      <  7:355888/149‑1 (MQ=255)
                                                                                                                                                    gTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAAt  <  8:76013/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTGGTTCGGGTTGAAAT  >  NC_002947/6018935‑6019231

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