Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,846 T→A 6.1% intergenic (‑135/+45) aroQ‑I ← / ← dsbD‑I type II 3‑dehydroquinate dehydratase/thiol/disulfide interchange protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947650,8460TA 6.1% 249.8 / 4.8 82intergenic (‑135/+45)aroQ‑I/dsbD‑Itype II 3‑dehydroquinate dehydratase/thiol/disulfide interchange protein 1
Reads supporting (aligned to +/- strand):  ref base T (45/32);  new base A (3/2);  total (48/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 = 9.45e-01

CGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCGCA  >  NC_002947/650703‑650993
                                                                                                                                                |                                                                                                                                                   
cGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATc                                                                                                                                                >  4:71704/1‑149 (MQ=255)
    gCCATCTGCGGATTCCTTCGATTTGTAGGTCAAATCGACACAGCGCGGCGAGTGTGCCGGAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCgg                                                                                                                                             <  6:256058/148‑1 (MQ=255)
             ggATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCAtt                                                                                                                                    >  6:207222/1‑148 (MQ=255)
              gATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCAttt                                                                                                                                   >  8:79476/1‑148 (MQ=255)
                   cTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAt                                                                                                                             >  8:306119/1‑149 (MQ=255)
                    ttGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAtt                                                                                                                            >  4:249308/1‑149 (MQ=255)
                    ttGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAtt                                                                                                                            >  6:284312/1‑149 (MQ=255)
                         tttGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGTCAcata                                                                                                                                      >  8:35590/1‑131 (MQ=255)
                         tttGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATTTAAATATCCGGTCAcata                                                                                                                                      <  7:35590/134‑4 (MQ=255)
                             tAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGc                                                                                                                   >  3:170740/1‑149 (MQ=255)
                              aGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCt                                                                                                                  >  3:318568/1‑149 (MQ=255)
                                ggCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTcg                                                                                                                >  6:263050/1‑149 (MQ=255)
                                ggCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTcg                                                                                                                >  5:297561/1‑149 (MQ=255)
                                  cAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTcgcg                                                                                                              <  3:71704/149‑1 (MQ=255)
                                        gACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAg                                                                                                         <  6:1386/148‑1 (MQ=255)
                                           acaGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGaa                                                                                                          <  2:10976/144‑1 (MQ=255)
                                           acaGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGaa                                                                                                          >  1:10976/1‑144 (MQ=255)
                                           acaGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTta                                                                                                      >  4:134608/1‑148 (MQ=255)
                                             aGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGtt                                                                                                       <  3:193284/145‑1 (MQ=255)
                                             aGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGtt                                                                                                       >  4:193284/1‑145 (MQ=255)
                                              gcgcGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGa                                                                                                  <  5:284312/149‑1 (MQ=255)
                                              gcgcGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGa                                                                                                  <  5:263050/149‑1 (MQ=255)
                                               cgcgGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGAc                                                                                                 <  7:157976/149‑1 (MQ=255)
                                               cgcgGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGAc                                                                                                 <  4:42060/149‑1 (MQ=255)
                                                  ggcgAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGa                                                                                              >  7:260230/1‑149 (MQ=255)
                                                  ggcgAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGa                                                                                              <  8:260230/149‑1 (MQ=255)
                                                      agggtgCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGGGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACttt                                                                                          <  7:73680/146‑1 (MQ=255)
                                                            ccGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCcgcg                                                                                    >  5:147817/1‑149 (MQ=255)
                                                            ccGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCcgcg                                                                                    >  4:94684/1‑149 (MQ=255)
                                                            ccGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCcgcg                                                                                    >  1:298943/1‑149 (MQ=255)
                                                             cGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACAGACTTTGCcgcg                                                                                    >  4:48178/1‑148 (MQ=255)
                                                               tAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCAc                                                                                 <  7:79476/149‑1 (MQ=255)
                                                                 aaGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCAcgc                                                                               >  1:300261/1‑149 (MQ=255)
                                                                   gCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCAcgcgc                                                                             <  4:170740/149‑1 (MQ=255)
                                                                   gCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCAcgcgc                                                                             <  4:318568/149‑1 (MQ=255)
                                                                           aaGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAAt                                                                     >  7:29187/1‑149 (MQ=255)
                                                                             gtgtCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCt                                                                   >  4:80055/1‑149 (MQ=255)
                                                                                  cAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGAAGATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgcgc                                                               >  1:98193/1‑149 (MQ=255)
                                                                                  cAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGtt                                                                                                       >  4:267930/1‑108 (MQ=255)
                                                                                  cAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGtt                                                                                                       <  3:267930/108‑1 (MQ=255)
                                                                                  cAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTgcgc                                                               <  8:262016/148‑1 (MQ=255)
                                                                                   aGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGcgcg                                                              >  7:256265/1‑148 (MQ=255)
                                                                                   aGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGcgcg                                                              >  5:324530/1‑148 (MQ=255)
                                                                                   aGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGa                                                             <  2:300261/149‑1 (MQ=255)
                                                                                     ttCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAAt                                                           >  7:149014/1‑149 (MQ=255)
                                                                                        ccGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCa                                                        >  8:111675/1‑149 (MQ=255)
                                                                                          ggCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCAcc                                                      <  3:57980/149‑1 (MQ=255)
                                                                                           gCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCg                                                     <  1:96623/149‑1 (MQ=255)
                                                                                            cAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCgg                                                    <  3:249308/149‑1 (MQ=255)
                                                                                              aTAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCa                                                  <  5:207222/149‑1 (MQ=255)
                                                                                               tAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCAt                                                 <  7:306119/149‑1 (MQ=255)
                                                                                               tAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCAt                                                 <  6:297561/149‑1 (MQ=255)
                                                                                                 gCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCAtt                                                >  3:46916/1‑148 (MQ=255)
                                                                                                     gCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTAtct                                           >  6:333131/1‑149 (MQ=255)
                                                                                                       aCGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCg                                         >  6:308045/1‑149 (MQ=255)
                                                                                                        cGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCAttt                                               >  4:139231/1‑142 (MQ=255)
                                                                                                        cGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCAttt                                               <  3:139231/142‑1 (MQ=255)
                                                                                                        cGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGc                                        <  3:48178/149‑1 (MQ=255)
                                                                                                         gATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCg                                                     >  2:313425/1‑135 (MQ=255)
                                                                                                         gATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCg                                                     <  1:313425/135‑1 (MQ=255)
                                                                                                           tGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCg                                                     <  8:23475/133‑1 (MQ=255)
                                                                                                           tGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCg                                                     >  7:23475/1‑133 (MQ=255)
                                                                                                               cGAGGTACGCGCGAAGTATGTGACCGGATATTTAAATACGGTCACTCATTTCAGGATTTTTACGAGCTCGCGGGAAGTTATAGACAGACTTTGCCGCGCACGCGCGTCAGAATCTGCGAGAATTCCCCGGCAGTTATCTCGCCAATCAc                                 >  4:83951/1‑149 (MQ=255)
                                                                                                                gAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAACTCACCGGCATTTATCTCGCCAATCAcc                                <  2:298943/149‑1 (MQ=255)
                                                                                                                  ggTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCAccc                               >  4:149324/1‑148 (MQ=255)
                                                                                                                   gTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCAcc                                >  8:62045/1‑146 (MQ=255)
                                                                                                                   gTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCAcc                                <  7:62045/146‑1 (MQ=255)
                                                                                                                      cgcgcgAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATc                           <  8:256265/148‑1 (MQ=255)
                                                                                                                         gcgAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCg                       >  2:161439/1‑149 (MQ=255)
                                                                                                                           gAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGAATTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCgg                      >  1:240130/1‑148 (MQ=255)
                                                                                                                             aGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTaa                   <  3:134608/149‑1 (MQ=255)
                                                                                                                             aGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTaa                   <  6:147817/149‑1 (MQ=255)
                                                                                                                                 agtgACCGGATTTAAATATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTc               >  7:297693/2‑149 (MQ=255)
                                                                                                                                   tgACCGGATTTAAATATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGc             >  4:169369/1‑149 (MQ=255)
                                                                                                                                   tgACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGc             >  5:112799/1‑149 (MQ=255)
                                                                                                                                    gACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCgat            >  6:136016/1‑147 (MQ=255)
                                                                                                                                     aCCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCTGTAATTTCGCTg           >  4:181191/1‑149 (MQ=255)
                                                                                                                                     aCCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGc             <  7:259083/147‑1 (MQ=255)
                                                                                                                                     aCCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGc             >  8:259083/1‑147 (MQ=255)
                                                                                                                                        ggATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCg        >  3:242958/1‑149 (MQ=255)
                                                                                                                                         gATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGt       >  7:186003/1‑149 (MQ=255)
                                                                                                                                             atTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCg    >  8:278784/1‑149 (MQ=255)
                                                                                                                                              tttAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCg    >  6:161522/1‑148 (MQ=255)
                                                                                                                                                tAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCgca  <  8:29187/148‑1 (MQ=255)
                                                                                                                                                tAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCgca  <  3:80055/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CGTTGCCATCTGCGGATTCCTTGGATTTGTAGGGCAATTCGACACAGCGCGGCGAGTGTGCCGTAAAGCAACGGCAAGTGTCCAGTTCCCGGCAATAGCCAGCACGATGTCCGAGGTTCGCGCGAAGTATGTGACCGGATATTTAAATCCGGTCACTCATTTCAGGATTTTTCCGAGCTCGCGGGAAGTTATAGACCGACTTTGCCGCGCACGCGCGTCAGAATCTGCGCGAATTCACCGGCATTTATCTCGCCAATCACCCGATCAGCGGTAATTTCGCTGCCGTTCGCA  >  NC_002947/650703‑650993

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: