Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 838,760 T→A 6.2% *314L (TAA→TTA)  prs ← ribose‑phosphate pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947838,7600TA 6.2% 199.0 / 4.7 64*314L (TAA→TTA) prsribose‑phosphate pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base A (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

GCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGGTGTTGGTCACCACCAGCTCGTCCA  >  NC_002947/838612‑838907
                                                                                                                                                    |                                                                                                                                                   
gcTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTa                                                                                                                                                    <  1:280607/150‑1 (MQ=255)
                  aaTCCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGAt                                                                                                                                  <  1:1128956/150‑1 (MQ=255)
                           aTTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCCTCGCGGAACATCGCGCTCATCGActcttc                                                                                                                         >  2:582861/1‑150 (MQ=255)
                                ggttCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACATCCGTTGGGGCGGGGCGCTTCACATCAACACGTCTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGc                                                                                                                    <  2:583764/148‑1 (MQ=255)
                                     tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGa                                                                                                                                   >  2:818716/1‑130 (MQ=255)
                                     tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGa                                                                                                                                   <  1:818716/130‑1 (MQ=255)
                                     tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATg                                                                                                               >  2:1645405/1‑150 (MQ=255)
                                                cccccGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGa                                                                                                                 >  1:221801/1‑137 (MQ=255)
                                                cccccGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACAGCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGa                                                                                                                 <  2:221801/137‑1 (MQ=255)
                                                 ccccGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCt                                                                                                   >  1:1497510/1‑150 (MQ=255)
                                                      aGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGc                                                                                              <  2:547253/150‑1 (MQ=255)
                                                         gCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGAc                                                                                           <  1:1153095/150‑1 (MQ=255)
                                                           gTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGCTGCAGCTGACCGCTTCAGCGACaa                                                                                         >  1:1318067/1‑150 (MQ=255)
                                                            tttGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCGGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCATCGACAAc                                                                                        >  2:739526/1‑150 (MQ=255)
                                                                           aaCGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAg                                                                         >  1:891527/1‑150 (MQ=255)
                                                                                 gATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGTCg                                                                   <  1:1645405/150‑1 (MQ=255)
                                                                                  aTGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCgg                                                                  >  1:240398/1‑150 (MQ=255)
                                                                                   tGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGa                                                                                <  1:1122538/135‑1 (MQ=255)
                                                                                   aggCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGa                                                                                >  2:1122538/2‑135 (MQ=255)
                                                                                       aaaaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGt                                                                                                                       >  2:772956/1‑92 (MQ=255)
                                                                                       aaaaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGt                                                                                                                       <  1:772956/92‑1 (MQ=255)
                                                                                        aaaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGc                                                                                  >  1:736196/1‑128 (MQ=255)
                                                                                        aaaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGc                                                                                  <  2:736196/128‑1 (MQ=255)
                                                                                        aaaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCCTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACgg                                                            <  1:582861/150‑1 (MQ=255)
                                                                                         aaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGaa                                                                                                                                              >  2:1083310/1‑67 (MQ=255)
                                                                                         aaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGaa                                                                                                                                              <  1:1083310/67‑1 (MQ=255)
                                                                                         aaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGc                                                                                  <  1:619135/127‑1 (MQ=255)
                                                                                         aaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGc                                                                                  >  2:619135/1‑127 (MQ=255)
                                                                                              ccccgccccAACCTGTGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGa                                                                 <  2:1445148/139‑1 (MQ=255)
                                                                                              ccccgccccAACCTGTGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGa                                                                 >  1:1445148/1‑139 (MQ=255)
                                                                                               cccgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAgg                                                     >  2:1554884/1‑150 (MQ=255)
                                                                                               cccgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAgg                                                     >  2:1356989/1‑150 (MQ=255)
                                                                                                ccgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTAAGCGGGACac                                                                                                                                           <  1:621981/63‑2 (MQ=255)
                                                                                                ccgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTAAGCGGGACac                                                                                                                                           >  2:621981/1‑62 (MQ=255)
                                                                                                     ccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGag                                               <  2:1308166/150‑1 (MQ=255)
                                                                                                      cAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCg                                                                   <  1:577578/129‑1 (MQ=255)
                                                                                                      cAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCg                                                                   >  2:577578/1‑129 (MQ=255)
                                                                                                      cAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGagc                                              >  2:72796/1‑150 (MQ=255)
                                                                                                          acaGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCg                                          >  1:20265/1‑150 (MQ=255)
                                                                                                             ggTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCa                                       <  2:240398/150‑1 (MQ=255)
                                                                                                              gTTGGGGCGGGGCGCTTCACATCAACACGTGTTCCGCTAAGCGGGACac                                                                                                                                           >  1:646528/1‑48 (MQ=21)
                                                                                                              gTTGGGGCGGGGCGCTTCACATCAACACGTGTTCCGCTAAGCGGGACac                                                                                                                                           <  2:646528/49‑2 (MQ=21)
                                                                                                              gTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGt                                                           >  2:1283776/1‑129 (MQ=255)
                                                                                                              gTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGt                                                           <  1:1283776/129‑1 (MQ=255)
                                                                                                               ttGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGCACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGa                                     <  1:589501/150‑1 (MQ=255)
                                                                                                                    gcggggcgCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCg                                <  2:20265/150‑1 (MQ=255)
                                                                                                                    gcggggcgCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCg                                >  2:1224502/1‑150 (MQ=255)
                                                                                                                             ttCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGt                                                           <  1:595396/114‑1 (MQ=255)
                                                                                                                             ttCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGt                                                           >  2:595396/1‑114 (MQ=255)
                                                                                                                             ttCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGgtgt                       <  1:1421402/150‑1 (MQ=255)
                                                                                                                                  aTCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCt                                                   >  2:1520887/1‑117 (MQ=255)
                                                                                                                                  aTCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCt                                                   <  1:1520887/117‑1 (MQ=255)
                                                                                                                                       cacGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAAc                            >  2:118966/1‑135 (MQ=255)
                                                                                                                                       cacGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAAc                            <  1:118966/135‑1 (MQ=255)
                                                                                                                                              cccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGGTGTTGGTCACCACCAGCTCg      >  2:1246477/1‑150 (MQ=255)
                                                                                                                                              cccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGGTGTTGGTCACCACCAGCTCg      >  2:434673/1‑150 (MQ=255)
                                                                                                                                               ccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGc                                                                        >  1:59874/1‑83 (MQ=255)
                                                                                                                                               ccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGc                                                                        <  2:59874/83‑1 (MQ=255)
                                                                                                                                               ccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGc                                 <  1:280277/122‑1 (MQ=255)
                                                                                                                                               ccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGc                                 >  2:280277/1‑122 (MQ=255)
                                                                                                                                               ccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGGTGTTGGTCACCACCAGCTCGt     <  1:1246477/150‑1 (MQ=255)
                                                                                                                                                 gCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATa                                                               <  2:1099806/90‑1 (MQ=255)
                                                                                                                                                 gCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATa                                                               >  1:1099806/1‑90 (MQ=255)
                                                                                                                                                  cTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGGTGTTGGTCACCACCAGCTCGTCCa  <  2:1318067/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGGTGTTGGTCACCACCAGCTCGTCCA  >  NC_002947/838612‑838907

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: