Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 87,493 T→A 15.0% N98I (AAT→ATT)  PP_0076 ← choline betaine‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294787,4930TA15.0% 157.2 / 20.2 60N98I (AAT→ATT) PP_0076choline betaine‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (28/23);  new base A (5/4);  total (33/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 = 8.95e-01

GATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAA  >  NC_002947/87348‑87638
                                                                                                                                                 |                                                                                                                                                 
gATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGc                                                                                                                                               >  2:1158757/1‑150 (MQ=255)
 aTTCCGTACAGTTTTTTTTCGAACTGTCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCg                                                                                                                                              >  2:964871/1‑150 (MQ=255)
 aTTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCg                                                                                                                                              <  2:432093/150‑1 (MQ=255)
                  gTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGaa                                                                                                                             <  2:833884/150‑1 (MQ=255)
                  gTCGAACTGCCCGGCGTGTTTGTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGCCACTGTGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGaa                                                                                                                             <  2:586192/150‑1 (MQ=255)
                         tGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCAGCGCATCGTCCGACACAGCTAGTGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTtgc                                                                                                                      <  1:142299/150‑1 (MQ=255)
                                 gtgtTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGc                                                                                                                                               <  2:202716/117‑1 (MQ=255)
                                 gtgtTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGc                                                                                                                                               >  1:202716/1‑117 (MQ=255)
                                              tCAGCGCAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCg                                                                                                 >  2:566023/1‑150 (MQ=255)
                                              tCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCg                                                                                                 <  1:66713/150‑1 (MQ=255)
                                                aGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAg                                                                                               >  1:167298/1‑150 (MQ=255)
                                                       tCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAg                                                                                           <  1:166887/147‑1 (MQ=255)
                                                       tCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAg                                                                                           >  2:166887/1‑147 (MQ=255)
                                                        ccTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGt                                                                                       <  1:225487/150‑1 (MQ=255)
                                                         cTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTc                                                                                      >  2:762340/1‑150 (MQ=255)
                                                           tCACCCCGGCATTCCACACATAGGCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGa                                                                                    <  2:460735/150‑1 (MQ=255)
                                                               cccGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTg                                                                                >  1:1102164/1‑150 (MQ=255)
                                                                        ccACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCtt                                                                       <  2:833199/150‑1 (MQ=255)
                                                                        ccACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCtt                                                                       >  1:775486/1‑150 (MQ=255)
                                                                          acacaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAt                                                                                                                                                   <  2:65393/72‑1 (MQ=255)
                                                                          acacaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAt                                                                                                                                                   >  1:65393/1‑72 (MQ=255)
                                                                           cacaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCa                                                                                                        >  1:989155/1‑114 (MQ=255)
                                                                           cacaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCa                                                                                                        <  2:989155/114‑1 (MQ=255)
                                                                            acaTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAgg                                                                   >  1:153885/1‑150 (MQ=255)
                                                                             caTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGa                                                                     <  1:1109870/147‑1 (MQ=255)
                                                                             caTAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGa                                                                     >  2:1109870/1‑147 (MQ=255)
                                                                                    gggCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTa                                                           <  1:1158757/150‑1 (MQ=255)
                                                                                        aCTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACtt                                                                                                                                      <  1:205744/71‑1 (MQ=255)
                                                                                        aCTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACtt                                                                                                                                      >  2:205744/1‑71 (MQ=255)
                                                                                         cTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAg                                                                                           >  1:389118/1‑113 (MQ=255)
                                                                                         cTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAg                                                                                           <  2:389118/113‑1 (MQ=255)
                                                                                          tGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAgatgat                                                     >  2:903834/1‑150 (MQ=255)
                                                                                           gCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATc                                                    <  1:762340/150‑1 (MQ=255)
                                                                                                    aaCTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGcc                                           >  2:761660/1‑150 (MQ=255)
                                                                                                      cTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAATGGCGATgt                                                                                                                                          >  1:427255/1‑52 (MQ=255)
                                                                                                      cTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCAATGGCGATgt                                                                                                                                          <  2:427255/53‑2 (MQ=255)
                                                                                                           ttCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGa                                                                                                >  1:259243/1‑90 (MQ=255)
                                                                                                           tccccTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGa                                                                                                <  2:259243/88‑1 (MQ=255)
                                                                                                               cTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCa                                                                                                            >  2:725886/1‑74 (MQ=255)
                                                                                                               cTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCa                                                                                                            <  1:725886/74‑1 (MQ=255)
                                                                                                               cTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGAt                                >  1:1119673/1‑150 (MQ=255)
                                                                                                               cTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGAt                                >  1:255197/1‑150 (MQ=255)
                                                                                                               cTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGAt                                >  2:346561/1‑150 (MQ=255)
                                                                                                                  cAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCtt                             >  2:727303/1‑150 (MQ=255)
                                                                                                                                 ctgctgCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAg              >  2:825591/1‑150 (MQ=255)
                                                                                                                                  tgctgcACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGTCTGGCGATCTTc                            >  2:498429/1‑135 (MQ=255)
                                                                                                                                  tgctgcACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTc                            <  1:498429/135‑1 (MQ=255)
                                                                                                                                  tgctgcACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAgg             <  1:771420/150‑1 (MQ=255)
                                                                                                                                      gttcaTCGCCAATGGCGATGTACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCttt         >  1:1098814/4‑150 (MQ=255)
                                                                                                                                      gcACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCttt         >  1:475567/1‑150 (MQ=255)
                                                                                                                                       ttcaTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCa               <  2:507655/141‑1 (MQ=255)
                                                                                                                                       ttcaTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCa               >  1:507655/3‑143 (MQ=255)
                                                                                                                                        tcaTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTcc       <  2:25586/149‑1 (MQ=255)
                                                                                                                                         caTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCa      <  2:475567/150‑1 (MQ=255)
                                                                                                                                         caTCGCAATTGGCGATGAACAGTGCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGCCGTCCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCa      <  1:566023/150‑1 (MQ=255)
                                                                                                                                          aTCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATAAAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAg     >  1:219754/1‑150 (MQ=255)
                                                                                                                                          aTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTc                            <  1:154809/127‑1 (MQ=255)
                                                                                                                                          aTCGCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTc                            >  2:154809/1‑127 (MQ=255)
                                                                                                                                             gCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCTGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTaa  >  1:535031/1‑150 (MQ=255)
                                                                                                                                             gCCAATGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTaa  >  2:238163/1‑150 (MQ=255)
                                                                                                                                               cATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAgatgat                                                     >  2:301543/1‑97 (MQ=255)
                                                                                                                                               cATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAgatgat                                                     <  1:301543/97‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GATTCCGTACAGTTTTTTGTCGAACTGCCCGGCGTGTTTTTGCAGGTCAGCGAAGTCCTTCACCCCGGCATTCCACACATAGTCGGGCACTGCGAGGGTGAACTCGGTTCCCTCCAGGTTGCGCGACAGCTGCTGCACATCGCCATTGGCGATGAACTTGTCATGGAAGCCTTGCTGCGCCGGCATCCAGTTTCCGAGGAAGGCGTCGACCTGGCCGTCCTTGAGGCCGCCGTAGATGATCGGTACGGCCAGGCTGTCGATCTTCACCTGGTAACCCAGGCTTTCCAGTAA  >  NC_002947/87348‑87638

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