Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,668,719 G→C 8.5% D55H (GAC→CAC)  rpsP → 30S ribosomal protein S16

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,668,7190GC 8.5% 185.1 / 7.1 59D55H (GAC→CAC) rpsP30S ribosomal protein S16
Reads supporting (aligned to +/- strand):  ref base G (30/24);  new base C (3/2);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

ATTCAGATTTGCAGTACCAGCTGCTGTATCAACTGCAGCACCGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGCCGTCTGAGCGTGTTGCTCAGCTGCTG  >  NC_002947/1668572‑1668859
                                                                                                                                                    |                                                                                                                                             
aTTCAGATTTGCAGTACCAGCTGCTGTATCAACTGCAGCACCGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGa                                                                                                                                              >  4:130027/1‑149 (MQ=255)
     gATTTGCAGTACCAGCTGCTGTATCAACTGCAGCACCGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGcc                                                                                                                                         <  3:285260/149‑1 (MQ=255)
      aTTTGCAGTACCAGCTGCTGTATCAACTGCAGCACCGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCg                                                                                                                                        >  3:233858/1‑149 (MQ=255)
                ccAGCTGCTGTATCAACTGCAGCACCGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGa                                                                                                                               <  3:266458/148‑1 (MQ=255)
                ccAGCTGCTGTATCAACTGCAGCACCGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAg                                                                                                                              >  7:92078/1‑149 (MQ=255)
                       ctgTATCAACTGCAGCACCGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAg                                                                                                                       >  3:63366/1‑149 (MQ=255)
                       ctgTATCAACTGCAGCACCGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAg                                                                                                                       >  6:165107/1‑149 (MQ=255)
                              aaCTGCAGCACCGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTttctt                                                                                                                <  3:130027/149‑1 (MQ=255)
                                cTGCAGCACCGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCCGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGTGCGTAGGTttcttc                                                                                                              >  6:262772/1‑148 (MQ=255)
                                   cagcaCCGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAAcc                                                                                                          <  7:179538/149‑1 (MQ=255)
                                      caCCGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCaa                                                                                                       <  5:165107/149‑1 (MQ=255)
                                        ccGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCaa                                                                                                            >  4:121214/1‑142 (MQ=255)
                                        ccGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCaa                                                                                                            <  3:121214/142‑1 (MQ=255)
                                                  gAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACGACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGgcg                                                                                           >  7:306939/1‑149 (MQ=255)
                                                    aCGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGtt                                                                                                                                 <  6:343550/110‑1 (MQ=255)
                                                    aCGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGtt                                                                                                                                 >  5:343550/1‑110 (MQ=255)
                                                     cGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTGTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAAGCGCTGCTGGCGcc                                                                                         >  6:290668/1‑148 (MQ=255)
                                                          ttCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCCGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGt                                                                                    <  5:262772/148‑1 (MQ=255)
                                                            cacaTGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGtt                                                                                   >  8:222162/1‑147 (MQ=255)
                                                              caTGGTAACCATTCGTCTGGCCCGTTGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCaa                                                                                                            <  4:64973/120‑1 (MQ=255)
                                                              caTGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGAACGTGACCAACTCGAGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCaa                                                                                                            >  3:64973/1‑120 (MQ=255)
                                                               aTGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCt                                                                              >  7:129387/1‑149 (MQ=255)
                                                                  gTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTc                                                                           <  8:92078/149‑1 (MQ=255)
                                                                       cATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTc                                                                       >  1:342350/1‑148 (MQ=255)
                                                                        aTTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAg                                                                                >  2:119371/1‑138 (MQ=255)
                                                                        aTTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAg                                                                                <  1:119371/138‑1 (MQ=255)
                                                                        aTTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCaa                                                                     >  7:130462/1‑149 (MQ=255)
                                                                        aTTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGCCaa                                                                     <  3:208796/149‑1 (MQ=255)
                                                                                 gCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCgg                                                                         <  1:246474/137‑1 (MQ=255)
                                                                                  cccGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAgcgc                                                           <  4:63366/149‑1 (MQ=255)
                                                                                             tCGAAAAAGCGCCCATTCTACCAGCTGTCCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCgg                                                                         <  5:63041/124‑1 (MQ=255)
                                                                                             tCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTg                                                >  6:181332/1‑149 (MQ=255)
                                                                                             ctgAAAAAGCGACCATTCTACCACCTGACCGTTCCCAACTCGCGTAACGCCAGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACt                                                 <  7:318448/146‑1 (MQ=255)
                                                                                               gAAAAAGCGCCCATTCTACCACCTGACCGTTGACCAACTCGCGTAACGCCCGTGCCGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTgg                                               >  1:26634/1‑149 (MQ=255)
                                                                                                aaaaaGCGCCCATTCTACCCCCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCACCCCAAGCGCTGCTGGCGCCTAAGTTAAGCAGTCGGGACACCACGAGCTCGTACCCTACTgcta                                             >  8:318448/1‑146 (MQ=255)
                                                                                                    aGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGc                                         <  7:222162/149‑1 (MQ=255)
                                                                                                     gcgcCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGcc                                                                                         <  5:228960/100‑1 (MQ=255)
                                                                                                     gcgcCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGcc                                                                                         >  6:228960/1‑100 (MQ=255)
                                                                                                       gcCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGa                                           >  6:201874/1‑144 (MQ=255)
                                                                                                       gcCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGa                                           <  5:201874/144‑1 (MQ=255)
                                                                                                       gcCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAg                                      <  2:342350/149‑1 (MQ=255)
                                                                                                         ccATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTgg                                               <  8:308918/138‑1 (MQ=255)
                                                                                                         ccATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTgg                                               >  7:308918/1‑138 (MQ=255)
                                                                                                                accaccTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAg                             >  5:154103/1‑149 (MQ=255)
                                                                                                                   accTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGCCg                          >  7:155941/1‑149 (MQ=255)
                                                                                                                    ccTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGCCGt                         >  5:134832/1‑149 (MQ=255)
                                                                                                                      tGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGCCGTCt                       >  1:27073/1‑149 (MQ=255)
                                                                                                                        aCCGTGACCAACTCGCGTAACGCCCGTGACTGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGCCGTCTGa                     <  8:306939/149‑1 (MQ=255)
                                                                                                                        aCCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGt                                   >  6:261680/1‑135 (MQ=255)
                                                                                                                        aCCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGt                                   <  5:261680/135‑1 (MQ=255)
                                                                                                                          cGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGCCGTCTGAGc                   >  3:155413/1‑149 (MQ=255)
                                                                                                                              gACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGCCGTCTGAGCgtg                >  7:225635/1‑149 (MQ=255)
                                                                                                                               aCCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGAGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGCCGTCTGAGCgtg                <  2:27073/148‑1 (MQ=255)
                                                                                                                                     caacgaaaCGGCCGTCACGGGCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGCCGTCTGAGCGTGTTGCt           >  6:141932/7‑147 (MQ=255)
                                                                                                                                         gaaaCGGCCGTCACGGGCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGCCGTCTGAGCgtgt               <  3:136016/137‑1 (MQ=255)
                                                                                                                                         gaaaCGGCCGTCACGGGCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGCCGTCTGAGCgtgt               >  4:136016/3‑139 (MQ=255)
                                                                                                                                              gCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGCCGTCTGAGCGTGTTGCTCAGctgctg  >  4:340785/1‑147 (MQ=255)
                                                                                                                                                ccGTCACGGGCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGcc                           <  5:325721/120‑1 (MQ=255)
                                                                                                                                                ccGTCACGGGCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGcc                           >  6:325721/1‑120 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
ATTCAGATTTGCAGTACCAGCTGCTGTATCAACTGCAGCACCGACTATAGGAACGATGTTCACATGGTAACCATTCGTCTGGCCCGTGGCGGCTCGAAAAAGCGCCCATTCTACCACCTGACCGTGACCAACTCGCGTAACGCCCGTGACGGCCGTTTCGTTGAGCGCGTAGGTTTCTTCAACCCAATCGCTGCTGGCGCCGAAGTTAAGCTGTCGGTCAACCAAGAGCGCGTCACCTACTGGCTGAGCCAGGGTGCACAGCCGTCTGAGCGTGTTGCTCAGCTGCTG  >  NC_002947/1668572‑1668859

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