Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,999,490 T→C intergenic (‑345/‑134) ygeV ← / → ygeW sigma‑54‑dependent Fis family transcriptional regulator/knotted carbamoyltransferase YgeW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,999,4900TC97.1% 234.7 / ‑1.1 68intergenic (‑345/‑134)ygeV/ygeWsigma‑54‑dependent Fis family transcriptional regulator/knotted carbamoyltransferase YgeW
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (29/37);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

CCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGTCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAGTGAAATAAACCACTTTGCCTGTCATTCCACTACCGGGACTTTATGATGAA  >  NZ_CP009273/2999355‑2999628
                                                                                                                                       |                                                                                                                                          
ccccTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCtt                                                                                                                                         >  2:885395/1‑139 (MQ=255)
      tATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAg                                                                                                                                   >  2:187350/1‑139 (MQ=255)
            aTTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCcgcgc                                                                                                                             >  2:810660/1‑139 (MQ=255)
            aTTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCcgc                                                                                                                               <  1:932228/137‑1 (MQ=255)
             tttCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCcgcgcg                                                                                                                            <  2:242821/139‑1 (MQ=255)
             tttCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCcgcgcg                                                                                                                            >  2:896665/1‑139 (MQ=255)
             tttCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCcgcgcg                                                                                                                            >  2:504285/1‑139 (MQ=255)
                  gATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATaa                                                                                                                       >  2:1021144/1‑139 (MQ=255)
                   aTACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGAt                                                                                                                         <  2:1050789/136‑1 (MQ=255)
                   aTACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGAt                                                                                                                         <  1:261629/136‑1 (MQ=255)
                   aTACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATaa                                                                                                                       <  2:738926/138‑1 (MQ=255)
                   aTACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATaa                                                                                                                       <  2:266120/138‑1 (MQ=255)
                   aTACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATaa                                                                                                                       <  2:343082/138‑1 (MQ=255)
                         agcagcAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGTCTTTTGCAGCCGCGCGGATAATTCTgg                                                                                                                 >  2:379654/1‑138 (MQ=255)
                         agcagcAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTgg                                                                                                                 >  2:1002649/1‑138 (MQ=255)
                         agcagcAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTgg                                                                                                                 >  1:930126/1‑138 (MQ=255)
                          gcagcaAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGca                                                                                                               <  2:337068/139‑1 (MQ=255)
                          gcagcaAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGca                                                                                                               <  2:84275/139‑1 (MQ=255)
                          gcagcaAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGca                                                                                                               <  2:1253595/139‑1 (MQ=255)
                          gcagcaAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGca                                                                                                               <  1:1103050/139‑1 (MQ=255)
                              caAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACt                                                                                                           >  2:327010/1‑139 (MQ=255)
                                aataaataACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTa                                                                                                         <  1:596558/139‑1 (MQ=255)
                                         cgcgATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTc                                                                                                <  2:114778/139‑1 (MQ=255)
                                         cgcgATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTc                                                                                                <  1:327010/139‑1 (MQ=255)
                                         cgcgATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTc                                                                                                <  2:824549/139‑1 (MQ=255)
                                         cgcgATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTc                                                                                                <  1:187350/139‑1 (MQ=255)
                                           cgATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACACCCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTccc                                                                                              <  1:524014/139‑1 (MQ=255)
                                             aTAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAg                                                                                            >  2:121716/1‑139 (MQ=255)
                                             aTAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAg                                                                                            >  2:991660/1‑139 (MQ=255)
                                             aTAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAg                                                                                            >  1:383518/1‑139 (MQ=255)
                                              tAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGc                                                                                                                >  1:575363/1‑118 (MQ=255)
                                              tAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGc                                                                                                                <  2:575363/118‑1 (MQ=255)
                                                    aaGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGc                                                                                     <  1:810660/139‑1 (MQ=255)
                                                    aaGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGc                                                                                     <  2:668453/139‑1 (MQ=255)
                                                                   gTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAAt                                                                      <  1:744624/139‑1 (MQ=255)
                                                                        ttACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACt                                                                                                           <  1:765700/97‑1 (MQ=255)
                                                                        ttACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACt                                                                                                           >  2:765700/1‑97 (MQ=255)
                                                                         tACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTa                                                                                                   >  1:741119/1‑104 (MQ=255)
                                                                         tACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTa                                                                                                   <  2:741119/104‑1 (MQ=255)
                                                                                  ccAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTca                                                       <  1:504285/139‑1 (MQ=255)
                                                                                    aGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGTCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTcaca                                                     <  1:379654/139‑1 (MQ=255)
                                                                                    aGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTcaca                                                     <  1:247516/139‑1 (MQ=255)
                                                                                     gAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAg                                                    <  1:1002649/139‑1 (MQ=255)
                                                                                     gAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAg                                                    <  1:121716/139‑1 (MQ=255)
                                                                                        gCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTc                                                                                                <  1:460859/92‑1 (MQ=255)
                                                                                        gCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTc                                                                                                >  2:460859/1‑92 (MQ=255)
                                                                                             cATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGaaa                                                                              >  2:327868/1‑105 (MQ=255)
                                                                                             cATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGaaa                                                                              <  1:327868/105‑1 (MQ=255)
                                                                                                     tgtgATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAGTGAAATAAACCACtt                                     <  1:448294/138‑1 (MQ=255)
                                                                                                      gtgATCAACCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAGTGAAATAAACCACTTTg                                   >  2:397505/1‑139 (MQ=255)
                                                                                                             aCCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGca                                                                                                               <  2:829227/56‑1 (MQ=255)
                                                                                                             aCCCCACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGca                                                                                                               >  1:829227/1‑56 (MQ=255)
                                                                                                              ccccACACATTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAGTGAAATAAACCACTTTGCCTGTCAt                           >  1:642298/1‑139 (MQ=255)
                                                                                                                 cacacaTTATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAGTGAAATAAACCACTTTGCCTGTCATTcc                        <  1:1162033/139‑1 (MQ=255)
                                                                                                                        tATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTc                                                                                                >  2:428307/1‑60 (MQ=255)
                                                                                                                        tATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTc                                                                                                <  1:428307/60‑1 (MQ=255)
                                                                                                                        tATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGc                                                                                     >  1:696112/1‑71 (MQ=255)
                                                                                                                        tATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGc                                                                                     <  2:696112/71‑1 (MQ=255)
                                                                                                                        tATCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAGTGAAATAAACCACTTTGCCTGTCATTCCACTACCg                 >  1:1078785/1‑139 (MQ=255)
                                                                                                                         aTCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAGTGAAATAAACCACTTTGCCTGTCATTCCACTACCgg                >  2:548721/1‑139 (MQ=255)
                                                                                                                         aTCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAGTGAAATAAACCACTTTGCCTGTCATTCCACTACCgg                >  1:75475/1‑139 (MQ=255)
                                                                                                                         aTCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAGTGAAATAAACCACTTTGCCTGTCATTCCACTACCgg                >  2:99666/1‑139 (MQ=255)
                                                                                                                         aTCCGTCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAGTGAAATAAACCACTTTGCCTGTCATTCCACTACCg                 <  1:991660/138‑1 (MQ=255)
                                                                                                                              tCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTgg                                                                                                                 >  2:1135255/1‑37 (MQ=255)
                                                                                                                              tCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTgg                                                                                                                 <  1:1135255/37‑1 (MQ=255)
                                                                                                                              tCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGt                                                                                                     >  2:340256/1‑49 (MQ=255)
                                                                                                                              tCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGt                                                                                                     <  1:340256/49‑1 (MQ=255)
                                                                                                                              tCAAATTAGCCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAGTGAAATAAACCACTTTGCCTGTCATTCCACTACCGGGACtt           >  1:11297/1‑139 (MQ=255)
                                                                                                                                       tCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAAGCACTTTGGCAAGTCACAGTTAAATAAACCTCTTTGCCTGTCATTCCACTGCCGGGACTTTATGATGaa  <  2:1179439/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                          
CCCCTTTATCTTATTTCTGATACTTAGCAGCAAATAAATAACGCGATAAAAAAAGCCAAACGTTTTCGTATTTTACAAACAACCAGAAGCTGGCATCAATTTGTGATCAACCCCACACATTATCCGTCAAATTAGTCTTTTGCAGCCGCGCGGATAATTCTGGCACACTTATTGTTAGTCCCAGGTATAGCTGTGAAAACACCAATCACTTTGGCAAGTCACAGTGAAATAAACCACTTTGCCTGTCATTCCACTACCGGGACTTTATGATGAA  >  NZ_CP009273/2999355‑2999628

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