Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,305,853 C→T 35.8% D126N (GAC→AAC)  ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,305,8530CT35.8% 45.7 / 45.0 57D126N (GAC→AAC) ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (16/20);  new base T (4/16);  total (21/36)
Fisher's exact test for biased strand distribution p-value = 8.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01

CGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCC  >  CP009273/2305706‑2305999
                                                                                                                                                    |                                                                                                                                                  
cGTTTTTCCCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCg                                                                                                                                                   <  1:460992/149‑1 (MQ=255)
cGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCg                                                                                                                                                   <  1:650833/149‑1 (MQ=255)
cGTTTTGACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCg                                                                                                                                                   <  1:568436/149‑1 (MQ=255)
  tttttACCCTGGGACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTTGGt                                                                                                                                                 <  2:315792/149‑1 (MQ=255)
      tACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTACAg                                                                                                                                             >  2:936709/1‑149 (MQ=255)
         ccTGGTACTGAACAGCAAAGGTCAGGCCGTCAACCAGACCGAAGAAGTCAGCGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACCTCGGTCCACGaa                                                                                                                                          <  2:607763/149‑1 (MQ=255)
           tGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGCAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTAACCACCGAATTCTGGCAGTACGTCGTTCAAGGAAGt                                                                                                                                        >  2:276250/1‑149 (MQ=255)
              tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTAc                                                                                                                                     <  2:483496/149‑1 (MQ=255)
              tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTAc                                                                                                                                     <  1:214808/149‑1 (MQ=255)
               aCTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTTGGTCCAGGAAGTTACg                                                                                                                                    <  1:781767/149‑1 (MQ=255)
                  gAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTTGGTCCAGGAAGTTACGTCa                                                                                                                                 <  1:534368/149‑1 (MQ=255)
                   aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCAt                                                                                                                                <  2:250065/149‑1 (MQ=255)
                       gCAAAGTTCAGGCCGTCAACCAGAACGAAGAAGTCAGTGTAACGTTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATGCTGGCAGTACGTCGGTCCAGGAAGTTACGTCAAAaac                                                                                                                            >  2:519638/1‑149 (MQ=255)
                         aaaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTTGGTCCAGGAAGTTACGTCATAaacaa                                                                                                                          <  1:752134/149‑1 (MQ=255)
                          aggTTCAGGCCGGCAACCAGACCGAAGAAGTCAGTGTTACGGTCGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGGAGTTGTCAGAACCGTAGGTGTCACCTCCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaacaac                                                                                                                         <  2:205801/147‑1 (MQ=255)
                            gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCAAAT‑GAATGAAGTTATAAAAACAGTCGGTGACACCACAGAAAACTGCCAGAACGAGAGTCCACGAACTACCGCCATAAACAACGc                                                                                                                       >  2:854513/1‑149 (MQ=255)
                             ttCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGGGACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCACGAAGTTACGTCATAAACAACGcc                                                                                                                      <  1:653049/149‑1 (MQ=255)
                               cAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGt                                                                                                                    <  2:714227/149‑1 (MQ=255)
                                  gCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCCCGCCCCT‑GCCTGAAGTTGTCATAACCGTAGGTGTCACCACCGAATTCTGCCAGTAAGTCGGTCCCGGAAGTGACGTCATAAACAAGGCCGTAGt                                                                                                                 >  2:1123443/1‑149 (MQ=255)
                                     gTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGTTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTAc                                                                                                              <  2:440332/149‑1 (MQ=255)
                                       cAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGa                                                                                                            <  1:145091/149‑1 (MQ=255)
                                       cAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGa                                                                                                            <  1:253034/149‑1 (MQ=255)
                                       acgccAGCCCGAAGAAGTCAGGGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGa                                                                                                            <  2:538978/146‑1 (MQ=255)
                                                 gaagaaGTCAGTGTTACGGACGGGCGCGAAGCCGTTACAAC‑CAGTTGGCATGACATTGTCAACAACAAATTAGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                                                                                                  >  2:875822/1‑149 (MQ=255)
                                                   agaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGaaa                                                                                                <  1:332519/149‑1 (MQ=255)
                                                   agaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGaaa                                                                                                <  1:955276/149‑1 (MQ=255)
                                                    gaagGCAGTGTTACGGGACGTCGCGAAGCCCTTACCAGGCTGGC‑GCATGGAGGTGTCAGAACCGTAGGTGTCACCACCGAATTCGGGCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                                                                                               <  2:116146/149‑1 (MQ=255)
                                                    gaagCCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGACAAAGTAACCACCGAATTTTTTCATTACGTTGGTCCAGTAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                                                                                               <  1:600394/149‑1 (MQ=255)
                                                              ttACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAGCCCACATc                                                                                     >  2:862390/1‑149 (MQ=255)
                                                               tagggcaGGTCGCGAAGCCGTAGCCACGTGGCT‑GCATGAATTGGTCAGAACCGTAGGTGTCACCACCGAATTCTGTCAGTACGCCGGTCCAGGAAGGTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGACCCcacgttc                                                                                    <  2:325119/143‑5 (MQ=255)
                                                                       gTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGCAttt                                                                            >  2:1094388/1‑149 (MQ=255)
                                                                         cgcgAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCa                                                                          <  1:845174/149‑1 (MQ=255)
                                                                                   ttACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGaa                                                                >  1:151345/1‑149 (MQ=255)
                                                                                    tACCACGCGGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGCATTTCAGACCTGCGAAt                                                               <  2:1226385/149‑1 (MQ=255)
                                                                                      ccACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTCTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACTTCATAAACAACGGCGGAGTTACGACCGTAGTCGAAAGAACCCACACCCTGGAATTTCAGACCTGCGAATGc                                                             <  1:176158/149‑1 (MQ=255)
                                                                                      ccACGCTGCG‑GCATGAAGTTGTCAGCACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAACCAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGc                                                             <  2:646944/149‑1 (MQ=255)
                                                                                        aCGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACAGCCTGGAATTTCAGACCTGCGACTGCca                                                           >  2:1022222/1‑149 (MQ=255)
                                                                                          gctgct‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCcaca                                                         >  1:800703/1‑149 (MQ=255)
                                                                                          gctgct‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCcaca                                                         >  1:76633/1‑149 (MQ=255)
                                                                                                  cATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGGCCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTcc                                                  >  2:229165/1‑149 (MQ=255)
                                                                                                  cATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTcc                                                  <  1:399603/149‑1 (MQ=255)
                                                                                                              aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTgtttt                                      >  1:1198365/1‑149 (MQ=255)
                                                                                                               gAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTgttttc                                     >  1:224986/1‑149 (MQ=255)
                                                                                                                  ccGTAGGGTTCACCACCGAATCCTGGCAGTCCGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTgttttcgtt                                  <  2:278276/149‑1 (MQ=255)
                                                                                                                    gTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAAAGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGATTGCCACACGGGTCTAGGAGTTgttttcgtttt                                >  1:812399/1‑149 (MQ=255)
                                                                                                                     tAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCCGGCGTTgttttcgttttc                               >  1:997737/1‑149 (MQ=255)
                                                                                                                     tAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGACgctttcgttttc                               >  2:1091252/1‑149 (MQ=255)
                                                                                                                        gtgtCACCACCGAATTCTGGCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAgc                            >  2:336904/1‑149 (MQ=255)
                                                                                                                         tctCACCACCGAATTCTGGCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAgcg                           <  1:441145/147‑1 (MQ=255)
                                                                                                                              ccaccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGtt                      >  1:166928/1‑149 (MQ=255)
                                                                                                                                accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTCCAGCGCTGTTGc                    >  2:529130/1‑149 (MQ=255)
                                                                                                                                        cTGGCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATc            <  1:569695/149‑1 (MQ=255)
                                                                                                                                         tGGCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACCCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCt           <  1:1029741/149‑1 (MQ=255)
                                                                                                                                           ggaGTACGTCGGTGCAGGAAGTTACGTCATAAACAAACGCGTAGCCACGACCGGAGACGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAAGGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGAGCTGTTGCCCTGGATCTGa         <  2:566814/147‑1 (MQ=255)
                                                                                                                                           gcagcaCGTTGGGCCAGGAAGTTACGTCATAAACCACGCCGTAGTGACGACGCTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCGGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGa         <  2:335928/149‑1 (MQ=255)
                                                                                                                                           gCAGTACGTTGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGa         <  1:23014/149‑1 (MQ=255)
                                                                                                                                                  gTCGGTCCAGGAAGTTAGGTCATAAACAACGCCGTAGTTATGACCGTAGTCGAAAGAACTAACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTccc  <  1:613161/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCT‑GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCC  >  CP009273/2305706‑2305999

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑