Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,849,778 A→G 6.2% S345G (AGT→GGT)  hemN → oxygen‑independent coproporphyrinogen III dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,849,7780AG 6.2% 185.5 / 2.1 64S345G (AGT→GGT) hemNoxygen‑independent coproporphyrinogen III dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (27/33);  new base G (2/2);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

CCCTGCCCGACGACGAACTGGCCATTGCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGCTGGATTTCGAACC  >  NC_002947/4849629‑4849917
                                                                                                                                                     |                                                                                                                                           
cccTGCCCGACGACGAACTGGCCATTGCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAAcagca                                                                                                                                             <  2:299760/150‑1 (MQ=255)
  cTGCCCGACGACGAACTGGCCATTGCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGt                                                                                                                                           <  1:776348/150‑1 (MQ=255)
     cccGACGACGAACTGGCCATTGCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGAGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGAc                                                                                                                                        <  2:95739/150‑1 (MQ=255)
       cgacgacgaACTGGCCATTGCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCt                                                                                                                                      <  1:1065296/150‑1 (MQ=255)
         acgacgaACTGGCCATTGCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGCCCAGTACCGGCAGAAAGGCACAGACCTCa                                                                                                                                    >  2:603496/1‑150 (MQ=255)
                 cTGGCCATTGCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCGGTGACCTCAACACCTAc                                                                                                                            >  2:1330287/1‑150 (MQ=255)
                   ggCCATTGCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCa                                                                                                                          <  1:564281/150‑1 (MQ=255)
                   ggCCATTGCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCAACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCa                                                                                                                          <  2:158855/150‑1 (MQ=255)
                    gCCATTGCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTAcccc                                                                                                                         <  1:860233/150‑3 (MQ=255)
                         tGCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAg                                                                                                                         <  2:1249084/145‑1 (MQ=255)
                         tGCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAg                                                                                                                         >  1:1249084/1‑145 (MQ=255)
                          gCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAAcagca                                                                                                                                             >  1:1396823/1‑124 (MQ=255)
                          gCCCAGGAAGAAGGCACCCTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTCCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAAcagca                                                                                                                                             <  2:1396823/124‑1 (MQ=255)
                            ccAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGt                                                                                                                                           <  1:1676327/124‑1 (MQ=255)
                            ccAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGt                                                                                                                                           >  2:1676327/1‑124 (MQ=255)
                              aGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGt                                                                                                               <  1:240627/150‑1 (MQ=255)
                               ggAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTc                                                                                                              <  1:421349/150‑1 (MQ=255)
                                gaagaaGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTcc                                                                                                             <  1:93522/150‑1 (MQ=255)
                                   gaagGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAAc                                                                                                          >  1:1423636/1‑150 (MQ=255)
                                                         tCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCCGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGcgcg                                                                                    <  2:350586/150‑1 (MQ=255)
                                                             gggCTACACAACACACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACTTGTACTGCCAGAACAGCAGTGACCTAAACACCTACAAGGACAGCCTGTCGAACGCCCAGCTGGCCACATAGcgcgttgt                                                                                >  1:477477/1‑146 (MQ=255)
                                                               gCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGc                                                                              <  2:1345356/150‑1 (MQ=255)
                                                                  acacCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCtgt                                                                           <  2:341261/150‑1 (MQ=255)
                                                                      cacccacGGGCACTGCGACCTGATTGGCCTGGGTGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCaa                                                                       >  1:1112573/1‑150 (MQ=255)
                                                                      cacccacGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCaa                                                                       >  2:495171/1‑150 (MQ=255)
                                                                      cacccacGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCaa                                                                       >  2:1924948/1‑150 (MQ=255)
                                                                      cacccacGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGAGCGGCTTGCTGTGCaa                                                                       >  2:1349374/1‑150 (MQ=255)
                                                                        cccacGGGCACTGCGACCTGATTGGCCTGGTCTTGTCGGCGATCAGCCAGGTCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAAcc                                                                     >  1:44053/1‑150 (MQ=255)
                                                                         ccacGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAacac                                                                                                                                <  2:1535602/90‑1 (MQ=255)
                                                                         ccacGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAacac                                                                                                                                >  1:1535602/1‑90 (MQ=255)
                                                                          cacGGGCACTGCGCCCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGAGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCCGGACAGCTTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCac                                                                   <  2:1849691/150‑1 (MQ=255)
                                                                            cGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCacga                                                                 <  2:1930843/150‑1 (MQ=255)
                                                                               gCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCAcgacga                                                              >  2:373716/1‑150 (MQ=255)
                                                                                 aCTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGAcc                                                            >  2:179098/1‑150 (MQ=255)
                                                                                   tGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGc                                                          >  2:1858408/1‑150 (MQ=255)
                                                                                           gATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTg                                                                                                 <  1:743821/103‑1 (MQ=255)
                                                                                           gATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTg                                                                                                 >  2:743821/1‑103 (MQ=255)
                                                                                                      gCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCaaaggcca                                       <  2:1423636/150‑8 (MQ=255)
                                                                                                           tCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAgcgc                                                                                     <  2:1618560/99‑1 (MQ=255)
                                                                                                           tCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAgcgc                                                                                     >  1:1618560/1‑99 (MQ=255)
                                                                                                            cGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAAc                                 <  1:1924948/150‑1 (MQ=255)
                                                                                                             ggCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACt                                >  1:833029/1‑150 (MQ=255)
                                                                                                                      ccAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCa                                                                                                                          <  1:190143/51‑1 (MQ=255)
                                                                                                                      ccAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCa                                                                                                                          >  2:190143/1‑51 (MQ=255)
                                                                                                                        agagCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCAtt                     <  2:348727/150‑1 (MQ=255)
                                                                                                                        aGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGt                                                                                                               >  1:665238/1‑60 (MQ=255)
                                                                                                                        aGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGt                                                                                                               <  2:665238/60‑1 (MQ=255)
                                                                                                                         gATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATATGCCAttt                    >  2:982839/1‑150 (MQ=255)
                                                                                                                              gCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGc               >  1:842683/1‑150 (MQ=255)
                                                                                                                              gCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGc               <  1:208293/150‑1 (MQ=255)
                                                                                                                               cGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGCt              <  1:373716/150‑1 (MQ=255)
                                                                                                                                gACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGCTg             <  2:44053/150‑1 (MQ=255)
                                                                                                                                  ccTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGCTGGa           >  2:1073208/1‑150 (MQ=255)
                                                                                                                                   cTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGCTGGAt          >  1:1522852/1‑150 (MQ=255)
                                                                                                                                      tACTGCCAGAACAGCGGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACgacgac                                                             >  1:981701/1‑96 (MQ=255)
                                                                                                                                      tACTGCCAGAACAGCGGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACgacgac                                                             <  2:981701/96‑1 (MQ=255)
                                                                                                                                         tGCCAGAACAGCAGTGACCTCAACACCTCCCAGGACAGCATGTCCAACGCCAAGCTGGTCCCACCGCACCGCTTGCCGTGCAACCACGACTTCCTCATACTCCGGGCATTGATCCACCAACTGATCTGCCATTTCGAGCTGGATTTCGaa    <  2:220219/150‑1 (MQ=255)
                                                                                                                                         tGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGa                 >  2:1029393/1‑137 (MQ=255)
                                                                                                                                         tGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGa                 <  1:1029393/137‑1 (MQ=255)
                                                                                                                                         tGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGCTGGATTTCGaa    <  1:1657833/150‑1 (MQ=255)
                                                                                                                                          gCCAGAACAGCGGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGCTGGATTTCGAAc   <  1:1330287/150‑1 (MQ=255)
                                                                                                                                          gCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGAGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGCTGGATTTCGAAc   <  1:1349374/150‑1 (MQ=255)
                                                                                                                                           ccAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGt                                            <  1:347862/108‑1 (MQ=255)
                                                                                                                                           ccAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGt                                            >  2:347862/1‑108 (MQ=255)
                                                                                                                                           ccAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTAGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGCTGGATTTCGAAcc  >  2:269392/1‑150 (MQ=255)
                                                                                                                                              gAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGCTGGAtt         <  1:814093/140‑1 (MQ=255)
                                                                                                                                              gAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGCTGGAtt         >  2:814093/1‑140 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
CCCTGCCCGACGACGAACTGGCCATTGCCCAGGAAGAAGGCACCTTGCAGCGCAACTTCCAGGGCTACACCACCCACGGGCACTGCGACCTGATTGGCCTGGGCGTGTCGGCGATCAGCCAGATCGGCGACCTGTACTGCCAGAACAGCAGTGACCTCAACACCTACCAGGACAGCCTGTCCAACGCCCAGCTGGCCACCCAGCGCGGCTTGCTGTGCAACCACGACGACCGCATACGCCGGGCAGTGATCCAGCAACTGATCTGCCATTTCGAGCTGGATTTCGAACC  >  NC_002947/4849629‑4849917

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: