Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,207,303 A→G 6.0% N79D (AAC→GAC)  PP_4588 → nitroreductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,207,3030AG 6.0% 192.3 / 2.8 66N79D (AAC→GAC) PP_4588nitroreductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (34/28);  new base G (2/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

CGAAGAAACCCTGCTGAGCTTCCTCGAAGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTGGGGGCACCTGGCCCTGCAGGCCAGCATCA  >  NC_002947/5207161‑5207450
                                                                                                                                              |                                                                                                                                                   
cGAAGAAACCCTGCTGAGCTTCCTCGAAGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTggg                                                                                                                                               >  5:150603/1‑149 (MQ=255)
 gaagaaACCCTGCTGAGCTTCCTCGAAGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGgc                                                                                                                                              >  6:166943/1‑149 (MQ=255)
 gaagaaACCCTGCTGAGCTTCCTCGAAGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACATGAGCCTGATGTTGCCCTTCAATCGCAACTGGgc                                                                                                                                              >  2:301775/1‑149 (MQ=255)
         ccTGCTGAGCTTCCTCGAAGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAg                                                                                                                                       <  5:4422/148‑1 (MQ=255)
         ccTGCTGAGCTTCCTCGAAGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAg                                                                                                                                       >  6:4422/1‑148 (MQ=255)
                    tCCTCGAAGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCgctg                                                                                                                           >  6:312365/1‑149 (MQ=255)
                     ccTCGAAGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCgctgg                                                                                                                          >  3:38348/1‑149 (MQ=255)
                       tCGAAGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCgctggtg                                                                                                                        >  2:219722/1‑149 (MQ=255)
                       tCGAAGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTACTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCgctggtg                                                                                                                        >  1:299687/1‑149 (MQ=255)
                          aaGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCgctggtgctg                                                                                                                     <  7:135618/149‑1 (MQ=255)
                          aaGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTACTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCgctggtgctg                                                                                                                     <  2:299687/149‑1 (MQ=255)
                          aaGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCAGCGACACGCCGAACTGGGAGCGTTTCCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCgctggtgctg                                                                                                                     >  7:271838/1‑149 (MQ=255)
                              cgcccgcTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCgctg                                                                                                                           <  8:321065/139‑1 (MQ=255)
                              cgcccgcTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCgctg                                                                                                                           >  7:321065/1‑139 (MQ=255)
                                       ggCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCaa                                                                                                          <  5:312365/147‑1 (MQ=255)
                                               cGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGtt                                                                                                 <  1:301775/148‑1 (MQ=255)
                                                 gCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTcgc                                                                                              >  1:27641/1‑149 (MQ=255)
                                                   cTATAACTCGCAGCCTTGGCGCTTCCTTTATGCTCGCCGCGACACTCCGAACTTGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATAGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCAtt                                                                                                               <  8:287060/130‑1 (MQ=255)
                                                   cTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCgg                                                                                            >  5:19104/1‑149 (MQ=255)
                                                    tataACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCgg                                                                                            >  4:176124/1‑148 (MQ=255)
                                                      taACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCgg                                                                                            <  3:176124/146‑1 (MQ=255)
                                                         cTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCgg                                                                                      <  1:317352/149‑1 (MQ=255)
                                                          tCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCggg                                                                                     >  4:119694/1‑149 (MQ=255)
                                                            gCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGgcg                                                                                   >  5:69724/1‑149 (MQ=255)
                                                               gCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTa                                                                                <  6:19104/149‑1 (MQ=255)
                                                                 cTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCgctggtgc                                                                                                                       <  3:119745/108‑1 (MQ=255)
                                                                 cTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCgctggtgc                                                                                                                       >  4:119745/1‑108 (MQ=255)
                                                                  ttGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACg                                                                              >  8:389792/1‑148 (MQ=255)
                                                                      cgcTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGCTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGgaaga                                                                         <  6:413043/149‑1 (MQ=255)
                                                                       gcTTCCTTTATCCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCgctggtgctg                                                                                                                     >  7:366576/1‑104 (MQ=255)
                                                                          tCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTc                                                                                                                                  >  8:250781/1‑88 (MQ=255)
                                                                           ccTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGaaaaa                                                                     <  1:213817/148‑1 (MQ=255)
                                                                              ttATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCgg                                                                 >  3:250452/1‑149 (MQ=255)
                                                                                   cccGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCtg                                                            >  2:153389/1‑149 (MQ=255)
                                                                                    ccgccgCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGACGCTACGGAAGAAAAACCGGCGCtgt                                                           <  1:219722/149‑1 (MQ=255)
                                                                                     cgccgcGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCtgtg                                                          <  6:149643/149‑1 (MQ=255)
                                                                                     cgccgcGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCtgtg                                                          >  5:149643/1‑149 (MQ=255)
                                                                                         gcgACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGcac                                                      <  5:305060/149‑1 (MQ=255)
                                                                                         gcgACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGca                                                       >  1:8912/1‑148 (MQ=255)
                                                                                              acGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCtt                                                 >  3:102259/1‑149 (MQ=255)
                                                                                               cGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTc                                                <  7:257751/149‑1 (MQ=255)
                                                                                               cGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTc                                                <  6:150603/149‑1 (MQ=255)
                                                                                                gCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGATGAAAAACCGGCGCTGTGGCACACCTTc                                                >  2:100271/1‑148 (MQ=255)
                                                                                                gCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTc                                                >  4:45987/1‑148 (MQ=255)
                                                                                                       tGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGTGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCtt                                                 >  4:160073/1‑140 (MQ=255)
                                                                                                       tGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGTGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCtt                                                 <  3:160073/140‑1 (MQ=255)
                                                                                                         ggAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACggg                                        <  4:250452/147‑1 (MQ=255)
                                                                                                           aGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGATGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGctct                                    <  1:100271/149‑1 (MQ=255)
                                                                                                                       gCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCGAAGACCACGTTCGCGGCCCCGGTCGCTACGGAAGAAAAACCGGCGCTGTGGCTCACCTTCGACACGGTCTCTGCCTGGGGGCAc                        <  8:271838/149‑1 (MQ=255)
                                                                                                                        ccTGCTGGTGCCCTTCAATCGCGACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCg                                                                  <  1:231475/106‑1 (MQ=255)
                                                                                                                        ccTGCTGGTGCCCTTCAATCGCGACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCg                                                                  <  1:225071/106‑1 (MQ=255)
                                                                                                                        ccTGCTGGTGCCCTTCAATCGCGACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCg                                                                  >  2:231475/1‑106 (MQ=255)
                                                                                                                        ccTGCTGGTGCCCTTCAATCGCGACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCg                                                                  >  2:225071/1‑106 (MQ=255)
                                                                                                                        ccTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTGGGGGCAcc                       <  7:389792/149‑1 (MQ=255)
                                                                                                                                     ttCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTgggg                            >  5:267578/1‑131 (MQ=255)
                                                                                                                                     ttCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTgggg                            <  6:267578/131‑1 (MQ=255)
                                                                                                                                     ttCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTGGGGGCACCTGGCCGTGCAGGc          >  7:169790/1‑149 (MQ=255)
                                                                                                                                     ttCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTGGGGGCACCTGGCCCTGCAGGc          >  7:315763/1‑149 (MQ=255)
                                                                                                                                     ttCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTGGGGGCACCTGGCCCTGCAGGc          >  5:398362/1‑149 (MQ=255)
                                                                                                                                        aaTCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTGGGGGCACCTGG‑CCTGCAgggccag      <  3:119694/149‑6 (MQ=255)
                                                                                                                                           cGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTGGGGGCACCTGGCCCTGCAGGCCAGCAt    >  3:369552/1‑149 (MQ=255)
                                                                                                                                             cAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTggg                             >  7:160541/1‑122 (MQ=255)
                                                                                                                                             cAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTggg                             <  8:160541/122‑1 (MQ=255)
                                                                                                                                             cAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTGGGGGCACCTGGcc                  >  2:134396/1‑133 (MQ=255)
                                                                                                                                             cAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTGGGGGCACCTGGcc                  <  1:134396/133‑1 (MQ=255)
                                                                                                                                             cAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTGGGGGCACCTGGCCCTGCAGGCCAGCATCa  <  6:69724/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CGAAGAAACCCTGCTGAGCTTCCTCGAAGCCGCCCGCTGGGCACCTTCGGCCTATAACTCGCAGCCTTGGCGCTTCCTTTATGCCCGCCGCGACACGCCGAACTGGGAGCGTTACCTGAGCCTGCTGGTGCCCTTCAATCGCAACTGGGCGCAACAGGCGTCGGCGCTGGTGCTGGTCATTTCCAAGACCACGTTCGCGGCCCCGGGCGCTACGGAAGAAAAACCGGCGCTGTGGCACACCTTCGACACGGGCTCTGCCTGGGGGCACCTGGCCCTGCAGGCCAGCATCA  >  NC_002947/5207161‑5207450

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: ‑