Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,019,070 A→C 7.0% C44G (TGC→GGC)  dadX ← alanine racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,019,0700AC 7.0% 269.9 / 10.3 86C44G (TGC→GGC) dadXalanine racemase
Reads supporting (aligned to +/- strand):  ref base A (40/40);  new base C (3/3);  total (43/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.30e-01

GGCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTG  >  NC_002947/6018927‑6019217
                                                                                                                                               |                                                                                                                                                   
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                                                                                                                                            >  8:23904/1‑144 (MQ=255)
         cgcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGgc                                                                                                                                            <  7:23904/144‑1 (MQ=255)
         cgcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGgc                                                                                                                                            <  7:29506/144‑1 (MQ=255)
          gcCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTg                                                                                                                                      <  4:307211/149‑1 (MQ=255)
             aGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGcc                                                                                                                                   <  8:115507/149‑1 (MQ=255)
                  gCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAgg                                                                                                                              >  8:181233/1‑149 (MQ=255)
                  gCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAgg                                                                                                                              <  2:221942/149‑1 (MQ=255)
                  gCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAg                                                                                                                               <  6:210628/148‑1 (MQ=255)
                     ggCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCAt                                                                                                                           >  7:127319/1‑149 (MQ=255)
                        ttCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCgg                                                                                                                        >  6:129160/1‑149 (MQ=255)
                        ttCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCgg                                                                                                                        >  5:91194/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                                                                                                                                   >  8:266885/1‑120 (MQ=38)
                                          cagcagTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCACGGCGCCGTGCGCTGCGCCCTGGcc                                                                                                                                   <  7:266885/120‑1 (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                                                                               >  6:275695/1‑149 (MQ=255)
                                                                 cggcggCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACgg                                                                               <  8:286011/149‑1 (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                                                                                                        <  7:2664/63‑1 (MQ=255)
                                                                                                                              cTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAg                                                                                                        >  8:2664/1‑63 (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)
                                                                                                                                               |                                                                                                                                                   
GGCGATCAGCGCCAGTTCGCTGGCTTCAAAAAAACCTTCCAGCAGCAGTACCGGGGCTTTGATACCGGCGGCGCGCAGTTCCAGCGCCTCTTCGATGCAGGCCACGGCAAAGCCGTCGGCTTCGGCCTCGAGGGCCAGGGCGCAGCGCACGGCGCCGTGGCCGTAGGCATCGGCCTTGATCACGGCGAGGGCTTTGGCACCCGTCAGTTCACGGGCCAGACGGTAGTTGTGGCGCAGGGCCTGGAGGTCGATCAGGGCGCGGGCGGGACGCATAGCGGCTGCCTTATGGTG  >  NC_002947/6018927‑6019217

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