Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,381,985 C→T 6.2% I509I (ATC→ATT CJ019_05063 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,381,9850CT 6.2% 211.0 / 3.2 64I509I (ATC→ATTCJ019_05063hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/26);  new base T (2/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

GTGGAGTTTGGCGGCTCGCTCAAGGCCGAGCACGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCGGTGTGCCCGTCCAAGGG  >  ACB122/5381851‑5382119
                                                                                                                                      |                                                                                                                                      
gTGGAGTTTGGCGTCTCGCTCAAGGCCGAGCACGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATc                                                                                                                                        <  4:293120/135‑1 (MQ=255)
    aGTTTGGCGGCTCGCTCAAGGCCGAGCACGGCACCGGGCGCAACCTGGCGCCCTTCCTGGAGCTGGAGTGGGGGCACGACGCCGACCAGCGGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCAACCTcc                                                                                                                                    >  2:207432/1‑134 (MQ=255)
          gcggcTCGCTCAAGGCCGAGCACGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCg                                                                                                                              <  8:8873/135‑1 (MQ=255)
           cggcTCGCTCAAGGCCGAGCACGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCg                                                                                                                              <  4:99221/134‑1 (MQ=255)
           cggcTCGCTCAAGGCCGAGCACGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGGGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCg                                                                                                                              <  5:252611/134‑1 (MQ=255)
            ggcTCGCTCATGGCCGAGCACGGCACCGGGCGCAACATGGCGCCCGTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGAc                                                                                                                            <  5:102216/135‑1 (MQ=255)
            ggcTCGCTCAAGGCCGAGCACGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGAc                                                                                                                            <  5:297416/135‑1 (MQ=255)
                cgctcAAGGCCGAGCACGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCGGGAGTGGGGGCACGACGCCTACCAGGTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACgtggt                                                                                                                       <  8:394637/136‑1 (MQ=255)
                   tcAAGGCCGAGCACGGCACCGGGCGCAACATGGCGCCCGTCGTGGAGCTGGAGTGGGGGGACGCCGCCTACCCGCTGATGTGGAAGCTCAAGCGCCTGCTCGCCCCCAACGGCATCCTCCACCCCGCCggggggg                                                                                                                     >  5:244468/1‑132 (MQ=255)
                     aaGGCCGAGCACGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGa                                                                                                                  <  8:131993/136‑1 (MQ=255)
                     aaGGCCGAGCACGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGa                                                                                                                  <  8:138083/136‑1 (MQ=255)
                     aaGGCCGAGCACGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGa                                                                                                                  <  7:94336/136‑1 (MQ=255)
                             gCACGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGa                                                                                                           <  4:2122/135‑1 (MQ=255)
                               aCGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCCGCGGATGTGGAAGCTCCAGCGCCTGCTCGACCCCCACGGCATCCTCAACCCCGCCGGGGTGTTGAGCGAAGccc                                                                                                         >  6:129132/1‑135 (MQ=255)
                               aCGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGCAGCTCCAGCCCCTGCGCGACCCCAACCGCATCCTCAACCCCGCCGTGGTGTGGCGCGAAGAcc                                                                                                         >  1:233346/1‑135 (MQ=255)
                               aCGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATTCTCaa                                                                                                                                   <  1:59305/109‑1 (MQ=255)
                               aCGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATTCTCaa                                                                                                                                   <  3:260970/109‑1 (MQ=255)
                               aCGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATTCTCaa                                                                                                                                   >  2:59305/1‑109 (MQ=255)
                               aCGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATTCTCaa                                                                                                                                   >  4:260970/1‑109 (MQ=255)
                                      cGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGatat                                                                                                  >  4:47229/1‑135 (MQ=255)
                                           gcAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATACAcc                                                                                             >  5:368957/1‑135 (MQ=255)
                                           gcAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCCACCGGGTGTTGAGCGAAGACCCGGATATCCACCt                                                                                            >  8:3154/1‑136 (MQ=255)
                                           gcAACATGGCGCCCGTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCCACCCCCACCGGGGGTTGAGCGAAGCACCGGATATCCa                                                                                               >  4:87323/1‑133 (MQ=255)
                                                aTGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGCCCTGGTGTTGAGCGAAGACCCGGATATCCACCTGaa                                                                                         >  2:193061/1‑134 (MQ=255)
                                                     gcCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACTGCATCCTCAACCCCGACCTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAAcc                                                                                    >  1:244576/1‑134 (MQ=255)
                                                     gcCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAAcc                                                                                    >  3:293981/1‑134 (MQ=255)
                                                      cccTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGACATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCt                                                                                   >  3:193966/1‑134 (MQ=255)
                                                       ccTTCTTGGAGCTGGAGTCGGGCCACGACGCCTACCAGCTCATGTGGAATCTCAAGCGCCTGCTCGCCCCCAAAGCCATCCTCAACCCCGACGTGGTGTTGACCGAATCCCCCCATAACCACCTGAAAAACCTCa                                                                                 >  2:345250/1‑135 (MQ=255)
                                                       ccTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCCCCTGAAAAACCt                                                                                   >  6:334537/1‑133 (MQ=255)
                                                                             gCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACa                                                            >  5:386175/1‑134 (MQ=255)
                                                                             gCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACa                                                            >  2:73162/1‑134 (MQ=255)
                                                                               actacgCCAACCAGCTGATGTGCAAGCTAAAGCGCCTGCTCGACCCCAACGGCAGCCTCACCCCCGACGTGGTTTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACaaaa                                                         <  4:293981/132‑1 (MQ=255)
                                                                               acgacgCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACaaaa                                                         >  6:367570/1‑135 (MQ=255)
                                                                                 gacgCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAAtc                                                       <  2:8105/135‑1 (MQ=255)
                                                                                  acgCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAAtcgt                                                     >  7:138069/1‑136 (MQ=255)
                                                                                         ccAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACaa                                               >  6:190059/1‑135 (MQ=255)
                                                                                           cgcggATGTGGAAGCGCACGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGCATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGt                                             <  1:207432/131‑1 (MQ=255)
                                                                                           aGCTGATGTGGAAGCTCAAGCGCCTGCTCGCCCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTg                                            >  8:30633/1‑136 (MQ=255)
                                                                                                aTGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAAcccc                                                                                                                               <  2:357508/48‑1 (MQ=255)
                                                                                                aTGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAAcccc                                                                                                                               >  1:357508/1‑48 (MQ=255)
                                                                                                aTGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCAt                                         >  4:309462/1‑134 (MQ=255)
                                                                                                 tgtgGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCAt                                         >  4:77830/1‑133 (MQ=255)
                                                                                                   tgGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGa                                                                                                                  >  1:118878/1‑58 (MQ=255)
                                                                                                   tgGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGa                                                                                                                  <  2:118878/58‑1 (MQ=255)
                                                                                                   tgGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAg                                     >  3:188082/1‑135 (MQ=255)
                                                                                                   tgGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAg                                     >  3:128118/1‑135 (MQ=255)
                                                                                                   tgGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGt                                    >  7:20197/1‑136 (MQ=255)
                                                                                                   tgGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAACTGCATCGAg                                     >  6:331018/1‑135 (MQ=255)
                                                                                                          tCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCg                                                                                                                              <  2:302732/39‑1 (MQ=255)
                                                                                                          tCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCg                                                                                                                              >  1:302732/1‑39 (MQ=255)
                                                                                                                gcCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGc                         >  1:134862/1‑134 (MQ=255)
                                                                                                                      tCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCg                                       <  6:225580/114‑1 (MQ=255)
                                                                                                                      tCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCg                                       >  5:225580/1‑114 (MQ=255)
                                                                                                                      tCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCgg                  >  4:31495/1‑135 (MQ=255)
                                                                                                                        taCCCCAACTGAATCCTCAACCCGGACTTTGTGTTGAGCGAAGACCCGTATATCCAACTTAAAAACGTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGAGAACCGgtg                <  2:244576/134‑1 (MQ=255)
                                                                                                                           cccAACGGCAGCCTCACCCCCGACGGGGTGTGGAGCGAAGCCCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCGGTGTGc             <  5:391611/135‑1 (MQ=255)
                                                                                                                             cAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCGGTGTGccc           >  6:322632/1‑135 (MQ=255)
                                                                                                                              aaCGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCGGTGTGCCCGt         <  8:138069/136‑1 (MQ=255)
                                                                                                                              aaCGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCGGTGTGCCCGt         <  8:62509/136‑1 (MQ=255)
                                                                                                                               aCGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCGGTGTGCCCGt         <  3:309462/135‑1 (MQ=255)
                                                                                                                               aCGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCGGTGTGCCCGt         <  3:129868/135‑1 (MQ=255)
                                                                                                                               aCGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCGGTGTGCCCGt         <  2:134862/135‑1 (MQ=255)
                                                                                                                               aCGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCGGTGTGCCCGt         <  1:6625/135‑1 (MQ=255)
                                                                                                                                 ggCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCGGTGTGCCCGTc        >  4:255210/1‑134 (MQ=255)
                                                                                                                                     tCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCGGTGTGCCCGTCCAAggg  <  7:135574/136‑1 (MQ=255)
                                                                                                                                      ccTCACCCCCGACGTGGTGTTGAGCGAAGACCCGGATAGCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCGGTGTGCCCGTCCAAggg  <  5:129132/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GTGGAGTTTGGCGGCTCGCTCAAGGCCGAGCACGGCACCGGGCGCAACATGGCGCCCTTCGTGGAGCTGGAGTGGGGGCACGACGCCTACCAGCTGATGTGGAAGCTCAAGCGCCTGCTCGACCCCAACGGCATCCTCAACCCCGACGTGGTGTTGAGCGAAGACCCGGATATCCACCTGAAAAACCTCAAGCCGCTGCCTGCCGCCGACAAAATCGTCGACAAGTGCATCGAGTGCGGCTTCTGCGAACCGGTGTGCCCGTCCAAGGG  >  ACB122/5381851‑5382119

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