Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,551,703 A→T 6.2% T358S (ACC→TCC)  pydA → NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,551,7030AT 6.2% 186.7 / 2.2 64T358S (ACC→TCC) pydANADP‑dependent dihydropyrimidine dehydrogenase subunit PreA
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base T (2/2);  total (33/31)
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.28e-01

GCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATGGAGGCACAGGATACCG  >  NC_002947/4551558‑4551847
                                                                                                                                                 |                                                                                                                                                
gCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCt                                                                                                                                               >  4:23765/1‑149 (MQ=255)
  aTGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCg                                                                                                                                             >  7:16999/1‑149 (MQ=255)
   tGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGTACACCTCGc                                                                                                                                            <  2:136298/149‑1 (MQ=255)
   tGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGc                                                                                                                                            <  2:136296/149‑1 (MQ=255)
     ccACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGTGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCAc                                                                                                                                          <  1:118956/149‑1 (MQ=255)
       aCTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACAcc                                                                                                                                                >  8:25378/1‑141 (MQ=255)
       aCTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACAcc                                                                                                                                                <  7:25378/141‑1 (MQ=255)
       aCTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCAc                                                                                                                                          >  8:217817/1‑147 (MQ=255)
            cGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGa                                                                                                                                   >  2:284162/1‑149 (MQ=255)
              aaGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGAtt                                                                                                                                 <  6:259397/149‑1 (MQ=255)
                 cGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGcc                                                                                                                              <  1:284162/149‑1 (MQ=255)
                  gTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCa                                                                                                                             <  3:23765/149‑1 (MQ=255)
                   ttCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAg                                                                                                                            >  5:105308/1‑149 (MQ=255)
                   ttCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAg                                                                                                                            <  7:217817/149‑1 (MQ=255)
                      cgcgGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCAc                                                                                                                         >  3:37214/1‑149 (MQ=255)
                                cGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCt                                                                                                               >  7:86184/1‑149 (MQ=255)
                                 ggTGGGGCATACCACGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGa                                                                                                             <  8:180622/149‑1 (MQ=255)
                                      ggCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGc                                                                                                         <  3:138413/149‑1 (MQ=255)
                                       gCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCa                                                                                                        <  4:37214/149‑1 (MQ=255)
                                       gCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCa                                                                                                        <  6:105308/149‑1 (MQ=255)
                                             caccGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATg                                                                                                  >  1:295881/1‑149 (MQ=255)
                                              accGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGTGTCCTCGCAc                                                                                                                                          <  5:212544/108‑1 (MQ=255)
                                              accGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGTGTCCTCGCAc                                                                                                                                          >  6:212544/1‑108 (MQ=255)
                                              accGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGc                                                                                                 <  2:153169/149‑1 (MQ=255)
                                                     ggAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGc                                                                                                                                     >  2:33954/1‑106 (MQ=255)
                                                     ggAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGc                                                                                                                                     <  1:33954/106‑1 (MQ=255)
                                                     ggAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAgt                                                                                          >  7:55175/1‑149 (MQ=255)
                                                      gAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAgtg                                                                                         <  2:295881/149‑1 (MQ=255)
                                                         aTACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGa                                                                                      >  1:157928/1‑149 (MQ=255)
                                                             cTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCg                                                                                   >  5:94150/1‑148 (MQ=255)
                                                             cTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCg                                                                                   >  8:193096/1‑148 (MQ=255)
                                                                gACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAaga                                                                               <  5:327942/149‑1 (MQ=255)
                                                                   aTCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAgg                                                                             >  5:114352/1‑148 (MQ=255)
                                                                         tACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCg                                                                       >  4:244154/1‑148 (MQ=255)
                                                                           cAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAg                                                                    <  8:86184/149‑1 (MQ=255)
                                                                           cAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAg                                                                    >  4:120144/1‑149 (MQ=255)
                                                                             aGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGc                                                                  >  3:24372/1‑149 (MQ=255)
                                                                              gtcggtcgCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCt                                                                 <  8:16999/149‑1 (MQ=255)
                                                                                   tcgCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAAc                                                            <  8:55175/149‑1 (MQ=255)
                                                                                   tcgCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAAc                                                            >  1:189321/1‑149 (MQ=255)
                                                                                     gcgcACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGa                                                                                  >  5:65629/1‑125 (MQ=255)
                                                                                     gcgcACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGa                                                                                  <  6:65629/125‑1 (MQ=255)
                                                                                       gcACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCtgt                                                        <  3:127021/149‑1 (MQ=255)
                                                                                        cacaTTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCtgtg                                                       <  2:157928/149‑1 (MQ=255)
                                                                                         acaTTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGc                                                      >  4:220189/1‑149 (MQ=255)
                                                                                          ccttGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAATGCTGACGGCACGCATTCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGcc                                                     <  1:203914/147‑1 (MQ=255)
                                                                                                         gTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCgg                                      >  3:310509/1‑149 (MQ=255)
                                                                                                              ttGGCTGCGGGCGCTGGCACATTGCGTGCGAGGACACCTCGCACCAGGAGATTGCCAGCACGCTGAAGGCTGACGTCACGCATGACTACAGTGTGATCGAAGAGGAAAGCg                                                                       >  6:61328/1‑111 (MQ=255)
                                                                                                                 gCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGaaaa                               <  3:244154/148‑1 (MQ=255)
                                                                                                                    gCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTg                            >  1:184065/1‑148 (MQ=255)
                                                                                                                          gcTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGa                      >  4:73337/1‑148 (MQ=255)
                                                                                                                            tGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAAt                   >  1:207236/1‑149 (MQ=255)
                                                                                                                               cacaTTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATg                  >  2:241548/1‑147 (MQ=255)
                                                                                                                               cacaTTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATg                  >  1:140125/1‑147 (MQ=255)
                                                                                                                                acaTTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATGGAg               >  5:123140/1‑149 (MQ=255)
                                                                                                                                    ctggtgcgAGGTGTCCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCa                          <  4:29503/131‑1 (MQ=255)
                                                                                                                                    ctggtgcgAGGTGTCCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCa                          >  3:29503/4‑134 (MQ=255)
                                                                                                                                      cgtgcgAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATGGAGGCACAg         <  3:73337/149‑1 (MQ=255)
                                                                                                                                        tgcgAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCt                                                                 >  1:217/1‑91 (MQ=255)
                                                                                                                                        tgcgAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCt                                                                 <  2:217/91‑1 (MQ=255)
                                                                                                                                        tgcgAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATGGAGGCACAGGa       >  8:95925/1‑149 (MQ=255)
                                                                                                                                            aGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATGGAGGCACAGGATAcc   <  6:114352/149‑1 (MQ=255)
                                                                                                                                            aGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATGGAGGCACAGGATAcc   <  1:128959/149‑1 (MQ=255)
                                                                                                                                             ggACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATGGAGGCACAGGATACCg  <  4:310509/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATGGAGGCACAGGATACCG  >  NC_002947/4551558‑4551847

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: