Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,026,285 T→C 7.1% intergenic (‑39/+140) PP_3550 ← / ← PP_3551 transcriptional regulator EmrR/LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,026,2850TC 7.1% 167.6 / 2.5 57intergenic (‑39/+140)PP_3550/PP_3551transcriptional regulator EmrR/LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (24/28);  new base C (2/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

CCTCTGACTCCAGGTGCCAGTCAAGAATGCGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGGCCGCCCCTGCCAGGGCCTGCGTAAGGGCGGTTAAC  >  NC_002947/4026138‑4026429
                                                                                                                                                   |                                                                                                                                                
ccTCTGACTCCAGGTGCCAGTCAAGAATGCGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTa                                                                                                                                                <  1:770108/150‑1 (MQ=255)
 ctctGACTCCAGGTGCCAGTCAAGAATGCGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAg                                                                                                                                               <  1:688704/150‑1 (MQ=255)
 ctctGACTCCAGGTGCCAGTCAAGAATGCGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAg                                                                                                                                               >  1:454604/1‑150 (MQ=255)
    tGACTCCAGGTGCCAGTCAAGAATGCGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCa                                                                                                                                            <  1:844844/150‑1 (MQ=255)
     gACTCCAGGTGCCAGTCAAGAATGCGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAg                                                                                                                                           >  2:850967/1‑150 (MQ=255)
          cAGGTGCCAGTCAAGAATGCGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAAt                                                                                                                                      <  1:399216/150‑1 (MQ=255)
                ccAGTCAAGAATGCGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGACAATTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCaaa                                                                                                                                <  1:573954/150‑1 (MQ=255)
                   gTCAAGAATGCGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGcc                                                                                                                             >  1:387640/1‑150 (MQ=255)
                    tCAAGAATGCGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCt                                                                                                                            >  2:11164/1‑150 (MQ=255)
                         aaTGCGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCAtt                                                                                                                       <  1:516920/150‑1 (MQ=255)
                          aTGCGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTa                                                                                                                      <  1:81610/150‑1 (MQ=255)
                             cGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGc                                                                                                                   >  2:582951/1‑150 (MQ=255)
                                   ttGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTAt                                                                                                             >  1:89784/1‑150 (MQ=255)
                                            gCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATACTc                                                                                                    <  1:850967/150‑1 (MQ=255)
                                              tGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCgt                                                                                                  <  1:286933/150‑1 (MQ=255)
                                                 gcgGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCtt                                                                                                                           <  1:813031/122‑1 (MQ=255)
                                                 gcgGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCtt                                                                                                                           >  2:813031/1‑122 (MQ=255)
                                                     ccAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCt                                                                                                                                                  >  2:788809/1‑95 (MQ=255)
                                                     ccAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCt                                                                                                                                                  <  1:788809/95‑1 (MQ=255)
                                                     ccAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTc                                                                                                    >  2:258277/1‑141 (MQ=255)
                                                     ccAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTc                                                                                                    <  1:258277/141‑1 (MQ=255)
                                                       aaGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCTCTGGTTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGc                                                                                         >  1:225604/1‑150 (MQ=255)
                                                       aaGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGc                                                                                         >  1:217991/1‑150 (MQ=255)
                                                          cAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCgg                                                                                      >  1:738778/1‑150 (MQ=255)
                                                                gCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCAc                                                                                >  2:246161/1‑150 (MQ=255)
                                                                          cAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCa                                                                      <  2:89784/150‑1 (MQ=255)
                                                                           aGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAg                                                                     <  1:11164/150‑1 (MQ=255)
                                                                              tCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCgt                                                                                                  <  1:219507/118‑1 (MQ=255)
                                                                              tCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCgt                                                                                                  >  2:219507/1‑118 (MQ=255)
                                                                              tCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCg                                                                                       <  2:660577/129‑1 (MQ=255)
                                                                              tCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCg                                                                                       >  1:660577/1‑129 (MQ=255)
                                                                                 gAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGgcgcgc                                                               <  2:387640/150‑1 (MQ=255)
                                                                                       tttCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCg                                                         <  1:246161/150‑1 (MQ=255)
                                                                                              ggAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGc                                                       >  2:526362/1‑145 (MQ=255)
                                                                                              ggAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGc                                                       <  1:526362/145‑1 (MQ=255)
                                                                                               gagaAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCAcc                                                                               <  1:717274/120‑1 (MQ=255)
                                                                                               gagaAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCAcc                                                                               >  2:717274/1‑120 (MQ=255)
                                                                                                  aaaTGGTGCATCGTTGACTGCTCAGATACGTGTGTAGAAATATTACTGCGGAGGCAGGTAATGTCAAAGCCTTCATGAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGa                                              <  2:692002/150‑1 (MQ=255)
                                                                                                         ccATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATc                                                  <  1:163675/139‑1 (MQ=255)
                                                                                                         ccATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATc                                                  >  2:163675/1‑139 (MQ=255)
                                                                                                               ttGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGgccgc                                 >  1:795752/1‑150 (MQ=255)
                                                                                                                tGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGgccgcc                                <  2:217991/150‑1 (MQ=255)
                                                                                                                          gTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGGCCGCCCCTGCCAggg                      <  1:582951/150‑1 (MQ=255)
                                                                                                                                tgtgtAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGccc                                                             >  1:195076/1‑105 (MQ=255)
                                                                                                                                tgtgtAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGccc                                                             <  2:195076/105‑1 (MQ=255)
                                                                                                                                  tgtAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGGCCGCCCCTGCCAGGGCCTGCGTa              >  1:428874/1‑150 (MQ=255)
                                                                                                                                   gtAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGGCCGCCCCTGCCAGGGCCTGCGTaa             >  2:134081/1‑150 (MQ=255)
                                                                                                                                     aGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGGCCGCCCCTGCCAGGGCCTGCGTAAgg           <  2:795752/150‑1 (MQ=255)
                                                                                                                                      gAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGGCCGCCCCTGCCAGGGCCTGCGTAAggg          >  1:725196/1‑150 (MQ=255)
                                                                                                                                         atatTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGcc                                                      >  2:622237/1‑103 (MQ=255)
                                                                                                                                         atatTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGcc                                                      <  1:622237/103‑1 (MQ=255)
                                                                                                                                           atTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTCCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGGCCGCCCCTGCCAGGGCCTGCGTAAGGGCGGtt     <  2:454604/150‑1 (MQ=255)
                                                                                                                                            ttACTGCCTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGGCCGCCCCTGCCAGGGCCTGCGTa              <  1:605875/140‑1 (MQ=255)
                                                                                                                                            ttACTGCCTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGGCCGCCCCTGCCAGGGCCTGCGTa              >  2:605896/1‑140 (MQ=255)
                                                                                                                                            ttACTGCCTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGGCCGCCCCTGCCAGGGCCTGCGTa              >  2:605875/1‑140 (MQ=255)
                                                                                                                                            ttACTGCCTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGGCCGCCCCTGCCAGGGCCTGCGTa              <  1:605896/140‑1 (MQ=255)
                                                                                                                                              aCTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGGCCGCCCCTGCCAGGGCCTGCGTAAGGGCGGTTAAc  <  2:725196/150‑1 (MQ=255)
                                                                                                                                                tGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGgc                                    <  2:673005/114‑1 (MQ=255)
                                                                                                                                                tGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGgc                                    >  1:673005/1‑114 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCTCTGACTCCAGGTGCCAGTCAAGAATGCGGTCTTTGAGGATCGCTGCGCGGCCAAGCAACATGCCGATGGCGCAGGTCTGAAAGTTTTCGGGGGAGAAATGGGCCATCGTTGACTGCTCAGTTAAGTGTGTAGAAATATTACTGCTTAGGCAGTTAATGTCAAAGCCTTCATTAGGCAGCTATGAAATATTCGTGTTTTCTGCCGGTCCCACCGCCCCCGCAGGCGCGCCCTGCGGCCCATCGCGACACACAAGGCCGCCCCTGCCAGGGCCTGCGTAAGGGCGGTTAAC  >  NC_002947/4026138‑4026429

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