Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,244,986 C→A 6.6% G12V (GGC→GTC)  aruH ← arginine‑‑pyruvate transaminase AruH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,244,9860CA 6.6% 230.6 / 2.0 72G12V (GGC→GTC) aruHarginine‑‑pyruvate transaminase AruH
Reads supporting (aligned to +/- strand):  ref base C (35/32);  new base A (3/2);  total (38/34)
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

CCAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCGGCCT  >  NC_002947/4244839‑4245133
                                                                                                                                                   |                                                                                                                                                    
ccAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGcc                                                                                                                                                     >  8:226855/1‑149 (MQ=255)
ccAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTATCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGcc                                                                                                                                                     >  1:270595/1‑149 (MQ=255)
      cgcAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGAt                                                                                                                                               <  7:131976/149‑1 (MQ=255)
       gcAGGCTGTCGTTGGCCGCCTCGACGTTAGGCGCGGGAGTGTCGAAGGCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCCCAGCCTTCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCATCGCCGGCGACGGCGCTg                                                                                                                                              <  7:54108/149‑1 (MQ=255)
           gCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGtt                                                                                                                                          >  2:58604/1‑149 (MQ=255)
           gCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGGCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGtt                                                                                                                                          <  8:261072/149‑1 (MQ=255)
                        gccTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGc                                                                                                                            <  2:270595/149‑1 (MQ=255)
                          cTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCg                                                                                                                           >  7:102267/1‑148 (MQ=255)
                           tCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCg                                                                                                                           >  1:62961/1‑147 (MQ=255)
                              acgaTAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAgcg                                                                                                                      <  5:203198/149‑1 (MQ=255)
                              acgaTAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAgcg                                                                                                                      >  3:54469/1‑149 (MQ=255)
                                     gcgcGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGt                                                                                                                                      >  3:63194/1‑127 (MQ=255)
                                     gcgcGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGt                                                                                                                                      <  4:63194/127‑1 (MQ=255)
                                         gggAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAAtct                                                                                                           >  2:143847/1‑149 (MQ=255)
                                            aGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCt                                                                                                        <  1:143847/149‑1 (MQ=255)
                                            aGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCt                                                                                                        <  6:161598/149‑1 (MQ=255)
                                             gtgtCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCtt                                                                                                       >  8:339740/1‑149 (MQ=255)
                                              tgtCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCtt                                                                                                       <  7:226855/148‑1 (MQ=255)
                                                  gAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAAt                                                                                                  <  8:102267/149‑1 (MQ=255)
                                                      tCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGc                                                                                               >  4:257822/1‑148 (MQ=255)
                                                          ggTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGc                                                                                                                            <  1:327658/115‑1 (MQ=255)
                                                          ggTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGc                                                                                                                            >  2:327658/1‑115 (MQ=255)
                                                           gTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAAttt                                                                                          >  6:116562/1‑148 (MQ=255)
                                                            tCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCa                                                                                        <  4:54469/149‑1 (MQ=255)
                                                             cccccaccGACACCAGCAGGATGTCCTTGTCCCCAGCCTGCATCGCCACGGCGCGGTATTGGATGTCCCACGCCGCATCGCCGTCGCCGGCGAAGCGTTGGGGGATCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAg                                                                                       <  2:62961/147‑1 (MQ=255)
                                                                  accGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcgac                                                                                  >  7:90595/1‑149 (MQ=255)
                                                                    cGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcgacga                                                                                <  5:158234/149‑1 (MQ=255)
                                                                          cagcagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCt                                                                          >  4:116980/1‑149 (MQ=255)
                                                                          cagcagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCt                                                                          >  2:10498/1‑149 (MQ=255)
                                                                           agcagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCtt                                                                         >  4:146917/1‑149 (MQ=255)
                                                                            gcagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTc                                                                        <  1:58604/149‑1 (MQ=255)
                                                                              agGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCgg                                                                      <  2:108616/149‑1 (MQ=255)
                                                                                gATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGt                                                                    >  4:63392/1‑149 (MQ=255)
                                                                                   gtccttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGt                                                                              >  2:201288/1‑136 (MQ=255)
                                                                                   gtccttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGt                                                                              <  1:201288/136‑1 (MQ=255)
                                                                                     ccttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGaa                                                               >  7:211476/1‑149 (MQ=255)
                                                                                      cttgtcctCAGCCTGCAGCGCCAGGGCGCGCTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAAt                                                              >  1:1393/1‑149 (MQ=255)
                                                                                       ttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGaa                                                               <  6:15568/147‑1 (MQ=255)
                                                                                       ttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGaa                                                               >  5:15568/1‑147 (MQ=255)
                                                                                       ttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAAtt                                                             >  4:300427/1‑149 (MQ=255)
                                                                                       ttgtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAAtt                                                             >  2:362801/1‑149 (MQ=255)
                                                                                         gtcctCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGc                                                           >  1:154799/1‑149 (MQ=255)
                                                                                             tCAGCCTGCAGCGCCAGGGCGAGGTAGTGGATGTACCAGGCCGAAGCGCCGTCGCCGGCGATGCGTTGGGAGAGCTTGGCGTAGCGCATCTCAATCTCCTTTCAATGGCAATTTCAGTagaa                                                                                  >  4:144035/1‑121 (MQ=255)
                                                                                               aGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATa                                                     <  3:63392/149‑1 (MQ=255)
                                                                                                 ccTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACg                                                   >  7:296717/1‑149 (MQ=255)
                                                                                                  cTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTc                                                                        >  7:218299/1‑127 (MQ=255)
                                                                                                  cTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTc                                                                        <  8:218299/127‑1 (MQ=255)
                                                                                                     cAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGgtgt                                               <  7:339740/149‑1 (MQ=255)
                                                                                                        cgcCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTAc                                             >  5:37532/1‑148 (MQ=255)
                                                                                                          ccAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTAGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACaaa                                          >  1:340872/1‑149 (MQ=255)
                                                                                                               gcgcGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGc                                     <  8:296717/149‑1 (MQ=255)
                                                                                                               gcgcGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGc                                     <  3:257822/149‑1 (MQ=255)
                                                                                                                 gcgGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGa                                   <  1:10498/149‑1 (MQ=255)
                                                                                                                          gATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGt                                                                                      >  2:312239/1‑89 (MQ=255)
                                                                                                                          gATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGt                                                                                      <  1:312239/89‑1 (MQ=255)
                                                                                                                          gATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACgcgc                          >  4:161576/1‑149 (MQ=255)
                                                                                                                                ccAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGTGGTTAGGCGTAGCGCATGGCAATCTCCTTTCAATGGCaa                                                                                             <  1:55575/77‑1 (MQ=255)
                                                                                                                                ccAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCaa                                                                                             >  2:55575/1‑76 (MQ=255)
                                                                                                                                ccAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTc                                                                        >  2:212516/1‑97 (MQ=255)
                                                                                                                                ccAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTc                                                                        <  1:212516/97‑1 (MQ=255)
                                                                                                                                ccAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGGATTGCCCCTTACGTGGTACAAATCGGCGAGAACGGCGCCGTTGa                    >  3:38984/1‑149 (MQ=255)
                                                                                                                                ccAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCAGGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTg                     <  4:328456/148‑1 (MQ=255)
                                                                                                                                  ccaacGCATCGCCGGCGACGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGAtt                  >  6:204041/5‑149 (MQ=255)
                                                                                                                                  aGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGAtt                  >  7:75922/1‑149 (MQ=255)
                                                                                                                                    aacGCATCGCCGGCGACGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAAc                              <  3:38355/133‑1 (MQ=255)
                                                                                                                                    aacGCATCGCCGGCGACGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTAGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAAc                              >  4:38355/3‑135 (MQ=255)
                                                                                                                                     ccGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATttct               <  8:75922/149‑1 (MQ=255)
                                                                                                                                     ccGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATttct               <  2:154799/149‑1 (MQ=255)
                                                                                                                                      cGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATttctt              <  5:116562/149‑1 (MQ=255)
                                                                                                                                       gCAGAGCGGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGAAGTGTATCTTCGGGTTGGTATTGCCCCATACGGTGTAGAAATCGGCGAGCAACGCGCCGTTGATttcttc             <  2:340872/149‑1 (MQ=255)
                                                                                                                                          tcgcCGGCGACGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCa           >  5:243001/2‑148 (MQ=255)
                                                                                                                                           cgcCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCAtc         <  3:300427/149‑1 (MQ=255)
                                                                                                                                              cGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCg      >  1:146434/1‑149 (MQ=255)
                                                                                                                                                  gCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCACCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCGGCCt  >  5:65612/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCAGCGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCGGCCT  >  NC_002947/4244839‑4245133

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: