Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,741,995 G→A 14.5% intergenic (‑20/+393) PP_3305 ← / ← PP_3306 TerC family membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,741,9950GA14.5% 140.6 / 16.4 62intergenic (‑20/+393)PP_3305/PP_3306TerC family membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/28);  new base A (4/5);  total (29/33)
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.64e-01

CCGATTCCCGCGGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAA  >  NC_002947/3741847‑3742142
                                                                                                                                                     |                                                                                                                                                   
ccGATTCCCGCGGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGt                                                                                                                                                    >  2:935570/1‑150 (MQ=255)
 cGATTCCCGCGGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTa                                                                                                                                                   <  2:147660/150‑1 (MQ=255)
      cccGCGGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTgc                                                                                                                                              <  2:283836/150‑2 (MQ=255)
       ccGCGGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGaa                                                                                                                                             >  1:774274/1‑150 (MQ=255)
        cgcgGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGaac                                                                                                                                            <  1:826884/150‑1 (MQ=255)
             tGGAGCACACCTAAGTCAAGCGCGAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTa                                                                                                                                       <  1:492623/150‑1 (MQ=255)
              ggAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTAt                                                                                                                                      <  2:774274/150‑1 (MQ=255)
                cgatctccTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGa                                                                                                                                    <  1:306405/144‑1 (MQ=255)
                aGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGa                                                                                                                                    <  2:815236/150‑1 (MQ=255)
                aGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGa                                                                                                                                    <  1:736291/150‑1 (MQ=255)
                     acCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGATATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCg                                                                                                                               >  1:349887/1‑150 (MQ=255)
                     acCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCg                                                                                                                               >  2:392146/1‑150 (MQ=255)
                      ccTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGa                                                                                                                              >  2:826513/1‑150 (MQ=255)
                       cTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTAATAGATGTACATCTACGAcacatcg                                                                                                                             >  2:306405/1‑144 (MQ=255)
                        tAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACa                                                                                                                            <  1:935570/150‑1 (MQ=255)
                          aGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATg                                                                                                                          <  1:826513/150‑1 (MQ=255)
                            tCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGACtataca                                                                                                                             >  2:1086822/1‑140 (MQ=255)
                            tCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGACtataca                                                                                                                             <  1:1086822/145‑6 (MQ=255)
                                gcgcCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTacga                                                                                                                                    >  2:709562/1‑131 (MQ=255)
                                gcgcCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGa                                                                                                                                    <  1:709562/134‑1 (MQ=255)
                                gcgcCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTAACAAGAACGTCTTGGTAGAGATACGTAATAAGGGCAGACATGAGCCCCTCCTGTATAGGAGTAGTTGAACAAATATGACTCCGACATGGCGGta                                                                                                                    >  2:1097291/1‑150 (MQ=255)
                                  gcCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGtata                                                                                                                  >  1:1285690/1‑150 (MQ=255)
                                       cAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGt                                                                                                             <  1:861750/150‑1 (MQ=255)
                                                   aacaGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAAc                                                                                                                >  1:1115371/1‑135 (MQ=255)
                                                   aacaGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAAc                                                                                                                <  2:1115371/135‑1 (MQ=255)
                                                         gAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCa                                                                                           <  1:1150205/150‑1 (MQ=255)
                                                                gAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATaa                                                                                    >  2:387711/1‑150 (MQ=255)
                                                                     ccaaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAAt                                                                               >  1:92521/1‑150 (MQ=255)
                                                                       aaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAAttt                                                                             >  1:513236/1‑150 (MQ=255)
                                                                                 gTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGcgcg                                                                   >  2:1073737/1‑150 (MQ=255)
                                                                                      tACCAAGACGTTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAg                                                              >  1:674825/1‑150 (MQ=255)
                                                                                      tACCAAGACGTTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAg                                                              >  1:242220/1‑150 (MQ=255)
                                                                                       actaaGAACGTCTTGTTAGAGATACGGACTAAGGGCAGACATTTGCCCCTCTTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGTTTCCGGTGTCACTAATAAAGAATTTATTGAGCGCGAGTAgg                                                             <  2:1002795/147‑1 (MQ=255)
                                                                                          aaGAACGTCTTGGTAGAGATACGTACTNAGGGCAGACATGNGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTg                                                          >  1:1355131/1‑150 (MQ=255)
                                                                                            gAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCa                                                        >  1:214235/1‑150 (MQ=255)
                                                                                                gTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGAc                                                                                                                                   <  2:434061/71‑1 (MQ=37)
                                                                                                gTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGAc                                                                                                                                   >  1:434061/1‑71 (MQ=38)
                                                                                                    ttGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATAGCGGTATAACCTTCAAAGATGATGAGCTGGGCCGATAGAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGaa                                                 <  1:387711/150‑1 (MQ=255)
                                                                                                      ggTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGaaaa                                               <  1:270110/150‑1 (MQ=255)
                                                                                                          gagaTACGTACTAAGGGCAGATATGAGCCCCGCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGCGACGCTGAACTGCACTAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACAtt                                           <  2:349887/150‑1 (MQ=255)
                                                                                                          gagaTACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACAtt                                           <  2:674825/150‑1 (MQ=255)
                                                                                                                cGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACaaa                                     <  2:1060910/150‑1 (MQ=255)
                                                                                                                      aaGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGCCTCCGACATGGCGGTATAACCGTCAGACGCGCCCGGCGTCAAGAAGCTACGACGAATCCGGCGCGATTGGTTTGCAAGTAGAAAACATTTACAAATCTTGt                               <  1:1097291/150‑1 (MQ=255)
                                                                                                                          gCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATc                                                                          >  2:622569/1‑103 (MQ=255)
                                                                                                                          gCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATc                                                                          <  1:622569/103‑1 (MQ=255)
                                                                                                                          gCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCAt                           <  2:92521/150‑1 (MQ=255)
                                                                                                                                aTGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAg                     >  1:226289/1‑150 (MQ=255)
                                                                                                                                         tcctGGATAGTCGTATTTGAACAACTATGCCTCCGAGATCGCGGTATAACCGTCAGAGGGGACCGGCTTCACGAATAACGAATTTATCGAGCGCGAGTACGCTGCAAGTAGAAAACATGTACAAATATTGTCCATGACAAGCTAAAACGc            <  2:242482/150‑1 (MQ=255)
                                                                                                                                          cctGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCt           <  2:1355131/150‑1 (MQ=255)
                                                                                                                                           ctGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCtt          <  2:1285690/150‑1 (MQ=255)
                                                                                                                                             gTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCtaga        >  1:242482/1‑147 (MQ=255)
                                                                                                                                              tataGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTa                                                    <  1:1167528/105‑1 (MQ=255)
                                                                                                                                              tataGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTa                                                    >  2:1167528/1‑105 (MQ=255)
                                                                                                                                              tataGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCa       <  2:242220/150‑1 (MQ=255)
                                                                                                                                                  gTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGa   >  2:444990/1‑150 (MQ=255)
                                                                                                                                                   tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTc                                                                                                            >  2:98138/1‑44 (MQ=255)
                                                                                                                                                   tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTc                                                                                                            <  1:98138/44‑1 (MQ=255)
                                                                                                                                                   tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCagag                                                                                                        <  1:185669/48‑1 (MQ=255)
                                                                                                                                                   tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCagag                                                                                                        >  2:185669/1‑48 (MQ=255)
                                                                                                                                                   tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTc                              >  1:459721/1‑122 (MQ=255)
                                                                                                                                                   tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTc                              <  2:459721/122‑1 (MQ=255)
                                                                                                                                                   tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGaa  >  2:1180952/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGATTCCCGCGGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAA  >  NC_002947/3741847‑3742142

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: