Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,823,005 T→A 6.0% H750L (CAC→CTC)  mutS ← DNA mismatch repair protein MutS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,823,0050TA 6.0% 208.4 / 2.8 67H750L (CAC→CTC) mutSDNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base T (29/34);  new base A (2/2);  total (31/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 = 9.54e-01

CAGGCTGGTGGTTTCCAGCCGGCCAAGGTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTGGGCCAGGCGCTCAGCGGCAGCCCAGGCCAGCGAC  >  NC_002947/1822857‑1823142
                                                                                                                                                    |                                                                                                                                          
cAGGCTGGTGGTTTCCAGCCGGCCAAGGTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATtgt                                                                                                                                            <  3:276026/149‑1 (MQ=255)
cAGGCTGGTGGTTTCCAGCCGGCCAAGGTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATtgt                                                                                                                                            <  1:336715/149‑1 (MQ=255)
     tggtggTTTCCAGCCGGCCAAGGTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCg                                                                                                                                       >  6:174290/1‑149 (MQ=255)
      ggtggtTTCCAGCCGGCCAAGGTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCgg                                                                                                                                      >  3:201001/1‑149 (MQ=255)
                 gCCGGCCAAGGTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAgg                                                                                                                           <  6:351095/149‑1 (MQ=255)
                  ccGGCCAAGGTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGt                                                                                                                          <  7:124790/149‑1 (MQ=255)
                  ccGGCCAAGGTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGt                                                                                                                          <  5:346978/149‑1 (MQ=255)
                      ccAAGGTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCAc                                                                                                                      <  8:56738/149‑1 (MQ=255)
                      ccAAGGTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCAc                                                                                                                      <  8:50407/149‑1 (MQ=255)
                          ggTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTg                                                                                                                                    <  8:156467/131‑1 (MQ=255)
                          ggTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTg                                                                                                                                    >  7:156467/1‑131 (MQ=255)
                          ggTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTg                                                                                                                  <  6:159533/149‑1 (MQ=255)
                          ggTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTg                                                                                                                  <  5:174290/149‑1 (MQ=255)
                              cTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGcca                                                                                                              <  2:236895/149‑1 (MQ=255)
                               tCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGtt                                                                                                                   >  5:346197/1‑143 (MQ=255)
                               tCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGtt                                                                                                                   <  6:346197/143‑1 (MQ=255)
                                     gCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCgg                                                                                                       >  8:169244/1‑149 (MQ=255)
                                     gCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCgg                                                                                                       <  1:40513/149‑1 (MQ=255)
                                         gCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGgctcg                                                                                                   >  3:89431/1‑149 (MQ=255)
                                         gCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGgctcg                                                                                                   <  1:5812/149‑1 (MQ=255)
                                           tGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGgctcgct                                                                                                >  3:379413/1‑149 (MQ=255)
                                            ggATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGgctcgctc                                                                                               >  4:152131/1‑149 (MQ=255)
                                                     ccGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGATTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGTCCACCAGCGGCTCGCTCTCCGgcag                                                                                       >  3:218795/1‑148 (MQ=255)
                                                          ggTAGGCCGGCCAACTGCGCAACGGACAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATTCGTTCATTGTGCTCGGTGGCGTTCATGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGt                                                                                 <  4:218795/149‑1 (MQ=255)
                                                          ggTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGt                                                                                 <  4:67261/149‑1 (MQ=255)
                                                            tACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCa                                                                               <  3:152131/149‑1 (MQ=255)
                                                               gCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGtt                                                                            >  5:95227/1‑149 (MQ=255)
                                                                    ccAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAag                                                                       <  4:255093/149‑1 (MQ=255)
                                                                    ccAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGAGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCa                                                                               <  7:141989/141‑1 (MQ=255)
                                                                    ccAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGAGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCa                                                                               <  7:147613/141‑1 (MQ=255)
                                                                    ccAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGAGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCa                                                                               >  8:141989/1‑141 (MQ=255)
                                                                    ccAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGAGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCa                                                                               >  8:147613/1‑141 (MQ=255)
                                                                      aaCTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAag                                                                       >  6:78775/1‑147 (MQ=255)
                                                                      aaCTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAag                                                                       >  8:186702/1‑147 (MQ=255)
                                                                        cTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAagtagt                                                                   <  4:379413/149‑1 (MQ=255)
                                                                           cgcAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGACTCTCCGGCAGCACGGTCAGTTCGAAGTAgtgg                                                                 <  1:55417/149‑1 (MQ=255)
                                                                           cgcAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAgtgg                                                                 <  8:223930/148‑1 (MQ=255)
                                                                             cAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAgtgggtg                                                              <  4:201001/149‑1 (MQ=255)
                                                                              aaCGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGgca                                                                                        <  1:171603/122‑1 (MQ=255)
                                                                              aaCGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGgca                                                                                        >  2:171603/1‑122 (MQ=255)
                                                                                  gCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGaa                                                         <  7:169244/149‑1 (MQ=255)
                                                                                   ccAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGaa                                                         >  2:276496/1‑148 (MQ=255)
                                                                                             gCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAg                                               <  1:252079/148‑1 (MQ=255)
                                                                                              cTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTTAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAgg                                              >  5:420381/1‑148 (MQ=255)
                                                                                              cTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAgg                                              >  1:77535/1‑148 (MQ=255)
                                                                                                  ggctggcAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACg                                          <  7:9902/148‑1 (MQ=255)
                                                                                                  ggctggcAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACg                                          >  8:138222/1‑148 (MQ=255)
                                                                                                    ctggcAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAg                                       >  6:242783/1‑149 (MQ=255)
                                                                                                        cAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTgg                                   <  4:89431/149‑1 (MQ=255)
                                                                                                         aGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGt                                                                             >  3:325704/1‑106 (MQ=255)
                                                                                                         aGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGt                                                                             <  4:325704/106‑1 (MQ=255)
                                                                                                         aGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTggg                                  >  7:285224/1‑149 (MQ=255)
                                                                                                          gggCCGGGTAGCACGTGGTGCAGGAACACGATGGGTTCATTGTGCTCGGTGGCGTTCAGGTGCTCGTTGGGCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTGGGc                                 <  4:333780/149‑1 (MQ=255)
                                                                                                             ccGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCCCGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTGGGCCAg                              >  2:330976/1‑149 (MQ=255)
                                                                                                             ccGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTGGGctgc                              >  3:2736/1‑146 (MQ=255)
                                                                                                               gggTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTGGGCCAGgc                            <  7:186702/149‑1 (MQ=255)
                                                                                                                ggTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTGGGCCAGgcg                           <  2:77535/149‑1 (MQ=255)
                                                                                                                   aGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGtt                                                                            <  3:262489/97‑1 (MQ=255)
                                                                                                                   aGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGtt                                                                            >  4:262489/1‑97 (MQ=255)
                                                                                                                     cACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGaa                                                                        <  6:228860/99‑1 (MQ=255)
                                                                                                                     cACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGaa                                                                        >  5:228860/1‑99 (MQ=255)
                                                                                                                        gtggtgCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGgcagc                                                                                      >  2:184746/1‑82 (MQ=255)
                                                                                                                        gtggtgCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGgcagc                                                                                      <  1:184746/82‑1 (MQ=255)
                                                                                                                           gtgCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAg                                                     >  6:237722/1‑112 (MQ=255)
                                                                                                                           gtgCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAACTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTGGGCCAGGCGCTCAGCGGCAg                >  3:272334/1‑149 (MQ=255)
                                                                                                                                 gAAGACGATGCGTTCATTGTGCTCGGTGGCGGTCAGGTGCACGTTGGCCTCCAGCGGCTCGCTCTCCGGCATCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTGGGCCAGGCGCTCAGCGGCAGCCCa            <  4:272334/147‑1 (MQ=255)
                                                                                                                                     acGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTGGGCCAGGCGCTCAGCGGCAGCCCAGGCCAg      >  2:359777/1‑149 (MQ=255)
                                                                                                                                     acGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTGGGCCAGGCGCTCAGCGGCAGCCCAGGCCAg      >  3:342773/1‑149 (MQ=255)
                                                                                                                                          gCGTTCATTGTGCTCGGTGGCGTTTAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTGGGCCAGGCGCTCAGCGGCAGCCCAGGCCAGCGAc  <  6:420381/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
CAGGCTGGTGGTTTCCAGCCGGCCAAGGTGCTCGCGGGCACGCTGGATCACTGCCGTCGGTACGCCGGCCAACTGCGCAACGGCCAGGCCGTAGCTTTGGCTGGCAGGGCCGGGTAGCACGTGGTGCAGGAACACGATGCGTTCATTGTGCTCGGTGGCGTTCAGGTGCACGTTGGCCACCAGCGGCTCGCTCTCCGGCAGCACGGTCAGTTCGAAGTAGTGGGTGGCGAACAAGGTGTAGGCACGCAGCTGGGCCAGGCGCTCAGCGGCAGCCCAGGCCAGCGAC  >  NC_002947/1822857‑1823142

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: