Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,565,754 A→T 6.7% C578S (TGC→AGC)  mutL ← DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,565,7540AT 6.7% 221.9 / 5.0 74C578S (TGC→AGC) mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base T (3/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

GTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAG  >  NC_002947/5565607‑5565895
                                                                                                                                                     |                                                                                                                                             
gTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCaa                                                                                                                                              <  5:215010/149‑1 (MQ=255)
         gATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTg                                                                                                                                     >  5:142256/1‑149 (MQ=255)
         gATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTg                                                                                                                                     >  2:273483/1‑149 (MQ=255)
          aTCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGc                                                                                                                                    >  3:183623/1‑149 (MQ=255)
          aTCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGc                                                                                                                                    >  6:34560/1‑149 (MQ=255)
          aTCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGATGGCGCGCACGGCACCGTGGCAAGCCATGGTGc                                                                                                                                    >  4:282311/1‑149 (MQ=255)
             gTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGa                                                                                                                                 <  6:142256/149‑1 (MQ=255)
              tCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGag                                                                                                                                >  3:122842/1‑149 (MQ=255)
                 aGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagcag                                                                                                                             <  5:34560/149‑1 (MQ=255)
                            ggTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCACGCACGGCACCGTGGCaa                                                                                                                                              >  5:243081/1‑121 (MQ=255)
                            ggTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCACGCACGGCACCGTGGCaa                                                                                                                                              <  6:243081/121‑1 (MQ=255)
                             gTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagcag                                                                                                                             <  4:164817/137‑1 (MQ=255)
                             gTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagcag                                                                                                                             >  3:164817/1‑137 (MQ=255)
                              tCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTgc                                                                                                                >  8:102678/1‑149 (MQ=255)
                              tCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGCCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTgc                                                                                                                <  6:247943/149‑1 (MQ=255)
                                cAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTgcgc                                                                                                              <  8:210854/149‑1 (MQ=255)
                                   gTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCTCGGCACCATGGCTTGCCACGGTGCCGtgcgc                                                                                                                             <  2:70377/131‑3 (MQ=38)
                                   gTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCt                                                                                                            >  3:66119/1‑148 (MQ=255)
                                   gTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATATCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCTCGGCACCATGGCTTGCCACGGTGCCGtgcgc                                                                                                                             >  1:70377/1‑129 (MQ=38)
                                     aGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGaa                                                                                                         >  5:17436/1‑149 (MQ=255)
                                     aGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGaa                                                                                                         >  4:310253/1‑149 (MQ=255)
                                      ggACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGaa                                                                                                         <  7:170573/148‑1 (MQ=255)
                                          gacCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCTCGGCACCATGGCTTGCCACGGTGCCgtg                                                                                                                                >  5:173557/1‑119 (MQ=39)
                                          gacCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCTCGGCACCATGGCTTGCCACGGTGCCgtg                                                                                                                                <  6:173557/121‑3 (MQ=39)
                                              aTGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCAc                                                                                                <  4:66119/149‑1 (MQ=255)
                                                    gCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAAGGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGt                                                                                            >  1:96344/1‑149 (MQ=255)
                                                    gCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGc                                                                                          >  8:123969/1‑149 (MQ=255)
                                                     cACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGACTGAATGCGATCACTGGTGcc                                                                                         <  4:183623/149‑1 (MQ=255)
                                                      aCTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCg                                                                                        <  7:123969/149‑1 (MQ=255)
                                                        tGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTa                                                                                      >  7:19724/1‑149 (MQ=255)
                                                         ggCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAtt                                                                                                                       <  3:6843/115‑1 (MQ=255)
                                                         ggCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAtt                                                                                                                       >  4:6843/1‑115 (MQ=255)
                                                          gCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACt                                                                                    >  5:151780/1‑149 (MQ=255)
                                                            cGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTc                                                                                   >  2:146723/1‑148 (MQ=255)
                                                                 cgcTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATa                                                                               <  4:122842/147‑1 (MQ=255)
                                                                  gcTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTg                                                                                           <  4:93985/134‑1 (MQ=255)
                                                                  gcTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTg                                                                                           >  3:93985/1‑134 (MQ=255)
                                                                  gcTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAgg                                                                            >  5:149187/1‑149 (MQ=255)
                                                                      ggTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCgg                                                                        <  1:273483/149‑1 (MQ=255)
                                                                      ggTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCgg                                                                        <  6:149187/149‑1 (MQ=255)
                                                                       gtgtTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCgg                                                                        >  5:282314/1‑148 (MQ=255)
                                                                       gtgtTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGACgtg                                                                       >  1:134127/1‑147 (MQ=255)
                                                                            ttCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCa                                                                  >  1:186143/1‑149 (MQ=255)
                                                                            ttCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCa                                                                  >  1:34438/1‑149 (MQ=255)
                                                                              ccATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGc                                                                       <  5:184540/142‑1 (MQ=255)
                                                                              ccATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGc                                                                       >  6:184540/1‑142 (MQ=255)
                                                                               cATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGc                                                                       >  1:292411/1‑141 (MQ=255)
                                                                               cATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGc                                                                       >  1:297148/1‑141 (MQ=255)
                                                                               cATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGc                                                                       <  2:297148/141‑1 (MQ=255)
                                                                               cATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGc                                                                       <  2:292411/141‑1 (MQ=255)
                                                                                 tGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGc                                                                   >  8:33536/1‑143 (MQ=255)
                                                                                 tGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGc                                                                   <  7:33536/143‑1 (MQ=255)
                                                                                  gTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCt                                                            >  2:2447/1‑149 (MQ=255)
                                                                                  gTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCAGCAAGCACGTCCt                                                            >  8:118973/1‑149 (MQ=255)
                                                                                   tCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCt                                                            <  6:282314/148‑1 (MQ=255)
                                                                                     gcgcAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGc                                                          <  2:193383/148‑1 (MQ=255)
                                                                                     gcgcAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCa                                                         >  6:99269/1‑149 (MQ=255)
                                                                                                cGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTg                                              <  3:310253/149‑1 (MQ=255)
                                                                                                cGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTg                                              <  7:102678/149‑1 (MQ=255)
                                                                                                 gTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTgg                                             <  1:2447/149‑1 (MQ=255)
                                                                                                 gTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTgg                                             <  6:17436/149‑1 (MQ=255)
                                                                                                  ttCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGc                                            <  5:222146/149‑1 (MQ=255)
                                                                                                  ttCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCAGCAAGCACGTCCTGCACCAGGCGGTTGGc                                            <  7:118973/149‑1 (MQ=255)
                                                                                                    cATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCt                                           <  6:54069/148‑1 (MQ=255)
                                                                                                        tcGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCa                                                                      >  7:220572/1‑117 (MQ=255)
                                                                                                        tcGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCa                                                                      <  8:220572/117‑1 (MQ=255)
                                                                                                          gggAATTGCCA‑‑GACGCCGGCTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAg                                       <  3:282311/146‑1 (MQ=255)
                                                                                                           ggAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTg                                   <  2:96344/149‑1 (MQ=255)
                                                                                                            gAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGc                                  <  1:286949/149‑1 (MQ=255)
                                                                                                                        aCGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAgg                        >  4:213888/1‑149 (MQ=255)
                                                                                                                            cGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGGTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCgg                    <  7:111701/149‑1 (MQ=255)
                                                                                                                                  gcgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCt               <  1:146723/148‑1 (MQ=255)
                                                                                                                                   cgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACAAGGCGGTTGGATTCAGCCTGATTGAGAAAGGCCGGGATCTgg             >  8:45587/1‑149 (MQ=255)
                                                                                                                                          ggCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCg      >  4:71158/1‑149 (MQ=255)
                                                                                                                                          ggCACCATGGCTTGCCACGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCg      >  8:34066/1‑149 (MQ=255)
                                                                                                                                            cACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGc     >  8:241685/1‑148 (MQ=255)
                                                                                                                                              ccGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAg  <  6:151780/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCA‑‑GACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAG  >  NC_002947/5565607‑5565895

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: