Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,551,701 A→T 6.3% D357V (GAC→GTC) ‡ 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,7010AT 6.3% 187.7 / 2.7 64D357V (GAC→GTC) ‡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.03e-01

GCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATGGAGGCACAGGATACCG  >  NC_002947/4551552‑4551847
                                                                                                                                                     |                                                                                                                                                  
gCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCA‑CGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGa                                                                                                                                                    >  7:180622/1‑149 (MQ=255)
      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)
                                       ggTGGGGCATACCA‑CGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGa                                                                                                             <  8:180622/149‑1 (MQ=255)
                                            ggCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGc                                                                                                         <  3:138413/149‑1 (MQ=255)
                                             gCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCa                                                                                                        <  6:105308/149‑1 (MQ=255)
                                             gCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCa                                                                                                        <  4:37214/149‑1 (MQ=255)
                                                   caccGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATg                                                                                                  >  1:295881/1‑149 (MQ=255)
                                                    accGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGTGTCCTCGCAc                                                                                                                                          >  6:212544/1‑108 (MQ=255)
                                                    accGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGTGTCCTCGCAc                                                                                                                                          <  5:212544/108‑1 (MQ=255)
                                                    accGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGc                                                                                                 <  2:153169/149‑1 (MQ=255)
                                                           ggAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGc                                                                                                                                     <  1:33954/106‑1 (MQ=255)
                                                           ggAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGc                                                                                                                                     >  2:33954/1‑106 (MQ=255)
                                                           ggAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAgt                                                                                          >  7:55175/1‑149 (MQ=255)
                                                            gAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAgtg                                                                                         <  2:295881/149‑1 (MQ=255)
                                                               aTACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGa                                                                                      >  1:157928/1‑149 (MQ=255)
                                                                   cTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCg                                                                                   >  8:193096/1‑148 (MQ=255)
                                                                   cTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCg                                                                                   >  5:94150/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                                                            >  1:189321/1‑149 (MQ=255)
                                                                                         tcgCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAAc                                                            <  8:55175/149‑1 (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                  >  1:140125/1‑147 (MQ=255)
                                                                                                                                     cacaTTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATg                  >  2:241548/1‑147 (MQ=255)
                                                                                                                                      acaTTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATGGAg               >  5:123140/1‑149 (MQ=255)
                                                                                                                                          ctggtgcgAGGTGTCCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCa                          >  3:29503/4‑134 (MQ=255)
                                                                                                                                          ctggtgcgAGGTGTCCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCa                          <  4:29503/131‑1 (MQ=255)
                                                                                                                                            cgtgcgAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATGGAGGCACAg         <  3:73337/149‑1 (MQ=255)
                                                                                                                                              tgcgAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCt                                                                 <  2:217/91‑1 (MQ=255)
                                                                                                                                              tgcgAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCt                                                                 >  1:217/1‑91 (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)
                                                                                                                                                     |                                                                                                                                                  
GCATGGGCATGCCACTGTCGAAGCGTTCCGCGGGCAGGCGGTGGGGCATACCACCGACTGGAAATACCTGGACATCAACTACAAGTCGGTCGCGCACATTGACCAGGAGGCGTGCATTGGCTGCGGGCGCTGCCACATTGCGTGCGAGGACACCTCGCACCAGGCGATTGCCAGCACGCTGAAGGCTGACGGCACGCATGCCTACAGTGTGATCGAAGAGGAATGCGTAGGCTGTAACCTGTGCCAGATCACCTGCCCGGTGGAAAACTGCATCGAAATGGAGGCACAGGATACCG  >  NC_002947/4551552‑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: ‑