Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,280,813 A→G 6.1% E317G (GAA→GGA)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,280,8130AG 6.1% 200.5 / 2.7 66E317G (GAA→GGA) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base G (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGC  >  NC_000913/1280672‑1280961
                                                                                                                                             |                                                                                                                                                    
taGTGCGTTACATAGTGTCCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCATAACTGTGCGTTCTGTGGCTGGCACAGAAACCGCGCACCGATGCGGCAATGTCGCTGGCGATGTGCCACGTAATTCTGCGTTAATTcc                                                                                                                                                 <  4:362124/146‑1 (MQ=255)
   tGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc                                                                                                                                            <  8:323598/149‑1 (MQ=255)
     cGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCg                                                                                                                                           <  1:16976/148‑1 (MQ=255)
       ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                                                                                                                                         >  3:309145/1‑148 (MQ=255)
       ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                                                                                                                                         >  1:1927/1‑148 (MQ=255)
          cAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCAGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAAcc                                                                                                                                     >  4:324528/1‑149 (MQ=255)
           aaaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAAcc                                                                                                                                     >  8:387262/1‑148 (MQ=255)
           aaaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAAcc                                                                                                                                     >  7:369446/1‑148 (MQ=255)
            aaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAAccc                                                                                                                                    >  1:90080/1‑148 (MQ=255)
             aGGGACCTAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCaa                                                                                                                                  <  6:370732/149‑1 (MQ=255)
             aGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCaa                                                                                                                                  <  4:83771/149‑1 (MQ=255)
               ggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc                                                                                                                                <  4:309145/149‑1 (MQ=255)
                 aCCAAAACTGTTGCCGTCACAACAGACTACGCTGAAATCTCCAAACTGTGGGATCTGTGGCTGGCACCGAAACAGCGCACCGATGCGGCAATGGCTCTGGAGATGGGACACGTAAAGCTGCGTGAATTACACGTCGACAACCCAAGc                                                                                                                                >  2:99641/1‑147 (MQ=255)
                   cAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTCGCCCCGAAACAGGGCACCGATGCGGCAATCGCGCTGGCGATGCGCCACGTAATGCAGCGTGAATTCCCCCTCGACAACCCAAGCCac                                                                                                                             >  2:261410/1‑147 (MQ=255)
                                 tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGt                                                                                                                            <  4:131935/135‑1 (MQ=255)
                                 tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGt                                                                                                                            >  3:131935/1‑135 (MQ=255)
                                     acCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGc                                                                                                          >  5:370443/1‑149 (MQ=255)
                                       cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGt                                                                                                        >  7:221020/1‑149 (MQ=255)
                                          aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGc                                                                                                     <  2:45841/149‑1 (MQ=255)
                                                 tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGGATTCCACCTCGACAACCCAAGCCAg                                                                                                                             >  7:57131/1‑118 (MQ=255)
                                                 tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGGATTCCACCTCGACAACCCAAGCCAg                                                                                                                             >  7:63561/1‑118 (MQ=255)
                                                 tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGGATTCCACCTCGACAACCCAAGCCAg                                                                                                                             <  8:57131/118‑1 (MQ=255)
                                                 tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGGATTCCACCTCGACAACCCAAGCCAg                                                                                                                             <  8:63561/118‑1 (MQ=255)
                                                 tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACAcc                                                                                               >  4:70465/1‑148 (MQ=255)
                                                        gCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGc                                                                                        <  8:369446/148‑1 (MQ=255)
                                                          cAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCCCCTCGACAACCCAAGCCAGTATTTCACCGCCTATGTGCGTCGCTACCCCGACATGCCGa                                                                                     >  5:166173/1‑149 (MQ=255)
                                                                 tgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTgg                                                                               <  1:145534/148‑1 (MQ=255)
                                                                 tgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATATCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTgg                                                                               <  2:90080/148‑1 (MQ=255)
                                                                   cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGa                                                                            <  4:233501/149‑1 (MQ=255)
                                                                        tgtgGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCCCCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTg                                                                       >  5:78340/1‑149 (MQ=255)
                                                                        tgtgGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCt                                                                        >  8:46592/1‑148 (MQ=255)
                                                                          tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCt                                                                        <  7:46592/146‑1 (MQ=255)
                                                                          tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTgg                                                                      <  5:19952/148‑1 (MQ=255)
                                                                           ggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTgg                                                                      >  8:172241/1‑147 (MQ=255)
                                                                            gctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGgaa                                                                    <  1:261410/148‑1 (MQ=255)
                                                                            gctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGgaa                                                                    <  7:387262/148‑1 (MQ=255)
                                                                                gcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCAc                                                                                                                     >  3:90284/1‑95 (MQ=255)
                                                                                gcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCAc                                                                                                                     <  4:90284/95‑1 (MQ=255)
                                                                                    cGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGGCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGAc                                                           <  1:23528/149‑1 (MQ=255)
                                                                                     gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACg                                                          >  2:85480/1‑149 (MQ=255)
                                                                                      aaaCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACgg                                                         >  5:101106/1‑149 (MQ=255)
                                                                                           gggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTtac                                                     >  3:151874/1‑148 (MQ=255)
                                                                                               aCCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGAAGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGc                                                >  1:264901/1‑149 (MQ=255)
                                                                                               aCCGATGCGGCAATGGCGCCGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGc                                                <  4:105634/149‑1 (MQ=255)
                                                                                                 cGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                                                                                             <  5:247231/102‑1 (MQ=255)
                                                                                                 cGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGa                                                                                             >  6:247231/1‑102 (MQ=255)
                                                                                                  gATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGc                                             <  1:85480/149‑1 (MQ=255)
                                                                                                   aTGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTg                                              >  2:10185/1‑147 (MQ=255)
                                                                                                    tGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAg                                           <  4:124418/149‑1 (MQ=255)
                                                                                                           aTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTAt                                    >  4:259241/1‑149 (MQ=255)
                                                                                                                gcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCa                                                                                                                      >  2:228258/1‑62 (MQ=255)
                                                                                                                gcTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCa                                                                                                                      <  1:228258/62‑1 (MQ=255)
                                                                                                                      gATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCgctg                         >  5:231162/1‑149 (MQ=255)
                                                                                                                       aTGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCgctgc                        <  3:324528/149‑1 (MQ=255)
                                                                                                                          ggCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGa                     <  7:172241/149‑1 (MQ=255)
                                                                                                                           gCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGAt                    <  1:10185/149‑1 (MQ=255)
                                                                                                                            ccACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCgctgct                       <  8:206757/145‑1 (MQ=255)
                                                                                                                            ccACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCgctgct                       >  7:206757/1‑145 (MQ=255)
                                                                                                                            ccACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATc                   >  7:157659/1‑149 (MQ=255)
                                                                                                                              aCGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTg                 >  8:350385/1‑149 (MQ=255)
                                                                                                                              aCGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTg                 >  2:308161/1‑149 (MQ=255)
                                                                                                                              aCGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGGCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTg                 <  3:70465/149‑1 (MQ=255)
                                                                                                                                   aggatgcGGGCATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGGTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGa            <  6:78340/145‑1 (MQ=255)
                                                                                                                                    tgctgcGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGAt           >  1:335725/1‑149 (MQ=255)
                                                                                                                                         cGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGGGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGGCGTATGCTGCGCGCTGCTGATCTGGTTGATgcg        <  1:308161/147‑1 (MQ=255)
                                                                                                                                           tGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGATGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTAGATGCGCTgg    >  6:35476/1‑149 (MQ=255)
                                                                                                                                             aaTTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGt                                      <  5:325082/113‑1 (MQ=255)
                                                                                                                                             aaTTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGc  <  6:101106/149‑1 (MQ=255)
                                                                                                                                             aaTTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTAAACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGt                                      >  6:325082/1‑113 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
AAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCTACACCGACATGCCGATGCTGGTGATGCTGGAAGAACGCGACGGTTACTACGCTGCAGGTCGTATGCTGCGCGCTGCTGATCTGGTTGATGCGCTGGGC  >  NC_000913/1280672‑1280961

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑