Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,195,985 G→A 8.8% S177L (TCG→TTG) ‡ ybgJ ← allophanate hydrolase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,195,9850GA 8.8% 194.1 / 6.0 68S177L (TCG→TTG) ‡ybgJallophanate hydrolase subunit 1
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base A (4/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TCCACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAA  >  NC_002947/5195838‑5196133
                                                                                                                                                    |                                                                                                                                                     
tCCACATCAGCACCGTTGGCCAGGTATTCGCTGCGCGATTCCGGCAGGTAGCGCACCCTGTCAGCGGGTTGC‑TAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGa                                                                                                                                                      <  2:225338/149‑1 (MQ=255)
         ccaccCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCgg                                                                                                                                             <  1:247323/149‑1 (MQ=255)
              cATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAgg                                                                                                                                        <  1:139613/149‑1 (MQ=255)
               aTGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGgc                                                                                                                                       >  5:87642/1‑149 (MQ=255)
                   ccAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGgccgcc                                                                                                                                   <  4:19127/149‑1 (MQ=255)
                        aaTTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTg                                                                                                                              <  6:87642/149‑1 (MQ=255)
                          ttCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGCGAAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTg                                                                                                                              >  1:274808/1‑148 (MQ=255)
                          ttCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTg                                                                                                                              <  2:274808/147‑1 (MQ=255)
                          ttCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGgc                                                                                                                            >  6:352181/1‑149 (MQ=255)
                               tGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCg                                                                                                                       >  2:122076/1‑149 (MQ=255)
                                 cgcgAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                      <  3:255220/148‑1 (MQ=255)
                                 cgcgAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                      <  5:324095/148‑1 (MQ=255)
                                 cgcgAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCg                                                                                                                       <  3:9191/147‑1 (MQ=255)
                                   cgAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCaa                                                                                                                   <  5:340696/149‑1 (MQ=255)
                                     aaaCCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATa                                                                                                                 <  5:306710/149‑1 (MQ=255)
                                         cGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGccgccg                                                                                                                                  >  4:281186/1‑128 (MQ=255)
                                         cGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGccgccg                                                                                                                                  <  3:281186/128‑1 (MQ=255)
                                          ggCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACTCTCACGGTCGAACAACCGAACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGTCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCAc                                                                                                            <  2:363560/149‑1 (MQ=255)
                                           gCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACg                                                                                                           >  7:238735/1‑149 (MQ=255)
                                              cGAAGCGCAGCCGGTCACCGGGTTGC‑AAGAGGGTGTATCGCTGACGCTCACGGTCGACCAACCGCCCCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTg                                                                                                        <  4:307681/149‑1 (MQ=255)
                                               gAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCACGGCAATACCCACGCTGc                                                                                                       >  1:281837/1‑149 (MQ=255)
                                                 ttgcaaCCCGGTGACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGccc                                                                                                     >  2:75241/6‑149 (MQ=255)
                                                       cccGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCcgccgc                                                                                                >  3:101085/1‑148 (MQ=255)
                                                           gTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCAcgc                                                                                           >  8:230145/1‑149 (MQ=255)
                                                              aCCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCTGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCt                                                                                        >  1:198360/1‑149 (MQ=255)
                                                              aCCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                      <  3:139173/119‑1 (MQ=255)
                                                              aCCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCgg                                                                                                                      >  4:139173/1‑119 (MQ=255)
                                                                       c‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCg                                                                                                                                              >  8:4043/1‑86 (MQ=255)
                                                                       c‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCg                                                                                                                                              <  7:4043/86‑1 (MQ=255)
                                                                            agGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACt                                                                           <  8:238735/149‑1 (MQ=255)
                                                                             ggCTGTAACCCTCCCGCTCACGGGCGAACCACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTc                                                                          <  2:312548/149‑1 (MQ=255)
                                                                                tGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACt                                                                           <  3:144417/145‑1 (MQ=255)
                                                                                tGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACt                                                                           >  4:144417/1‑145 (MQ=255)
                                                                                 gTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGAc                                                                      >  6:138965/1‑149 (MQ=255)
                                                                                  tAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGAc                                                                      >  6:35197/1‑148 (MQ=255)
                                                                                      ccTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAg‑gg                                                                >  4:296018/1‑149 (MQ=255)
                                                                                          actctcacgGTCGAACAACCGCAGCGGGGTGCGCCCGATCTGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTg                                                             <  2:281837/145‑1 (MQ=255)
                                                                                          acgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTg                                                             <  5:138965/148‑1 (MQ=255)
                                                                                           cgctcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGGGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGc                                                           <  1:75241/149‑1 (MQ=255)
                                                                                              tcacgGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAg                                                        >  7:271906/1‑149 (MQ=255)
                                                                                                      gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATa                                                                                                                 <  2:312005/85‑1 (MQ=255)
                                                                                                      gAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAGCGCCTGGCCAgccctcgt                                                 <  3:160767/149‑8 (MQ=255)
                                                                                                        acaacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCtcgtc                                               >  7:187365/1‑148 (MQ=255)
                                                                                                        acaacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCtcgtc                                               >  7:52850/1‑148 (MQ=255)
                                                                                                          aacCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGAc                                            >  3:165688/1‑149 (MQ=255)
                                                                                                            ccGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGAc                                            >  7:327532/1‑147 (MQ=255)
                                                                                                                 ccGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGccc                                     <  8:271906/149‑1 (MQ=255)
                                                                                                                  cGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCa                                    >  2:340149/1‑149 (MQ=255)
                                                                                                                           cccGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCg                            >  5:136944/1‑148 (MQ=255)
                                                                                                                              gATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAAc                        >  7:136928/1‑149 (MQ=255)
                                                                                                                                   ggTTCCAGCCACCCGGCGATACTGCCGGGTAGGCTGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                   <  2:198360/149‑1 (MQ=255)
                                                                                                                                   ggTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                   <  6:136944/149‑1 (MQ=255)
                                                                                                                                   ggTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGt                   <  8:327532/149‑1 (MQ=255)
                                                                                                                                     ttCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCtcgtc                                               >  3:115359/1‑119 (MQ=255)
                                                                                                                                     ttCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCtcgtc                                               <  4:115359/119‑1 (MQ=255)
                                                                                                                                     ttCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGc                 >  1:298075/1‑149 (MQ=255)
                                                                                                                                           ctaCCCGGCAGTATCGCCGGGTGGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGc                              <  2:75270/128‑1 (MQ=255)
                                                                                                                                           ctaCCCGGCAGTATCGCCGGGTGGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGc                              >  1:75270/3‑130 (MQ=255)
                                                                                                                                           ctaCCCGGCAGTATCGCCGGGTGGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGcc           >  8:65633/3‑149 (MQ=255)
                                                                                                                                           ctaCCCGGCAGTATCGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAgg                                                       >  1:273432/3‑105 (MQ=255)
                                                                                                                                           ctaCCCGGCAGTATCGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAgg                                                       <  2:273432/103‑1 (MQ=255)
                                                                                                                                           ccACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGTGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGc            >  6:368393/1‑148 (MQ=255)
                                                                                                                                           ccACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGTGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGc            <  5:368393/148‑1 (MQ=255)
                                                                                                                                             aCCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGcc           >  7:145124/1‑147 (MQ=255)
                                                                                                                                             aCCCGGCAGTATCGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGcc           >  8:310013/1‑147 (MQ=255)
                                                                                                                                               ccGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAgg       <  3:296018/149‑1 (MQ=255)
                                                                                                                                                  gCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAgg       <  5:323915/146‑1 (MQ=255)
                                                                                                                                                  gCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAgg       >  6:323915/1‑146 (MQ=255)
                                                                                                                                                    gATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCaaa  >  6:368010/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCCACATCACCACCCATGGCCAGGAATTCGCTGCGCGAAACCGGCACGAAGCGCACCCGGTCACCGGGTTGC‑AAGAGGCTGTAACCCTCACGCTCACGGTCGAACAACCGCACCGGGGTGCGCCCGATCAGGTTCCAGCCACCCGGCGATACTGCCGGGTAGGCCGCCGTCTGGCGCTCGGCAATACCCACGCTGCCCGCCGCCACGCGCTTGCGTGGGGTACTCAGACGAG‑GGCTGGCCAGGCGCTCGTCGACCAGGCCCATGAAGCCGAAACCCGGTGCGAAGCCCAGGGCAAA  >  NC_002947/5195838‑5196133

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