Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,019,715 G→C 9.9% D113E (GAC→GAG)  ompA ← outer membrane porin A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,019,7150GC 9.9% 177.9 / 4.4 58D113E (GAC→GAG) ompAouter membrane porin A
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base C (4/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTG  >  NC_000913/1019570‑1019861
                                                                                                                                                 |                                                                                                                                                  
aTTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCa                                                                                                                                                <  1:1317237/150‑1 (MQ=255)
   ccAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGgt                                                                                                                                             >  2:979046/1‑150 (MQ=255)
   ccAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGgt                                                                                                                                             <  1:1360407/150‑1 (MQ=255)
   ccAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGgt                                                                                                                                             >  1:3570/1‑150 (MQ=255)
       aCGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGgtcgtc                                                                                                                                         <  2:3570/150‑1 (MQ=255)
       aCGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGgtcgtc                                                                                                                                         <  2:2313/150‑1 (MQ=255)
        cGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCa                                                                                                                                        >  2:1624347/1‑150 (MQ=255)
           gTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTg                                                                                                                                     >  2:208567/1‑150 (MQ=255)
               cGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTg                                                                                                                                 <  1:114466/150‑1 (MQ=255)
                 ttttCAGGAGTGATGGCGTACTCAACACCGCCAGCGAAGACCGTAGAAACGCCGGTGTCGGGGTTTTTACCATAAACGTTGTATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTggg                                                                                                                               <  2:697338/149‑1 (MQ=255)
                 aTTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTggg                                                                                                                               <  2:1228293/150‑1 (MQ=255)
                  tttCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGt                                                                                                                              <  1:979046/150‑1 (MQ=255)
                            gATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAAc                                                                                                                           >  1:1178028/1‑143 (MQ=255)
                            gATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAAc                                                                                                                           <  2:1178028/143‑1 (MQ=255)
                                gcgTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCg                                                                                                                >  1:107718/1‑150 (MQ=255)
                                     cTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAg                                                                                                           <  1:925649/150‑1 (MQ=255)
                                       cAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGtt                                                                                                         <  2:1620985/150‑1 (MQ=255)
                                         acacCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGa                                                                                                       >  1:72225/1‑150 (MQ=255)
                                           acCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAAc                                                                                                     <  2:270406/150‑1 (MQ=255)
                                            ccGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACg                                                                                                    >  1:532746/1‑150 (MQ=255)
                                             cGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATCTCCATGTACTCAGTGAAATTATAAAAAAGTTTAGCGGTCAGTTGAACGc                                                                                                   >  2:394319/1‑150 (MQ=255)
                                              gCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACgac                                                                                                  >  1:72994/1‑148 (MQ=255)
                                              gCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGcc                                                                                                  >  1:1304418/1‑150 (MQ=255)
                                                 aGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTg                                                                                               >  1:1135597/1‑150 (MQ=255)
                                                  gCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGa                                                                                              >  1:48819/1‑150 (MQ=255)
                                                   cGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAg                                                                                             <  2:675645/150‑1 (MQ=255)
                                                      aGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGGGTAGATGTCCAGGTCGGCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCtt                                                                                          <  2:1225580/150‑1 (MQ=255)
                                                          cGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTa                                                                                      <  2:1135597/150‑1 (MQ=255)
                                                          cGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTa                                                                                      <  1:1593913/150‑1 (MQ=255)
                                                          cGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTa                                                                                      <  1:208567/150‑1 (MQ=255)
                                                           ggAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTAt                                                                                     >  1:943475/1‑150 (MQ=255)
                                                               aaaCGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCAc                                                                                 <  1:177453/150‑1 (MQ=255)
                                                                 aCGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATCTCCATGTACTCAGTAAAATTAAAAAAAAATTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCg                                                                               >  2:738652/1‑150 (MQ=255)
                                                                                tttttACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCt                                                                >  1:1220765/1‑150 (MQ=255)
                                                                                 ttttACCATAAACGTTGGATTTAGTCTCTGCACGCCATACCATGCCACCCAGACGAGAGTAGATGTCCAGGTCGTCAGAGAGAGGGAAACCCATATTATCGGTAAGTAGTACGCCCAGAGCAATGTATGCACCGTTTTCAACGCtcgctt                                                               >  2:1619072/1‑145 (MQ=255)
                                                                                   ttACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTg                                                             <  1:816380/150‑1 (MQ=255)
                                                                                      ccATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTTGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTAc                                                          >  1:371478/1‑150 (MQ=255)
                                                                                      ccATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATCTCCATGTACTCAGTAAAATTAAAAAAAAATTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTAc                                                          >  2:1296777/1‑150 (MQ=255)
                                                                                       cATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACg                                                         <  2:107718/150‑1 (MQ=255)
                                                                                           aaCGTTGGATTTAGTGTCTGCACGCCATTCCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCAt                                                     <  2:40149/150‑1 (MQ=255)
                                                                                                   aTTTAGTGTCTGCACGCCATACCATGCGACCGAGACGAGTGTAGATGTCCAGGTCGTCAGTGATGGGGTACCCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTa                                             <  2:1526301/150‑1 (MQ=255)
                                                                                                          gtCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTACTCAGTAATTTGTTAACCAAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTc                                      >  2:158603/1‑150 (MQ=255)
                                                                                                           tCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCg                                     <  2:532746/150‑1 (MQ=255)
                                                                                                              gCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTaa                                  >  1:94128/1‑150 (MQ=255)
                                                                                                                 cGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAAccc                               >  1:163888/1‑150 (MQ=255)
                                                                                                                      tACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAAcc                                >  1:680007/1‑144 (MQ=255)
                                                                                                                      tACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAAcc                                <  2:680007/144‑1 (MQ=255)
                                                                                                                      tACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTc                          <  2:1220765/150‑1 (MQ=255)
                                                                                                                      tACCATGCCACCCAGACGAGTGTAGATCTCCATGTACTCAGTAAAATTAAAAAAAAAATTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTc                          >  2:109894/1‑150 (MQ=255)
                                                                                                                        ccATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCaa                        >  1:425103/1‑150 (MQ=255)
                                                                                                                          ctGCCACCCAGACGCGCGGCGCTCTCCAGGTCGTCAGGGCTTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAg                      <  2:1582468/149‑1 (MQ=255)
                                                                                                                          aTGCCCCCCCGACGCGTGTCGATCTCCAGGTCGTCCGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAg                      <  2:1576783/150‑1 (MQ=255)
                                                                                                                            gCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGcc                    <  2:1304418/150‑1 (MQ=255)
                                                                                                                                     aCGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACgg           <  2:163888/150‑1 (MQ=255)
                                                                                                                                        aGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGa                                                 <  2:1155183/109‑1 (MQ=255)
                                                                                                                                        aGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGa                                                 >  1:1155183/1‑109 (MQ=255)
                                                                                                                                            tAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTAAGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACACATTTCAAAGCCAACATACGGGTTAAcc    >  2:626490/1‑150 (MQ=255)
                                                                                                                                              gATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACAGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTg  >  1:610282/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATTCCAGACGGGTAGCGATTTCAGGAGTGATCGCGTACTCAACACCGCCAGCGAAGACCGGAGAAACGCCGGTGTCGTGGTTTTTACCATAAACGTTGGATTTAGTGTCTGCACGCCATACCATGCCACCCAGACGAGTGTAGATGTCCAGGTCGTCAGTGATTGGGTAACCCAGTTTAGCGGTCAGTTGAACGCCCTGAGCTTTGTATGCACCGTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTG  >  NC_000913/1019570‑1019861

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: ‑