Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 447,159:1 +C 6.6% intergenic (+39/‑263) PP_0368 → / → PP_0369 acyl‑CoA dehydrogenase/diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947447,1591.C 6.6% 246.6 / 11.1 61intergenic (+39/‑263)PP_0368/PP_0369acyl‑CoA dehydrogenase/diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base . (32/25);  new base C (2/2);  total (34/27)
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.80e-01

CTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTCCTA  >  NC_002947/447015‑447302
                                                                                                                                                 |                                                                                                                                               
cTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGc                                                                                                                                             >  7:92506/1‑149 (MQ=255)
    cAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCAAAAACCCGCCTCCCGGTCAGGGCAGGCGGGt                                                                                                                                         <  3:40/149‑1 (MQ=255)
         cATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGttttt                                                                                                                                     >  7:42324/1‑148 (MQ=255)
         cATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGttttt                                                                                                                                     >  6:79156/1‑148 (MQ=255)
             cTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCgg                                                                                                                                           <  4:54231/138‑1 (MQ=255)
             cTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCgg                                                                                                                                           >  3:54231/1‑138 (MQ=255)
             cTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACg                                                                                                                                 <  8:42324/148‑1 (MQ=255)
              tGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTc                                                                                                                               >  4:159480/1‑149 (MQ=255)
                       cccgtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGttttt                                                                                                                                     <  2:235846/134‑1 (MQ=255)
                       cccgtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGttttt                                                                                                                                     >  1:235846/1‑134 (MQ=255)
                        ccgtacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCa                                                                                                                     <  5:79156/149‑1 (MQ=255)
                           tacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGttttt                                                                                                                                     <  6:53853/130‑1 (MQ=255)
                           tacccATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGttttt                                                                                                                                     >  5:53853/1‑130 (MQ=255)
                                       gACCATGCTGTGGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCg                                                                                                        <  5:222754/147‑1 (MQ=255)
                                       gACCATGCTGTAGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCg                                                                                                        >  6:222754/1‑147 (MQ=255)
                                        aCCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAg                                                                                                     <  5:322088/149‑1 (MQ=255)
                                        aCCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAg                                                                                                     <  1:209332/149‑1 (MQ=255)
                                          cATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGa                                                                                                   <  8:92506/149‑1 (MQ=255)
                                             gCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCgg                                                                                                                                           >  7:292808/1‑106 (MQ=255)
                                             gCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCgg                                                                                                                                           <  8:292808/106‑1 (MQ=255)
                                              cTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCgg                                                                                               <  5:33506/149‑1 (MQ=255)
                                                 tCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCggg                                                                                              <  2:322497/147‑1 (MQ=255)
                                                   gggCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTa                                                                                          >  3:226173/1‑149 (MQ=255)
                                                    ggCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGt                                                                                                              >  6:188945/1‑128 (MQ=255)
                                                    ggCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGt                                                                                                              <  5:188945/128‑1 (MQ=255)
                                                         caacaaCCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCcgc                                                                                     >  1:271837/1‑148 (MQ=255)
                                                                        gATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCctg                                                                     >  3:222586/1‑149 (MQ=255)
                                                                           ggACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAaga                                                                                >  4:64006/1‑135 (MQ=255)
                                                                           ggACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAaga                                                                                <  3:64006/135‑1 (MQ=255)
                                                                            gACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTgg                                                                 <  4:222586/149‑1 (MQ=255)
                                                                               gAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCa                                                              <  4:226173/149‑1 (MQ=255)
                                                                                        ttCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCAt                                                      >  3:156544/1‑148 (MQ=255)
                                                                                                    tACTGATACTGGCGCTGTACCCATAACCCGCTTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCAcc                                          <  6:250592/148‑1 (MQ=255)
                                                                                                     aCTGATACTGGCGCTGTACCCAAAACCCGCCAGCGCTGACCGAGAGGCGGGTTTTATACGTCTTCAGGTTCAGCGTGGACCCTCGGAGGAGCGGGGTTAACCGCGAAGAGTCCGGCCCTCCTGGACAATCTCAATAATCCCTTCACCg                                         >  5:250592/1‑148 (MQ=255)
                                                                                                       tGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggccg                                                                           <  6:91159/112‑1 (MQ=255)
                                                                                                       tGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAggccg                                                                           >  5:91159/1‑112 (MQ=255)
                                                                                                         aTACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGTCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTga                                     >  3:240109/1‑147 (MQ=255)
                                                                                                                 gcTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCCTCAGCCGTTCACCGTTGCCATCGCTAc                            >  4:272876/1‑149 (MQ=255)
                                                                                                                       cccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCAcc                                          >  1:22348/1‑129 (MQ=255)
                                                                                                                       cccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCAcc                                          <  2:22348/129‑1 (MQ=255)
                                                                                                                       cccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTAtt                      <  4:156544/149‑1 (MQ=255)
                                                                                                                       cccAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTAt                       >  1:340234/1‑148 (MQ=255)
                                                                                                                         cAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTcg                    >  2:80796/1‑149 (MQ=255)
                                                                                                                          aaaaCCCGCCTCCCGGTCAGGGCAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCAACAGCCGTTCACCGTTGCCAGCGCTACCGTATTcgc                   >  2:23538/1‑149 (MQ=25)
                                                                                                                           aaaCCCGCCTCCCGGTCAGGGCAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCa                                   <  5:286985/132‑1 (MQ=25)
                                                                                                                           aaaCCCGCCTCCCGGTCAGGGCAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCa                                   >  6:286985/1‑132 (MQ=25)
                                                                                                                              cccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCggg                >  3:96305/1‑148 (MQ=255)
                                                                                                                               ccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCctgc                                                                    >  8:241084/1‑95 (MQ=255)
                                                                                                                               ccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCctgc                                                                    <  7:241084/95‑1 (MQ=255)
                                                                                                                               ccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGt               >  1:268304/1‑148 (MQ=255)
                                                                                                                               ccGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGt               >  4:56154/1‑148 (MQ=255)
                                                                                                                                  ccTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAAc           <  4:96305/149‑1 (MQ=255)
                                                                                                                                   cTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAAcc          >  2:84937/1‑149 (MQ=255)
                                                                                                                                   cTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAAcc          <  7:32636/149‑1 (MQ=255)
                                                                                                                                      cccTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTg                                      >  1:269929/1‑118 (MQ=255)
                                                                                                                                      cccTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTg                                      <  6:101892/118‑1 (MQ=255)
                                                                                                                                      cccTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTg                                      <  2:269929/118‑1 (MQ=255)
                                                                                                                                      cccTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTg                                      >  5:101892/1‑118 (MQ=255)
                                                                                                                                        cTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTc     >  4:154670/1‑149 (MQ=255)
                                                                                                                                         tGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTcc    >  4:74244/1‑149 (MQ=255)
                                                                                                                                           acggggATGCGTGTTTTTTACGTCTTAAGGTTCAGCGTGGTCCCTCTTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCTTATTCGCGGGTAAACCCGCTCCTa  <  4:240109/146‑1 (MQ=21)
                                                                                                                                           aCCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTCCTa  >  5:94384/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CTTCCAGCGCATCCTGCCACGTACCCGTACCCATGTGGCGACCATGCTGTCGGGCGCCAACAACCTGATGGCGATGGACGAAGAGCACTTCGGCCTGTCGTACTGATACTGGCGCTGTACCCAAAACCCGCCTGCCCTGACCGGGAGGCGGGTTTTTTACGTCTTCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCGTATTCGCGGGTAAACCCGCTCCTA  >  NC_002947/447015‑447302

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: