Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,194,721 A→G F187F (TTT→TTC U712_10995 ← Anthranilate synthase component 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,194,7210AG98.4% 189.6 / ‑4.7 61F187F (TTT→TTCU712_10995Anthranilate synthase component 1
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (19/41);  total (19/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

CTGGTTTTTCATAATTAGAAGATACTGTCTCAAAGCTGGGTGTCTTATATGAATCAGCAGAAAGAAACAGTTCTTTTATGTTTTTTTGATCCATCATTTTTTCTATGAGATCTTGAAGCTCTTGATGATTTTGATGAAATACATCGATTTTTTCGTTTTCTGTTTCCTGTCCATTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTAGTTTCATGATCATATGCAATTAATGTCCGGCAAACAAACAGCATACACTTTTCCATTCC  >  CP006881/2194585‑2194852
                                                                                                                                        |                                                                                                                                    
ccgatttttCATAATTAGAAGATACTGTCTCAAAGCTGGGTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGgaa                                                                                                                                    <  4:392735/135‑1 (MQ=255)
ccgatttttCATAATTAGAAGATACTGTCTCAAAGCTGGGTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGgaa                                                                                                                                    <  1:354099/135‑1 (MQ=255)
                      tACTGTCTCAAAGCTGGGTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGTTTTTGAAGATACAGATGATTTTGATGGAATACATCCATTTTTTCTGttttt                                                                                                              >  4:267897/1‑138 (MQ=255)
                       aCTGTCTCAAAGCTGGGTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTg                                                                                                            >  4:287536/1‑139 (MQ=255)
                        cTGTCTCAAAGCTGGGTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGt                                                                                                           <  8:295377/139‑1 (MQ=255)
                        cTGTCTCAAAGCTGGGTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGt                                                                                                           <  3:65584/139‑1 (MQ=255)
                        cTGTCTCAAAGCTGGGTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGt                                                                                                           <  2:47849/139‑1 (MQ=255)
                         tGTCTCAAAGCTGGGTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGtt                                                                                                          >  5:119233/1‑139 (MQ=255)
                         tGTCTCAAAGCTGGGTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGtt                                                                                                          >  1:434897/1‑139 (MQ=255)
                           tctcAAAGCTGGGTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGttt                                                                                                         <  4:354680/138‑1 (MQ=255)
                           tctcAAAGCTGGGTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGttt                                                                                                         >  3:303000/1‑138 (MQ=255)
                               aaaGCTGGGTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCt                                                                                                      <  3:287536/137‑1 (MQ=255)
                                     gggTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCATCTCCagt                                                                                               >  2:245082/1‑136 (MQ=255)
                                      ggTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCagt                                                                                               >  8:317002/1‑135 (MQ=255)
                                      ggTGTCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCagt                                                                                               <  4:81931/137‑3 (MQ=255)
                                      ggTGTCTTATATGAATCAGCAGAAAGAAACAGTTCTTTTATGTTTTTTTGATCCATCATTTTTTCTATGAGATCTTGAAGCTCTTGATGATTTTGATGAAATACATCGATTTTTTCGTTTTCTGTTTCCTGTCCAtt                                                                                               <  6:307581/137‑1 (MQ=255)
                                         gtCTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGcc                                                                                          >  4:158631/1‑139 (MQ=255)
                                           cTTGTATGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCtt                                                                                        >  1:105018/1‑139 (MQ=255)
                                              gtatGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCtt                                                                                        <  6:300157/135‑1 (MQ=255)
                                              gtatGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCtt                                                                                        <  2:111982/135‑1 (MQ=255)
                                               tatGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTGTTTCCTCTCCAGTGAGCCTTGCata                                                                                   <  3:46234/139‑1 (MQ=255)
                                                atGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCGTGCtta                                                                                   >  6:87675/1‑139 (MQ=255)
                                                 tGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAGTACATACATTTTTTCGTTTTTTGTTTACTCTCCAGTGAGCCTTGCatat                                                                                  >  1:309692/1‑139 (MQ=255)
                                                 tGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCata                                                                                   <  7:324871/138‑1 (MQ=255)
                                                 tGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCata                                                                                   <  4:212807/138‑1 (MQ=255)
                                                 tGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCata                                                                                   <  3:74546/138‑1 (MQ=255)
                                                 tGAATCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTATCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCata                                                                                   <  4:164531/138‑1 (MQ=255)
                                                    aTCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGt                                                                                                           <  3:128409/111‑1 (MQ=255)
                                                    aTCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGt                                                                                                           >  4:128409/1‑111 (MQ=255)
                                                    aTCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTg                                                                                >  3:285953/1‑138 (MQ=255)
                                                     tcatcaGAAACAAACCTCTCTTTTATGTTTTTTTCGTCCATCATTTATTCAATGATATTTTGAAGCTACAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGa                                                                              <  2:327908/139‑1 (MQ=255)
                                                     tCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGa                                                                              <  7:366384/139‑1 (MQ=255)
                                                     tCAGCAGAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGa                                                                              <  1:150761/139‑1 (MQ=255)
                                                        gcagAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATa                                                                            >  1:173846/1‑138 (MQ=255)
                                                         cagAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATaaa                                                                          <  2:248665/139‑1 (MQ=255)
                                                         cagAAAGAAACAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATaaa                                                                          >  3:448841/1‑139 (MQ=255)
                                                            aaagaaaCAGCTCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTg                                                                       >  5:382886/1‑139 (MQ=255)
                                                                      ctCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGtttttg                                                             >  5:301578/2‑138 (MQ=255)
                                                                      ctCTTTTATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGtttttg                                                             <  4:133300/138‑2 (MQ=255)
                                                                           ttATGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGttt                                                         >  3:54322/1‑138 (MQ=255)
                                                                             aTGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCa                                                       <  1:245082/138‑1 (MQ=255)
                                                                             aTGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCa                                                       <  7:314623/138‑1 (MQ=255)
                                                                             aTGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCa                                                       <  5:436299/138‑1 (MQ=255)
                                                                             aTGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCa                                                       <  3:277760/138‑1 (MQ=255)
                                                                             aTGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCa                                                       <  2:17997/138‑1 (MQ=255)
                                                                             aTGTTTTTTTGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTAGACGTTTTTGGTTTCa                                                       <  2:19554/138‑1 (MQ=255)
                                                                                     ttGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATTTTGGATAAAGTGGACGTTTTTGGTTTCATGATCata                                               <  1:14417/138‑1 (MQ=255)
                                                                                     ttGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCata                                               <  1:318354/138‑1 (MQ=255)
                                                                                     ttGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCata                                               <  4:303000/138‑1 (MQ=255)
                                                                                     ttGGTCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCata                                               <  2:434897/138‑1 (MQ=255)
                                                                                       ggtCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCatacg                                             <  4:60680/136‑3 (MQ=255)
                                                                                        gtCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCATACGc                                            >  3:85370/2‑138 (MQ=255)
                                                                                         tCCATCATTTTTTCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTACAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCATACGCaa                                          <  3:267897/139‑1 (MQ=255)
                                                                                                   tttCAATGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCATACGCAATTAATGTcc                                 <  5:106504/138‑1 (MQ=255)
                                                                                                     tcaaTGAGATTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCATACGCAATTAATGTCCgg                               <  7:236612/135‑1 (MQ=255)
                                                                                                            gaTTTTGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCATACGCAATTAATGTCCGGcaaacaaa                       <  4:427611/139‑1 (MQ=255)
                                                                                                             attttGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCATACGCAATTAATGTCCGGcaaacaaa                       <  7:317002/135‑1 (MQ=255)
                                                                                                              ttttGAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCATACGCAATTAATGTCCGGCaaacaaac                      <  3:429691/136‑1 (MQ=255)
                                                                                                                  gAAGCTCCAGATGATTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCATACGCAATTAATGTCCGGCAAACAAACAGCa                  <  2:105018/138‑1 (MQ=255)
                                                                                                                            atgatTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCATACGCAATTAATGTCCGGCAAACAAACAGCATACACTTTTcc       <  3:248635/139‑1 (MQ=255)
                                                                                                                               atTTTGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCATACGCAATTAATGTCCGGCAAACAAACAGCATACACTTTTCCAt     <  4:285953/138‑1 (MQ=255)
                                                                                                                                  ttGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCATACGCAATTAATGTCCGGCAAACAAACAGCATACACTTTTCCAtgtc  <  3:158631/138‑4 (MQ=255)
                                                                                                                                  ttGATGGAATACATCCATTTTTTCGTTTTTTGTTTCCTCTCCAGTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTGGTTTCATGATCATACGCAATTAATGTCCGGCAAACAAACAGCATACACTTTTCCAtgtc  <  2:173846/138‑4 (MQ=255)
                                                                                                                                        |                                                                                                                                    
CTGGTTTTTCATAATTAGAAGATACTGTCTCAAAGCTGGGTGTCTTATATGAATCAGCAGAAAGAAACAGTTCTTTTATGTTTTTTTGATCCATCATTTTTTCTATGAGATCTTGAAGCTCTTGATGATTTTGATGAAATACATCGATTTTTTCGTTTTCTGTTTCCTGTCCATTGAGCCTTGCATATTGGATAAAGTGGACGTTTTTAGTTTCATGATCATATGCAATTAATGTCCGGCAAACAAACAGCATACACTTTTCCATTCC  >  CP006881/2194585‑2194852

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