Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC47.4% 72.2 / 36.0 57T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (29/1);  new base C (0/27);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 2.00e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGC  >  NC_000913/1376358‑1376653
                                                                                                                                                      |                                                                                                                                                    
ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAccc                                                                                                                                                   >  2:265887/1‑151 (MQ=255)
ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAccc                                                                                                                                                   >  1:305620/1‑151 (MQ=255)
   ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCg                                                                                                                                                <  2:13424/151‑1 (MQ=255)
      ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                             >  1:813982/1‑151 (MQ=255)
       cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                             >  2:613322/1‑150 (MQ=255)
               tCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                                    >  1:106229/1‑151 (MQ=255)
                 gccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgccg                                                                                                                                   >  1:720000/1‑150 (MQ=255)
                 gccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGt                                                                                                                                  >  1:303311/1‑151 (MQ=255)
                  ccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                                 >  1:784159/1‑151 (MQ=255)
                       tAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtc                                                                                                                            >  2:861496/1‑151 (MQ=255)
                            tGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                        >  1:509034/1‑150 (MQ=255)
                             ggACGGCGACCGCAAAAGGGTGAGTGATCCCCGGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGACCCGTGGTGTACCTCGACCCTTTAACCCGCTATCGGGATCATATGAAAAACCCCCCCGCCGATGCCCGCCGTTACAtcg                                                                                                                           <  1:728964/146‑1 (MQ=255)
                                 ggcgcccaaaaaaGGGGGAGGGATTCCCTGCGCGTACTGGAACGGGGCGCCGCGCGTCCCGGACCCGGGGTTTATCCCGAACCTTAAACCCCTTAACGGGATAAAATAAACACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGaaa                                                                                                                     <  2:737812/139‑1 (MQ=255)
                                  gCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                  >  2:410204/1‑150 (MQ=255)
                                     aCCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                              >  1:421566/1‑151 (MQ=255)
                                      ccGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                             >  2:198938/1‑151 (MQ=255)
                                      ccGAAAAAGGGGGGGTGTTCCCCTGCGCGTACTGGAACAGGCCCCCGCGCGTCCCGGAACCGTGGTGTACCTCAACCCGTAAACCCCCTATCAGAATAATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                               <  2:798584/149‑1 (MQ=255)
                                          aaaaTGGTGGTTGATTCCCTGCGCGTACTGAAACAGGTCGCCGCGCTCCCCGGACCCGTGTTTTATCTCGACCCGTTAACCCTCTATAAGAATCATATAACCACCCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                          <  2:626299/150‑1 (MQ=255)
                                          aaaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                          >  2:372016/1‑150 (MQ=255)
                                          aaaaTGGGGAGTGTTTCCCTGGGTTTACTGGAACAGCCCGCCGCGCGTCCCGAACCCGTGGTTTTTCCAGACCCTTTAACCCGTTACCAGAATTAAATTAAAACCCCCCCGGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                           <  1:633025/149‑1 (MQ=255)
                                           aaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                         >  1:877368/1‑150 (MQ=255)
                                           aaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                         >  1:36972/1‑150 (MQ=255)
                                              tGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                     >  1:541294/1‑151 (MQ=255)
                                                     tgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                              >  2:158800/1‑151 (MQ=255)
                                                          cccTGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGAACCGTGGTGTATCTCGAACCTTTAACCCGCTATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                          <  1:527714/150‑1 (MQ=255)
                                                          cccTGCGCGTACGGGAACAGTTCGCCGCGCGTCCCGGAGCCGTGTTGTTTCCCGACCCGTTAACCCGCTAGCAGTCTAATATGACCAACCCCCCCCCCGACGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTa                                                                                             <  2:128191/147‑1 (MQ=255)
                                                          cccGGCGCTTACTGAAACGGGTCGCCGCGCGTCCCGGACCCGGGGTTTATCCCAACCCTTTAACCCGCTATCAGGATAATAGGACAAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                         <  2:403514/151‑1 (MQ=255)
                                                            cTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                       >  1:812456/1‑151 (MQ=255)
                                                             tGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                       >  2:192701/1‑150 (MQ=255)
                                                                       ggAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                            >  1:781307/1‑151 (MQ=255)
                                                                             ggTCGCCGCGCTTCCCGGACCGGGGGTGTATCTCAAACCTTTAACCCGCTATCAGGATCATATAATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                         <  1:259856/148‑1 (MQ=255)
                                                                                cgccgcGCGTACGGGAACCGGGGTGTATCTGAACCCGTAAACCCGCTATAAGGATAATATGACCACCCCCCTCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                                       <  1:527691/147‑1 (MQ=255)
                                                                                  cggcgtgTCCCGGACCCGTGGTGTATCCCGCACGTTAACCCCGCTAACAGGATAATATAATACACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                                       <  2:258694/143‑1 (MQ=255)
                                                                                      gcgTCCCGGACCCGTGGTGTATCTCGACCCGTTAACCCGCTATCAGGATCATATGACCAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                  <  1:428165/146‑1 (MQ=255)
                                                                                        gTACCGGACCCGTGGTTTACCTCGACCCTTTAACCCGCTATCAGGATCATATGAACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                            <  1:748664/150‑1 (MQ=255)
                                                                                        gTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                              >  2:696803/1‑148 (MQ=255)
                                                                                             ggAACGGTGGTTTCTCAGGACCTGTTAACCCGCTATAAGGATCATATAAGCACCCCCCCCGCCTATTCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                         <  2:463231/132‑1 (MQ=255)
                                                                                             gaacccGTGGTGTATCCCAACCCTTAAACCCGCTATAAGAATAATATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                        <  2:187942/145‑1 (MQ=255)
                                                                                                     gttgtATCTAGACCCTTAAACCCGTTATCAGCATAATATAACCACCCCCCGCGAAGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                                <  2:63890/147‑1 (MQ=255)
                                                                                                     ggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                                >  2:655913/1‑149 (MQ=255)
                                                                                                       tttATCTCGACCCGTTAACCCGCTATCAGGATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                             <  1:601517/148‑1 (MQ=255)
                                                                                                       tgtATCTCGAACCGTTAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                             <  1:479543/150‑1 (MQ=255)
                                                                                                         tATCTCAACCCGTTAACCCGCTATCAGGATAATATGACAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                           <  1:647555/150‑1 (MQ=255)
                                                                                                          aTCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                         >  1:199910/1‑151 (MQ=255)
                                                                                                            cccGACCCGTAAAACCGCTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                        <  2:909437/148‑1 (MQ=255)
                                                                                                              cGAACCTTTAACCCGCTATCAGGATAATATCACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                     <  2:32338/151‑1 (MQ=255)
                                                                                                               gAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                    >  2:428827/1‑151 (MQ=255)
                                                                                                                     ttAACCCGCTATCAGGAATAAATAACAAACCCCCTCGCCAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACaa                                <  1:191751/149‑1 (MQ=255)
                                                                                                                             cccttaacAGGACAATATAAACACCCCCCCCGCAGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                           <  2:789088/142‑1 (MQ=255)
                                                                                                                             cGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                        >  1:192303/1‑151 (MQ=255)
                                                                                                                                 tATCAGAATCATATGACAAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                       <  1:701183/149‑1 (MQ=255)
                                                                                                                                     aGGATAAAATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                   <  1:629177/149‑1 (MQ=255)
                                                                                                                                         taatagGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc               <  1:192702/147‑1 (MQ=255)
                                                                                                                                         tCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCa              <  1:202669/150‑1 (MQ=255)
                                                                                                                                          cATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg            >  2:409734/1‑151 (MQ=255)
                                                                                                                                              tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc         >  1:323397/1‑150 (MQ=255)
                                                                                                                                                   aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg   >  1:897848/1‑151 (MQ=255)
                                                                                                                                                    acacCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGc  >  1:97240/1‑151 (MQ=255)
                                                                                                                                                    acacCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGc  >  2:728016/1‑151 (MQ=255)
                                                                                                                                                    acacCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGc  >  2:849414/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGC  >  NC_000913/1376358‑1376653

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