Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,706,643 T→A 6.4% L202Q (CTG→CAG)  PP_5009 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,706,6430TA 6.4% 190.8 / 4.0 62L202Q (CTG→CAG) PP_5009hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base A (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

AAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTCGACC  >  NC_002947/5706505‑5706792
                                                                                                                                          |                                                                                                                                                     
aaGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAgg                                                                                                                                            <  2:718650/150‑1 (MQ=255)
   aCAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGa                                                                                                                                         >  1:876466/1‑150 (MQ=255)
       ccAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACa                                                                                                                                     >  2:1382539/1‑150 (MQ=255)
       ccAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACa                                                                                                                                     >  1:844248/1‑150 (MQ=255)
            tctcGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCAGGCTGTCgtg                                                                                                                                             >  1:1210061/1‑135 (MQ=255)
            tctcGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCAGGCTGTCgtg                                                                                                                                             >  1:1211536/1‑135 (MQ=255)
            tctcGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCAGGCTGTCgtg                                                                                                                                             <  2:1211536/137‑3 (MQ=255)
            tctcGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCAGGCTGTCgtg                                                                                                                                             <  2:1210061/137‑3 (MQ=255)
                gcgcATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGa                                                                                                                                         <  2:651591/137‑1 (MQ=255)
                gcgcATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGa                                                                                                                                         >  1:651591/1‑137 (MQ=255)
                         ctggctggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGa                                                                                                                                              >  2:837760/1‑123 (MQ=255)
                         ctggctggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGa                                                                                                                                              <  1:837760/123‑1 (MQ=255)
                               ggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCgg                                                                                                             <  2:284530/150‑1 (MQ=255)
                                    acacCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGc                                                                                                                                           >  1:1303298/1‑115 (MQ=255)
                                    acacCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGc                                                                                                                                           <  2:1303298/115‑1 (MQ=255)
                                      acCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGa                                                                                                                                         <  1:632948/115‑1 (MQ=255)
                                      acCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGa                                                                                                                                         >  2:632948/1‑115 (MQ=255)
                                       ccATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAATCTTTGCAGATTATGCAAAAGCTGCAATGTACGGGCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGcacaca                                                                                                     <  1:804569/150‑1 (MQ=255)
                                            aaGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAg                                                                                                                                             <  1:166815/105‑1 (MQ=255)
                                            aaGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAg                                                                                                                                             >  2:166815/1‑105 (MQ=255)
                                            aaGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGc                                                                                                                                           >  1:876249/1‑107 (MQ=255)
                                            aaGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGc                                                                                                                                           <  2:876249/107‑1 (MQ=255)
                                            aaGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACATACCGtt                                                                                                <  1:518342/150‑1 (MQ=255)
                                                cTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTAccc                                                                                            >  2:111753/1‑150 (MQ=255)
                                                 tGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCt                                                                                           >  2:533726/1‑150 (MQ=255)
                                                   cATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTAc                                                                                         <  1:628355/150‑1 (MQ=255)
                                                        cTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGAt                                                                                    >  2:185137/1‑150 (MQ=255)
                                                         tCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGCACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGGTc                                                                                   >  1:352845/1‑150 (MQ=255)
                                                          cGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCt                                                                                  <  1:618016/150‑1 (MQ=255)
                                                           gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTa                                                                                 >  1:1198262/1‑150 (MQ=255)
                                                            agaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTAt                                                                                <  1:326315/150‑1 (MQ=255)
                                                             gaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTAtt                                                                               >  2:111539/1‑150 (MQ=255)
                                                                cctccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCt                                                                            <  2:844248/150‑1 (MQ=255)
                                                                 ctccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACATGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtc                                                                           >  2:102166/1‑150 (MQ=255)
                                                                 ctccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCggg                                                                                                            >  2:197602/1‑117 (MQ=255)
                                                                 ctccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCggg                                                                                                            <  1:197602/117‑1 (MQ=255)
                                                                     ctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtcatca                                                                       <  1:1382539/150‑1 (MQ=255)
                                                                        gCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTg                                                                              <  1:416998/140‑1 (MQ=255)
                                                                        gCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTg                                                                              >  2:416998/1‑140 (MQ=255)
                                                                         gccaGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATcaccac                                                                   >  2:174651/4‑150 (MQ=255)
                                                                                cTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgc                                                            >  1:1382036/1‑150 (MQ=255)
                                                                                cTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgc                                                            <  1:533726/150‑1 (MQ=255)
                                                                                      gatgatGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaaccgc                                                      >  2:471256/1‑150 (MQ=255)
                                                                                      gatgatGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaaccgc                                                      >  1:62168/1‑150 (MQ=255)
                                                                                        tgatgCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGTTGCAGCACACGCCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaaccgcaa                                                    <  1:174651/150‑1 (MQ=255)
                                                                                            gCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACaa                                                                                                                                    <  1:56851/66‑1 (MQ=255)
                                                                                            gCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACaa                                                                                                                                    >  2:56851/1‑66 (MQ=255)
                                                                                            gCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCacacac                                                                                                    >  2:640981/1‑98 (MQ=255)
                                                                                            gCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCacacac                                                                                                    <  1:640981/98‑1 (MQ=255)
                                                                                             cAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTgg                                               <  1:111753/150‑1 (MQ=255)
                                                                                              aaaaGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGt                                              >  2:1093412/1‑150 (MQ=255)
                                                                                                aaGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTc                                            <  1:1093412/150‑1 (MQ=255)
                                                                                                                ccTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTg                            >  2:353370/1‑150 (MQ=255)
                                                                                                                    cTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACAACATCACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGAGTCTGGAATGCAAGCCCTGcggt                        >  1:317907/1‑150 (MQ=255)
                                                                                                                        aCGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGc                                                                             <  1:722047/93‑1 (MQ=255)
                                                                                                                        aCGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGc                                                                             >  2:722047/1‑93 (MQ=255)
                                                                                                                               acacGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGCTTCGGTTTCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtcg             >  1:1170135/1‑150 (MQ=255)
                                                                                                                                  cGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcggttcggt                   >  2:224011/1‑141 (MQ=255)
                                                                                                                                  cGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcggttcggt                   <  1:224011/141‑1 (MQ=255)
                                                                                                                                  cGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCa          <  1:111539/150‑1 (MQ=255)
                                                                                                                                   gACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtc              >  1:1352486/1‑145 (MQ=255)
                                                                                                                                   gACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtc              <  2:1352486/145‑1 (MQ=255)
                                                                                                                                          tGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTCGAcc  >  2:290943/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
AAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTCGACC  >  NC_002947/5706505‑5706792

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