Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 72,897 G→A 6.0% R790H (CGC→CAC)  U712_00330 → Stage II sporulation protein E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00688172,8970GA 6.0% 213.5 / 3.0 67R790H (CGC→CAC) U712_00330Stage II sporulation protein E
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base A (2/2);  total (34/33)
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.95e-01

TGATGAGTGATGGCATTTTTGAGGGGCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATCCATTCCTGTTCCGGCAATCTTTCAAAACAAACAAGAAATTTCATAACGCTTCCGTAT  >  CP006881/72765‑73032
                                                                                                                                    |                                                                                                                                       
tgatgaGTGATGGCATTTTTGAGGGGCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGcac                                                                                                                                      >  8:28442/1‑136 (MQ=255)
tgatgaGTGATGGCATTTTTGAGGGGCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGcac                                                                                                                                      >  3:27387/1‑136 (MQ=255)
tgatgaGTGATGGCATTTTTGAGGGGCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGcac                                                                                                                                      >  3:303293/1‑136 (MQ=255)
    gaGTGATGGCATTTTTGAGGGGCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATc                                                                                                                                <  4:415725/138‑1 (MQ=255)
            gCATTTTTGAGGGGCCTAAGCATGTGGAAAATCCCGACTTATGGATGAGGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAAt                                                                                                                       >  2:73329/1‑139 (MQ=255)
              aTTTTTGAGGGGCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTg                                                                                                                     <  5:71484/139‑1 (MQ=255)
                 tttGCGGGGCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAg                                                                                                                   <  4:223601/138‑1 (MQ=255)
                 tttGCGGGGCCTAAGCATGTGGAAAAACACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAg                                                                                                                   <  5:383150/138‑1 (MQ=255)
                 tttGAGGGGCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACATATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAg                                                                                                                   <  4:331598/138‑1 (MQ=255)
                 tttGAGGGGCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAg                                                                                                                   <  1:93480/138‑1 (MQ=255)
                 tttGAGGGGCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAg                                                                                                                   <  5:441797/138‑1 (MQ=255)
                 tttGAGGGGCATAAGCATGATGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAg                                                                                                                   <  5:114096/138‑1 (MQ=255)
                   tgcggggCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGa                                                                                                                 <  3:68265/135‑1 (MQ=255)
                   tGAGGGGCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGa                                                                                                                 <  5:348370/138‑1 (MQ=255)
                         gCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACTatat                                                                                                           >  1:17094/1‑138 (MQ=255)
                                  tgtgGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACtgttgtt                                                                                                 <  4:134048/139‑1 (MQ=255)
                                  tgtgGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACtgttgtt                                                                                                 <  5:295838/139‑1 (MQ=255)
                                    tgGAAAATCACGACTTATGGATGTAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCACACAAGATCCGGCCAAATTGAGGACGATATGACtgttgttgt                                                                                               >  7:368517/1‑139 (MQ=255)
                                    tgGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAg                                                                                                                   <  5:43964/119‑1 (MQ=255)
                                    tgGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAg                                                                                                                   >  6:43964/1‑119 (MQ=255)
                                    tgGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGatat                                                                                                           >  3:429923/1‑127 (MQ=255)
                                    tgGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGatat                                                                                                           <  4:429923/127‑1 (MQ=255)
                                    tgGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACtgttgttgt                                                                                               >  6:211387/1‑139 (MQ=255)
                                    tgGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACtgttgttgt                                                                                               >  6:249600/1‑139 (MQ=255)
                                    tgGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACtgttgttgt                                                                                               >  8:56952/1‑139 (MQ=255)
                                    tgGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACtgttgtggt                                                                                               >  4:265199/1‑139 (MQ=255)
                                    tgGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCAGCCGCACAAGATCCGGCCAAATTGAGGACGATATGACtgttgttgt                                                                                               >  7:144636/1‑139 (MQ=255)
                                    tgGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGAGTGCCGATTTGCTTATGGAAGAAGTCATCCACACAAGATCCGGCCAAATTGAGGa                                                                                                                 >  3:227553/1‑121 (MQ=255)
                                    tgGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGACTGCCGATTTGCTTATGGAAGAAGTCATCCACACAAGATCCGGCCAAATTGAGGa                                                                                                                 <  4:227553/121‑1 (MQ=255)
                                            cACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGa                                                                                       <  1:187104/139‑1 (MQ=255)
                                             aCGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCTGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGa                                                                                       <  3:425259/138‑1 (MQ=255)
                                                aCTTATGGAAGCAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGa                                                                                   <  6:252277/139‑1 (MQ=255)
                                                 cTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGAt                                                                                  <  3:398180/139‑1 (MQ=255)
                                                 cTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGAt                                                                                  <  7:28442/139‑1 (MQ=255)
                                                 cTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGAt                                                                                  <  7:297272/139‑1 (MQ=255)
                                                 cTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGAt                                                                                  <  2:291688/139‑1 (MQ=255)
                                                 cTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGAt                                                                                  <  5:211387/139‑1 (MQ=255)
                                                 cTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGAt                                                                                  <  4:303293/139‑1 (MQ=255)
                                                           taGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCACACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATaca                                                                         <  8:368517/137‑1 (MQ=255)
                                                           aaGCGTAAAATGAAAGGGTTGAAAACAAATGCCCCGCAGGAGATTGCCGAGTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAATTTGAGGACGATATTCCTGGTGGTGTCGTCCGGAGTGATCCCAATacac                                                                        >  8:315847/1‑139 (MQ=255)
                                                           aaGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATacac                                                                        >  6:213858/1‑139 (MQ=255)
                                                           aaGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTGGCTTATGGAAGAAGTCATCCGCACAAGATCCGGGCAACTTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGCTGACAATacac                                                                        >  8:321733/1‑139 (MQ=255)
                                                            aGCGTAAAATGCAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATcaca                                                                              <  1:211835/132‑1 (MQ=255)
                                                            aGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATcaca                                                                              >  2:211835/1‑132 (MQ=255)
                                                            aGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATca                                                                                <  8:201346/130‑1 (MQ=255)
                                                            aGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATca                                                                                >  7:201346/1‑130 (MQ=255)
                                                                  aaaTGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTg                                                                 <  7:54811/139‑1 (MQ=255)
                                                                  aaaTGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTg                                                                 >  7:323168/1‑139 (MQ=255)
                                                                  aaaTGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATAGAGGACGATATGAATGTTGTTGTCGTCCGGATTGATCACAATACACAGAAgag                                                                 >  8:294683/1‑137 (MQ=255)
                                                                   aaTGAACGGGGTGAACACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACACGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGGCCGGATTGATCACAATACACCGAAGTgg                                                                <  6:258488/139‑1 (MQ=255)
                                                                   aaTGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTgg                                                                <  5:249600/139‑1 (MQ=255)
                                                                        aaGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAg                                                                                                                   >  4:155561/1‑83 (MQ=255)
                                                                        aaGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAg                                                                                                                   <  3:155561/83‑1 (MQ=255)
                                                                        aaGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGCTCACAATACACCGCAGTGGGCGTc                                                           >  1:205111/1‑139 (MQ=255)
                                                                        aaGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATc                                                           >  5:134177/1‑139 (MQ=255)
                                                                        aaGGGGTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATc                                                           <  7:437601/139‑1 (MQ=255)
                                                                                       aTGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATCCATTCCTGTTCCgg                                             >  5:140175/1‑138 (MQ=255)
                                                                                                       tGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATCCATTCCTGTTCCGGCAATCTTTCAaaacaaa                            >  4:126521/1‑139 (MQ=255)
                                                                                                         ccGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATCCATTCCTGTTCCGGCAATCTTTCAaaacaaac                           >  5:408142/1‑138 (MQ=255)
                                                                                                                       gaagaaGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATCCATTCCTGTTCCGGCAATCTTTCAAAACAAACAAGAAATTTCATac             <  6:408142/138‑2 (MQ=255)
                                                                                                                           aagTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATCCATTCCTGTTCCGGCAATCTTTCAAAACAAACAAGAAATTTCATAACGCtt        >  8:420097/1‑139 (MQ=255)
                                                                                                                           aagTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATCCATTCCTGTTCCGGCAATCTTTCAAAACAAACAAGAAATTTCATAACGCtt        >  4:8641/1‑139 (MQ=255)
                                                                                                                                aTCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGTGCATCCATTCCTGTTCCGGCAATCTTTCAAAACAAACAAGAAATTTCATAACGCTTCCGta   >  4:200175/1‑139 (MQ=255)
                                                                                                                                aTCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATCCATTCCTGTTCCGGCAATCTTTCAAAACAAACAAGAAATTTCATAACGCTTCCGta   >  2:205572/1‑139 (MQ=255)
                                                                                                                                aTCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATCCATTCCTGTTCCGGCAATCTTTCAAAACAAACAAGAAATTTCATAACGCTTCCGta   >  6:323203/1‑139 (MQ=255)
                                                                                                                                aTCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATCCATTCCTGTTCCGGCAATCTTTCAAAACAAACAAGAAATTTCATAACGCGTCCGta   >  1:342097/1‑139 (MQ=255)
                                                                                                                                  ccGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATCCATTCCTGTTCCGGCAATCTTTCAAAACAAACAAGAAATTTCATAACGCTTCCGtat  >  4:2704/1‑138 (MQ=255)
                                                                                                                                    |                                                                                                                                       
TGATGAGTGATGGCATTTTTGAGGGGCCTAAGCATGTGGAAAATCACGACTTATGGATGAAGCGTAAAATGAAAGGGTTGAAAACAAATGACCCGCAGGAGATTGCCGATTTGCTTATGGAAGAAGTCATCCGCACAAGATCCGGCCAAATTGAGGACGATATGACTGTTGTTGTCGTCCGGATTGATCACAATACACCGAAGTGGGCATCCATTCCTGTTCCGGCAATCTTTCAAAACAAACAAGAAATTTCATAACGCTTCCGTAT  >  CP006881/72765‑73032

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: