Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,725,646 T→G 6.3% K440N (AAA→AAC opgH ← glucans biosynthesis glucosyltransferase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,725,6460TG 6.3% 245.3 / 6.8 79K440N (AAA→AACopgHglucans biosynthesis glucosyltransferase H
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base G (3/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

GGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTG  >  NC_002947/5725498‑5725792
                                                                                                                                                     |                                                                                                                                                   
ggTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCt                                                                                                                                                     <  8:129646/149‑1 (MQ=255)
  tCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCTCACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCttt                                                                                                                                                   >  3:31277/1‑149 (MQ=255)
  tCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCttt                                                                                                                                                   >  4:197649/1‑149 (MQ=255)
  tCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCttt                                                                                                                                                   >  2:178303/1‑149 (MQ=255)
     cgcTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGcc                                                                                                                                                <  4:79358/149‑1 (MQ=255)
       cTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCg                                                                                                                                              <  4:335166/149‑1 (MQ=255)
       cTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCg                                                                                                                                              >  5:269081/1‑149 (MQ=255)
            gCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCgg                                                                                                                                             >  4:320701/1‑145 (MQ=255)
            gCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCgg                                                                                                                                             <  3:320701/145‑1 (MQ=255)
              tCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAATCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACgg                                                                                                                                       >  1:178201/1‑149 (MQ=255)
              tCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGACAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCACCgg                                                                                                                                       >  6:94756/1‑149 (MQ=255)
                gTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACgg                                                                                                                                       <  7:244919/147‑1 (MQ=255)
                 tCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCGGGCAAAGgcgc                                                                                                                                    <  7:116630/149‑1 (MQ=255)
                    aGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAg                                                                                                                                 >  2:200297/1‑149 (MQ=255)
                     gCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAgg                                                                                                                                <  3:197649/149‑1 (MQ=255)
                       aGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGc                                                                                                                                                 >  6:116137/1‑130 (MQ=255)
                       aGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGc                                                                                                                                                 <  5:116137/130‑1 (MQ=255)
                         gTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGgcg                                                                                                                             >  7:115033/1‑148 (MQ=255)
                            tGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAATCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAg                                                                                                                          <  2:178201/148‑1 (MQ=255)
                                 gcAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTc                                                                                                                    >  6:109682/1‑148 (MQ=255)
                                   aGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCGACACACCGCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAAGGCGCCTTTGCCCGGCAACGGGGCCAGGGCGCAGTGCTCGAt                                                                                                                 <  2:230821/149‑1 (MQ=255)
                                    gTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATg                                                                                                                <  6:158144/149‑1 (MQ=255)
                                      tcttcGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGaa                                                                                                              >  3:52472/1‑149 (MQ=255)
                                          cGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACgg                                                                                                           <  2:71203/148‑1 (MQ=255)
                                              gCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTc                                                                                                       <  7:124266/148‑1 (MQ=255)
                                               cTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCTGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCAt                                                                                                     <  4:261204/149‑1 (MQ=255)
                                                tGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCGGGCAAAGgc                                                                                                                                      <  7:138641/116‑1 (MQ=255)
                                                tGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCGGGCAAAGgc                                                                                                                                      >  8:138641/1‑116 (MQ=255)
                                                           tCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATcgc                                                                                         >  5:126231/1‑149 (MQ=255)
                                                            cGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATcgc                                                                                         <  7:323823/148‑1 (MQ=255)
                                                             gTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGt                                                                                       >  3:321224/1‑149 (MQ=255)
                                                             gTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGt                                                                                       >  6:190297/1‑149 (MQ=255)
                                                               aGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTtg                                                                                     >  3:208925/1‑149 (MQ=255)
                                                                 acAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTtgtg                                                                                   <  5:94756/148‑1 (MQ=255)
                                                                  cAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGTAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTgg                                                                                  <  2:363860/149‑1 (MQ=255)
                                                                  cAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTgg                                                                                  >  3:337806/1‑149 (MQ=255)
                                                                   aaTCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGc                                                                                 >  3:15525/1‑149 (MQ=255)
                                                                     tCCTCACACCCCAGCCGGCACGGCGCATCAGCACAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGccc                                                                               <  4:31277/149‑1 (MQ=255)
                                                                        acacacCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGccc                                                                               >  6:58703/1‑146 (MQ=255)
                                                                          acacCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCag                                                                            >  5:207112/1‑147 (MQ=255)
                                                                          acacCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCag                                                                            >  8:126600/1‑147 (MQ=255)
                                                                          acacCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGCATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCagta                                                                          <  5:117089/149‑1 (MQ=255)
                                                                            acCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCagtagt                                                                        <  6:269081/149‑1 (MQ=255)
                                                                              cccAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGgcg                                                                                                                             <  3:61161/95‑1 (MQ=255)
                                                                              cccAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGgcg                                                                                                                             >  4:61161/1‑95 (MQ=255)
                                                                              cccAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGc                                                                      <  1:200297/149‑1 (MQ=255)
                                                                               ccAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAgg                                                                                                                                <  4:34092/91‑1 (MQ=255)
                                                                               ccAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAgg                                                                                                                                >  3:34092/1‑91 (MQ=255)
                                                                                   ccGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTc                                                                 <  8:115033/149‑1 (MQ=255)
                                                                                    cGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCa                                                                <  5:109682/149‑1 (MQ=255)
                                                                                             gcATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGc                                                       <  3:73447/149‑1 (MQ=255)
                                                                                               aTCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCa                                                     <  5:141980/149‑1 (MQ=255)
                                                                                                  aGCTCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGaa                                                  <  6:223406/149‑1 (MQ=255)
                                                                                                   gcgcAGCTTCGACGAAGTCGTGGGTAAAGGATGGCACCGCCCAACGCGCATTTCCCCGCCAACGGCGCAAGGGCGCATTCATCGATGACCGGGTACATCCGGATGATCCCGTGTTGGCCCAAGAAGTGCGAc                                                                   >  7:53937/1‑132 (MQ=255)
                                                                                                   gcgcAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGTCTTTGCCCGGCAACGGCGCGAGGTCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGAc                                                                   <  8:53937/131‑1 (MQ=255)
                                                                                                   gcgcAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGAc                                                                   <  2:15361/131‑1 (MQ=255)
                                                                                                   gcgcAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGAc                                                                   >  1:15361/1‑131 (MQ=255)
                                                                                                    cgcAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGt                                                <  4:337806/149‑1 (MQ=255)
                                                                                                     gcAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTg                                               <  4:321224/149‑1 (MQ=255)
                                                                                                                gaAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTtgt                                                                                    <  3:168709/101‑1 (MQ=255)
                                                                                                                gaAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTtgt                                                                                    >  4:168709/1‑101 (MQ=255)
                                                                                                                       tGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAg                             <  4:52472/149‑1 (MQ=255)
                                                                                                                         ggAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTg                           >  3:265813/1‑149 (MQ=255)
                                                                                                                          gAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGcc                                                                                <  4:316340/95‑1 (MQ=255)
                                                                                                                          gAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGcc                                                                                >  3:316340/1‑95 (MQ=255)
                                                                                                                          gAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTgg                          >  5:4050/1‑149 (MQ=255)
                                                                                                                          gAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTgg                          >  4:27002/1‑149 (MQ=255)
                                                                                                                            aaaGGATGGCACCGGCGAACGCGCCTTTGCGCGGCAACGGCGCCAGGGCGCAGTGCGTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTgg                          >  1:101168/1‑149 (MQ=255)
                                                                                                                            aaaGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTggg                         <  3:27002/149‑1 (MQ=255)
                                                                                                                               ggATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCg                      <  4:265813/149‑1 (MQ=255)
                                                                                                                                     cACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACCcg                >  7:149833/1‑149 (MQ=255)
                                                                                                                                     cACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACAcg                >  3:46447/1‑149 (MQ=255)
                                                                                                                                      aCCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACAcgc               >  5:39798/1‑149 (MQ=255)
                                                                                                                                      aCCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACAcgc               >  1:114137/1‑149 (MQ=255)
                                                                                                                                      aCCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACAcgc               <  6:4050/149‑1 (MQ=255)
                                                                                                                                           ccctggcgcCGTTGCCGGGCAAAGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCggg            >  2:348675/6‑147 (MQ=255)
                                                                                                                                             aaCGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGc        >  4:188087/1‑149 (MQ=255)
                                                                                                                                             aaCGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGc        >  6:94744/1‑149 (MQ=255)
                                                                                                                                              tggcgcCGTTGCCGGGCAAAGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCTACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGc        >  8:306101/3‑148 (MQ=255)
                                                                                                                                                gcgcCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGAATCACCCAACTGCCAGAAGTGCAGACCGGCa                                     >  8:150379/1‑117 (MQ=255)
                                                                                                                                                 cgcCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAAc    <  5:190297/149‑1 (MQ=255)
                                                                                                                                                   ccTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAActg  <  4:46447/149‑1 (MQ=255)
                                                                                                                                                   ccTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAActg  <  1:178303/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGTCGCGCTTGAGCTCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTG  >  NC_002947/5725498‑5725792

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: