Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,245,2340GT43.9% 74.0 / 35.0 57L25F (TTG→TTTycgYuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (8/24);  new base T (25/0);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 1.55e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATG  >  NC_000913/1245089‑1245379
                                                                                                                                                  |                                                                                                                                                 
atatCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTgc                                                                                                                                              <  2:691811/151‑1 (MQ=255)
      gATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTAt                                                                                                                                        <  2:897844/151‑1 (MQ=255)
      gATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTAt                                                                                                                                        <  2:365557/151‑1 (MQ=255)
       aTAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTg                                                                                                                                      <  2:309676/151‑1 (MQ=255)
         aaaGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTg                                                                                                                                      <  2:482403/150‑1 (MQ=255)
          aaGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGGGCTTAtt                                                                                                                                       >  1:480808/1‑148 (MQ=255)
                gTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTTTTTAAAAttt                                                                                                                               >  1:611064/1‑150 (MQ=255)
                   cccGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTGAACATTTAGa                                                                                                                            >  2:988174/1‑150 (MQ=255)
                     cGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAAtt                                                                                                                         <  1:655666/151‑1 (MQ=255)
                      ggTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGGTGCGCTTATTTAAAATTTAGaa                                                                                                                           >  1:229895/1‑148 (MQ=255)
                          acaacaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCa                                                                                                                    <  2:189758/151‑1 (MQ=255)
                             acaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTGAACATTTAGAAATGAt                                                                                                                      >  1:46247/1‑146 (MQ=255)
                              cagcaGAATTCTTTTTTAATTATTGTTTCGTTGTTTTTGTTTTGATTCAGAAAATAAATATTTTAGACAGTCGTCAGTCAGGTGTTAAAAATGGGGCGAATAGTGATGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGAACAATc                                                                                                                 <  1:103109/147‑1 (MQ=255)
                               aacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTATTTTAAATTTTGAATTGATCaa                                                                                                                   >  2:951480/1‑147 (MQ=255)
                                acaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCtt                                                                                                              <  1:827671/151‑1 (MQ=255)
                                 caGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTTAACATTTTGAATTGATTAATTCtt                                                                                                              >  1:877396/1‑150 (MQ=255)
                                             ttAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGGTTATTTAACATTTTGAATTGATCAATCCTTCTTGCAGATTaa                                                                                                  >  2:1104773/1‑150 (MQ=255)
                                               aaTTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGcc                                                                                               <  1:494006/151‑1 (MQ=255)
                                                     tGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGAGCTTTTTTTATACTTATAATTGTTTTATCCCTTTCTCAGATTAATCCcaaa                                                                                           >  2:1099718/1‑149 (MQ=255)
                                                        ttCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTATTTAACATTTAAAATTTATTAATCCTTCTTGCCGATTAAGCCCAAATAacc                                                                                      >  2:1099082/1‑151 (MQ=255)
                                                        ttCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGGCTTTTTAACACTTTGAATTGATCAATACTTTTGGGAGAATAAACCaaaaaaacc                                                                                      >  1:406355/1‑146 (MQ=255)
                                                          cGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGGGCTTTTTTAACAATTTGAATTGATCAATCCTTCTTGGAGATTAAGACAAAAAAACCg                                                                                     >  1:223788/1‑150 (MQ=255)
                                                          cGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCg                                                                                     <  1:241810/150‑1 (MQ=255)
                                                                 tttGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGGGCTTATTTAACATTTAGAATTGATCAATCCCTCTTGCAGAGTAAGCCAAAATACCCGTTCg                                                                                 >  1:905486/1‑147 (MQ=255)
                                                                   tGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTATTTAACAATTATAATTGATCAATCCCTCTGGCAGATTAAAGCAAAAAACCCGTTTGGCaa                                                                             >  2:613142/1‑149 (MQ=255)
                                                                        tgATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCCGCGCCTATTGAATAATTTGCCTTAATAACACCTTTTGGCAGCTCAAACCAAACAAACCGTTTCGCAATGttt                                                                        >  1:1152849/1‑149 (MQ=255)
                                                                           ttCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCCTATTTAACATTTAGAATTTATCAAACCTTCTTGCAGATTAAACCAAAATAACCCTTCGGCAATGCTTCAGtt                                                                   >  2:1146773/1‑151 (MQ=255)
                                                                             cAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTa                                                                  <  1:483468/150‑1 (MQ=255)
                                                                             cAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTAc                                                                 <  2:297629/151‑1 (MQ=255)
                                                                              aGAAAATAAATATTTTAGACAGTAGACCGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGATTTTGGTGCGCGTTTTTAAAATTTAGAATTTTTTAATTCTTTTGGCAGGTTAAGTCAAAAAAACCGTTTGGGaaagatttcgtttc                                                                 >  1:819535/1‑138 (MQ=255)
                                                                               gAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGCTTTTGCTGGGACTAATGAAAAATTAGAATTGATCAAACCTTTTTGGGGAATAAGCCAAAATAACAGTTCGGCCATGCCTCacatta                                                                 >  1:517980/1‑144 (MQ=255)
                                                                                  aataaataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTCTTTTCTGCGCCTATTTAACATTTACAATTGATCAATTCTTTTGGGAGAAAAAGCCCAAATACCCCATCGGCCATGGTTccgatttgc                                                               >  2:826020/1‑140 (MQ=255)
                                                                                      aataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTTCGCTTTTTTAACATTTTGAATTGATAAATTCTTCTGGCCGATTAAGCCAAAAAAACCGTTCGGCAATGCTTCTGTTACGCGAAGAgg                                                        >  1:590943/1‑151 (MQ=255)
                                                                                           tttAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATg                                                    <  1:1151760/150‑1 (MQ=255)
                                                                                            ttAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACAATTAGAATTTAACAATCCTTCTGGCAGATTAAGCCAAAAAAACCGTTCGGCAATGCTTTAGTTACGCGGAGAGGAATgg                                                   >  1:844300/1‑150 (MQ=255)
                                                                                             tAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTTTTTAAAATTTAGAATTGTTTCAACCCTCTGGCAGATTAAACCAAAATACCCGTTCGGCAATGCTTTAGTTTCCCGGAGAGGGATgg                                                   >  1:1018562/1‑149 (MQ=255)
                                                                                             tAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGGTGCGCGTTTTGAACAATTGGCATTGATTAAATCTTCTTGCGGAATTAACCAAAATACCCGTTCGGCAAGGCTTCAGTTACGCGGAGAGGGAttg                                                   >  1:165593/1‑147 (MQ=255)
                                                                                             tAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGt                                                  <  1:757516/150‑1 (MQ=255)
                                                                                              aGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTGAACATTTAGAATTTATCAATCCTTCTTGCAGATTAAACCAAAATACCCCTTCGGCAAAGCTTCAGTTACGCGAAGAgggaagg                                                   >  2:730869/1‑143 (MQ=255)
                                                                                              aGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCt                                                <  2:702112/151‑1 (MQ=255)
                                                                                               gACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGGTTTTGGCGCGCTTTTTTAACAATTTAAATTGATCAATCCTTCTTGGAGATTAAGCCAAAAAAACCCTTCGGCCATGCTTCAGTTACGCGgagaaggat                                                     >  1:342690/1‑140 (MQ=255)
                                                                                               gACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTTTTTAACATTTAGAATTGATCAATCCTTTTGGCAGATTAAACCAAAAAACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGGATGGTCt                                                >  2:513620/1‑150 (MQ=255)
                                                                                                aCAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTTTTTAACATTTAGAATTGATCAATCCTTCTGGGAGATTAAACCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCt                                                >  2:425547/1‑149 (MQ=255)
                                                                                                aCAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGa                                              <  2:606391/151‑1 (MQ=255)
                                                                                                 cAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTTCGCTTATTGGACATTTTGAATTGGATAAACCTTTTGGGAGAATAAGCCAAAAAAACCGCTCTGGAATGCCTCAGTTACGCCAAGGGGAAtag                                                   >  2:486102/1‑143 (MQ=255)
                                                                                                 cAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTTTTTAAAATTTTGAATTGATCAATCCTTCTTGCAGATTAAGCCAAAAAAACCCTTCGGCAATGCTTCAGTTACGCGGAGAGGAAAgg                                                   >  1:867352/1‑145 (MQ=255)
                                                                                                 cAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGaa                                             <  2:508079/151‑1 (MQ=255)
                                                                                                  aGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGCTGCGCTTATTTAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAAAACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGaa                                             >  2:765118/1‑150 (MQ=255)
                                                                                                          gTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCtttttt                                    <  2:420143/151‑1 (MQ=255)
                                                                                                           tCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTTGGTGCGCTTTTTTAAAAATTTGAATTGATCAAATCTTCTGGGAGGTTAAACCAAAAAACCCGTTCGGGCAAGCTTCAGTTTAGCGAAGAGGAAAGGGCTTAAtttttttttt                                   >  1:932981/1‑151 (MQ=255)
                                                                                                            cAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTg                                  <  1:947398/151‑1 (MQ=255)
                                                                                                                    aaaaaTGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTaa                          <  2:642447/151‑1 (MQ=255)
                                                                                                                    aaaaaTGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTa                           <  1:820772/150‑1 (MQ=255)
                                                                                                                        aTGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATtc                      <  1:949308/151‑1 (MQ=255)
                                                                                                                                    gTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCt          >  2:723729/1‑151 (MQ=255)
                                                                                                                                        ggCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTgg        >  1:622786/1‑149 (MQ=255)
                                                                                                                                             gTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTAtg  <  1:1036714/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGTTTTGGCTGCGCTTATTGAACATTTAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATG  >  NC_000913/1245089‑1245379

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