Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 982,560 T→A intergenic (‑122/+481) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273982,5600TA100.0% 183.7 / NA 59intergenic (‑122/+481)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (31/28);  total (31/28)

CAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACA  >  CP009273/982418‑982705
                                                                                                                                              |                                                                                                                                                 
cAGGATATTCCGCTTCATCATTATTTATTACCCTCCTGGTTTTTTTTACGCCACCTGCCACTGCCGTCAATCAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTCGTGTCTTTATGTGTCTGCAGGCATCTTACCATTc                                                                                                                                             <  2:19809/149‑1 (MQ=255)
    atatTGCGCTTCATCATTATTTATTTCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCCGGCATCTTACCATTCAAAc                                                                                                                                         <  1:894146/149‑1 (MQ=255)
        tGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAAc                                                                                                                                     <  2:43136/149‑1 (MQ=255)
        tGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAAc                                                                                                                                     >  1:389525/1‑149 (MQ=255)
                 tcatTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACAATTAAAACTAACGTTTCGCTa                                                                                                                            >  1:1193990/1‑149 (MQ=255)
                   attattTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTAcc                                                                                                                          >  1:1015426/1‑149 (MQ=255)
                     tatttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTTTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGt                                                                                                                        >  1:1078656/1‑149 (MQ=255)
                         tattaCCCTCATGGTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAg                                                                                                                   <  1:962131/149‑1 (MQ=255)
                          attaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAg                                                                                                                   >  1:582090/1‑149 (MQ=255)
                           ttaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAAAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGt                                                                                                                  >  1:69811/1‑149 (MQ=255)
                             aCCCGCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGc                                                                                                                <  1:1031145/149‑1 (MQ=255)
                               ccTCGTGGTTTTTTTTATGACACCTGCGACTGCCGTCAATAAGTTCTTTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTa                                                                                                              <  2:307098/149‑1 (MQ=255)
                                 tCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACa                                                                                                            >  1:1103199/1‑149 (MQ=255)
                                         ttttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaa                                                                                                    >  1:608858/1‑149 (MQ=255)
                                         ttttttATGACACCTGCCACGGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagataaa                                                                                                    <  2:851636/149‑1 (MQ=255)
                                              tATGACACCTTCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAttt                                                                                               <  1:221315/149‑1 (MQ=255)
                                                tGACACCTTCCACTGCCGGCAATAAGTTCTTTCAATAAAAATTTACGTAACTATTGATGAGAGTTTGGTGTCTTTATGTTTCTGCAGGCATCTTACCATTCAGACTAACGTTTCGCTTCCGTGAAAGTGCTACAAATATAAAGATTTgg                                                                                             <  2:88764/149‑1 (MQ=255)
                                                  acacTTGTCACTGCCGTCAATAAGTTCTGTCAATTTAAATTTCCGGTATTCTTAGACATAATTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGtt                                                                                           <  1:540366/149‑1 (MQ=255)
                                                  acacCTGCCGCTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAGTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGtt                                                                                           <  2:152754/149‑1 (MQ=255)
                                                  acacCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGtt                                                                                           >  1:957292/1‑149 (MQ=255)
                                                      ccgCCACTGTCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaa                                                                                       <  2:1130202/147‑1 (MQ=255)
                                                       tGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAAAGTTTGGTGTCATTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaa                                                                                      >  1:39870/1‑149 (MQ=255)
                                                         ccACTGCCGTCAATAAGTTTTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaaa                                                                                    <  1:260962/149‑1 (MQ=255)
                                                           aCTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGa                                                                                  >  1:753860/1‑149 (MQ=255)
                                                                  acAATAAGTTCTGTCAATAAAAATTTACGGACCTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGATACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAata                                                                           >  1:420763/2‑149 (MQ=255)
                                                                   cAATAAGTTCTGTCAATTAAAATTTACGGATCTATTGATGAGGGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAatat                                                                          <  1:525867/149‑1 (MQ=255)
                                                                      tcagatCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATAAGTa                                                                       >  2:1128775/6‑149 (MQ=255)
                                                                       aaGTTCTGTCAATAAAAATTTACGTAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAAAGTTTCGATACCGTGAAAGTGATACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTaa                                                                      <  1:146003/149‑1 (MQ=255)
                                                                            ctctccaTAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCaaa                                                                 <  2:430884/143‑1 (MQ=255)
                                                                                 aaTAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTaa                                                            >  2:1224421/1‑149 (MQ=255)
                                                                                  aTAAAAATTTACGGAACTATTGATGAGAGTTTGGAGATTTTATGTGTCTGCAGGCAACTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATATAGTTTCAAAAAGAAAAAATATGTAACCAAAAGTaaa                                                           <  2:207991/149‑1 (MQ=255)
                                                                                    aaaaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAAt                                                         >  2:171113/1‑149 (MQ=255)
                                                                                      aaaaTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAAttt                                                       <  1:493955/149‑1 (MQ=255)
                                                                                            aCGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGaa                                                 >  1:595115/1‑149 (MQ=255)
                                                                                              ggAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAAATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAAGAACt                                               >  2:1046952/1‑149 (MQ=255)
                                                                                                 acgaTTGACTAGAATTTGGTGTCTTTATCTTTCTGCAGGCATCGTACCCTTCAAACTACCGTGTCGCTACCGTGAAAGTGCTACAAACATAAAGATTTGGAAAAAAAAAAAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtg                                            <  2:650461/146‑1 (MQ=255)
                                                                                                  cTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTAAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtgt                                           <  2:1036167/149‑1 (MQ=255)
                                                                                                   tATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAACTATGTAACCATAAGTAAAATTTAAGGAAcgtggtc                                          >  2:181619/1‑143 (MQ=255)
                                                                                                   tATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtgtg                                          >  2:187755/1‑149 (MQ=255)
                                                                                                   tATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtgtg                                          >  2:471840/1‑149 (MQ=255)
                                                                                                    aTTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGa                                         <  2:582090/149‑1 (MQ=255)
                                                                                                        atgaGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAAAATGTAACCAAAAGTAAAATTTAAGGAACTATGTGAacac                                     >  2:269480/1‑149 (MQ=255)
                                                                                                        atgaGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCACTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAacac                                     >  1:32805/1‑149 (MQ=255)
                                                                                                          gagaGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCg                                   <  1:1188398/149‑1 (MQ=255)
                                                                                                          gagaGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCg                                   >  2:396961/1‑149 (MQ=255)
                                                                                                            gagTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGTTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTc                                 <  1:314163/149‑1 (MQ=255)
                                                                                                            gagTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAATTAACTTTGTGACCACCGTc                                 <  2:576491/149‑1 (MQ=255)
                                                                                                                ttgcgctttttATGTAGCTGCAGGCATATAACCAGTCGATCTAACGTTGCGCTACCGTGAAAGTGCCACAAAGAGAAAGATTTTGTTTCAAAAAAAAAAACGATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCatat                             <  2:239246/141‑1 (MQ=255)
                                                                                                                 tGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGATTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATATAAGGAAATTTGTGAACACCGTCATAtt                            >  1:1078902/1‑149 (MQ=255)
                                                                                                                    tgtATTTATGTGTCTGAAGGCATCTTAAAATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTCTGTGAACACCGTCATATATcc                         >  2:957022/1‑149 (MQ=255)
                                                                                                                      tCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCAt                       <  2:957292/149‑1 (MQ=255)
                                                                                                                        tttATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATag                     >  2:414684/1‑149 (MQ=255)
                                                                                                                        tttATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGGGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCAAag                     >  1:654731/1‑149 (MQ=255)
                                                                                                                            tgtgTCTACACGCATCGTACCATTCAAACTACCCCTTCGCTAACTTGAAAATGCTACAAAGATAAAGATTTAGAAAAAAAAAAAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGAc                 <  2:1193990/149‑1 (MQ=255)
                                                                                                                             gtgtCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACg                >  1:1000653/1‑149 (MQ=255)
                                                                                                                              tgtCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGt               <  1:471840/149‑1 (MQ=255)
                                                                                                                                     aGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTAATATTTCCATAGAGACGTGATGata        >  2:668155/1‑149 (MQ=255)
                                                                                                                                     aGGCATCTTACCATTCAAACTAACGTTTCGCTACCATGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGGACGGTGTGAACACCGTCATATTTCCATAGAGACGTGATGata        >  2:1175815/1‑149 (MQ=255)
                                                                                                                                           cTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTATGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACa  >  1:65852/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
CAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACA  >  CP009273/982418‑982705

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