Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,019,064 G→C 6.8% L46V (CTG→GTG)  dadX ← alanine racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,019,0640GC 6.8% 286.1 / 8.3 88L46V (CTG→GTG) dadXalanine racemase
Reads supporting (aligned to +/- strand):  ref base G (41/41);  new base C (3/3);  total (44/44)
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.07e-01

AGGTCGTGCTCGGCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCC  >  NC_002947/6018916‑6019209
                                                                                                                                                    |                                                                                                                                                 
aGGTCGTGCTCGGCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAg                                                                                                                                                   <  5:20709/149‑1 (MQ=255)
 ggTCGTGCTCGGCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAgg                                                                                                                                                  <  5:167126/149‑1 (MQ=255)
 ggTCGTGCTCGGCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAgg                                                                                                                                                  <  2:72631/149‑1 (MQ=255)
           ggCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGgcgcagcgca                                                                                                                                        <  5:340218/149‑1 (MQ=255)
            gCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCAGCCTCGAGGGCCAGGGCGCAGCGCAc                                                                                                                                       <  1:187959/149‑1 (MQ=255)
              gATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACgg                                                                                                                                     >  5:117520/1‑149 (MQ=255)
                  aGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGcc                                                                                                                                 <  2:337274/149‑1 (MQ=255)
                   gcgcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCg                                                                                                                                >  4:63933/1‑149 (MQ=255)
                    cgcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGgc                                                                                                                                    >  8:29506/1‑144 (MQ=255)
                    cgcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGgc                                                                                                                                    <  7:23904/144‑1 (MQ=255)
                    cgcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGgc                                                                                                                                    <  7:29506/144‑1 (MQ=255)
                    cgcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGgc                                                                                                                                    >  8:23904/1‑144 (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)
                                   ttCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGgcg                                                                                                                                               <  3:37904/118‑1 (MQ=255)
                                   ttCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGgcg                                                                                                                                               >  4:37904/1‑118 (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                                                                                                              >  1:177312/1‑134 (MQ=255)
                                                    gcagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGcc                                                                                                              <  2:177312/134‑1 (MQ=255)
                                                    gcagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAgg                                                                                               >  5:364433/1‑149 (MQ=255)
                                                     cagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCACGGCGCCGTGCGCTGCGCCCTGGccctcgag                                                                                                                     >  4:152649/1‑120 (MQ=37)
                                                     cagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCACGGCGCCGTGCGCTGCGCCCTGGccctcgag                                                                                                                     <  3:152649/126‑7 (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                                                                                           >  6:242255/1‑141 (MQ=255)
                                                                gggCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCtt                                                                                           <  5:242255/141‑1 (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                                                                       <  8:286011/149‑1 (MQ=255)
                                                                            cggcggCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACgg                                                                       >  6:275695/1‑149 (MQ=255)
                                                                               cggcgCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCa                                                                                                       >  1:153003/1‑114 (MQ=255)
                                                                               cggcgCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCa                                                                                                       <  2:153003/114‑1 (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                                                                                        >  4:336448/1‑109 (MQ=255)
                                                                                                   tctTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTgg                                                                                        <  3:336448/109‑1 (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                                  >  7:76013/1‑149 (MQ=255)
                                                                                                                 cACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGa                                  >  4:233952/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                               >  1:133568/5‑126 (MQ=38)
                                                                                                                                           ctacggCCACGGCGCCGTGCGCTGCGCCCTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCa                               <  2:133568/122‑1 (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)
                                                                                                                                                    |                                                                                                                                                 
AGGTCGTGCTCGGCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCC  >  NC_002947/6018916‑6019209

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