Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,212,088 (C)8→7 intergenic (‑85/+172) iraM ← / ← ymgK anti‑adaptor protein IraM, inhibitor of sigma(S) proteolysis/protein YmgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,212,0810C.96.8% 252.2 / 2.6 62intergenic (‑78/+179)iraM/ymgKanti‑adaptor protein IraM, inhibitor of sigma(S) proteolysis/protein YmgK
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (39/21);  total (40/22)
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.78e-01

CTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGAT  >  NC_000913/1211946‑1212212
                                                                                                                                       |                                                                                                                                   
cTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccc                                                                                                                                 >  1:903913/1‑140 (MQ=255)
     aaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccccctca                                                                                                                            <  1:411594/140‑4 (MQ=255)
             aCCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAaaata                                                                                                                   >  1:576044/1‑140 (MQ=255)
                   tGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAaattaa                                                                                                             >  2:293719/1‑140 (MQ=255)
                        gTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTaaa                                                                                                        >  1:99132/1‑140 (MQ=255)
                           gAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                                                     <  2:576044/140‑1 (MQ=255)
                           gAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                                                     <  1:342270/140‑1 (MQ=255)
                            aaTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAatga                                                                                                    >  1:182065/1‑140 (MQ=255)
                            aaTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAatga                                                                                                    >  2:221391/1‑140 (MQ=255)
                            aaTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAatga                                                                                                    >  2:933622/1‑140 (MQ=255)
                            aaTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAatga                                                                                                    >  1:632353/1‑140 (MQ=255)
                              ttACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAatgatg                                                                                                  >  2:260586/1‑140 (MQ=255)
                                  cGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCtt                                                                                              >  2:194944/1‑140 (MQ=255)
                                     gtCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtata                                                                                           <  1:27218/140‑1 (MQ=255)
                                                tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATaaaataaaat                                                                                >  2:801912/1‑140 (MQ=255)
                                                tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATaaaataaaat                                                                                >  1:378318/1‑140 (MQ=255)
                                                      tCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAg                                                                          <  1:293719/140‑1 (MQ=255)
                                                       cATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGc                                                                         >  1:944873/1‑140 (MQ=255)
                                                        aTTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCa                                                                        <  1:221391/140‑1 (MQ=255)
                                                        aTTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCa                                                                        <  1:186407/140‑1 (MQ=255)
                                                          ttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAg                                                                      >  2:410211/1‑140 (MQ=255)
                                                              tATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAt                                                                  >  2:46838/1‑140 (MQ=255)
                                                              tATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAt                                                                  <  1:991110/140‑1 (MQ=255)
                                                               aTGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATGAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAtc                                                                 >  1:330377/1‑140 (MQ=255)
                                                                       tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGa                                                         >  1:964190/1‑140 (MQ=255)
                                                                       tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGa                                                         >  2:488291/1‑140 (MQ=255)
                                                                       tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGa                                                         >  2:646638/1‑140 (MQ=255)
                                                                       tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGa                                                         >  1:263783/1‑140 (MQ=255)
                                                                         ccACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATg                                                       <  1:315636/140‑1 (MQ=255)
                                                                            cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGtt                                                     >  2:755214/1‑140 (MQ=255)
                                                                            cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGa                                                                    <  2:354828/124‑1 (MQ=255)
                                                                            cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGa                                                                    >  1:354828/1‑124 (MQ=255)
                                                                            cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTa                                                    >  2:293229/1‑140 (MQ=255)
                                                                              gTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa                                                  <  2:964190/140‑1 (MQ=255)
                                                                                         ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt                                       >  1:755159/1‑140 (MQ=255)
                                                                                         ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt                                       >  2:22336/1‑140 (MQ=255)
                                                                                         ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt                                       >  1:188972/1‑140 (MQ=255)
                                                                                           ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAttt                                     >  2:360817/1‑140 (MQ=255)
                                                                                           ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAttt                                     >  1:920507/1‑140 (MQ=255)
                                                                                           ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAttt                                     >  2:697947/1‑140 (MQ=255)
                                                                                           ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAttt                                     >  1:945990/1‑140 (MQ=255)
                                                                                           ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAttt                                     >  1:395911/1‑140 (MQ=255)
                                                                                            tCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAtttt                                    >  1:317502/1‑140 (MQ=255)
                                                                                             cATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTtg                                   >  2:475107/1‑140 (MQ=255)
                                                                                               ttAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTtgtg                                 >  1:736982/1‑140 (MQ=255)
                                                                                                 aGGATGTTTATTTCTTGATTTTGCATATGAGTATATTATCCCCCCTCAAAAAAATAAATTAATTACAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGaa                               <  2:725533/140‑1 (MQ=255)
                                                                                                     tGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCa                           >  1:1013271/1‑140 (MQ=255)
                                                                                                     tGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCa                           >  2:543690/1‑140 (MQ=255)
                                                                                                       tttATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAAt                         >  2:478857/1‑140 (MQ=255)
                                                                                                       tttATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAAt                         <  1:194944/140‑1 (MQ=255)
                                                                                                       tttATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAAt                         >  2:137406/1‑140 (MQ=255)
                                                                                                       tttATTTCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAAt                         <  1:752364/140‑1 (MQ=255)
                                                                                                             tCTTGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAatta                   >  1:104434/1‑140 (MQ=255)
                                                                                                               ttGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAattatt                 >  1:367167/1‑140 (MQ=255)
                                                                                                                tGATTTTGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTg                >  1:360879/1‑140 (MQ=255)
                                                                                                                   ttttGCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACt             <  2:903913/140‑1 (MQ=255)
                                                                                                                     ttGCATATGAGTATATTACCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACtt            <  1:755214/140‑1 (MQ=255)
                                                                                                                       gCATATGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGt         <  2:360879/140‑1 (MQ=255)
                                                                                                                         atatGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGttt       <  2:755159/140‑1 (MQ=255)
                                                                                                                            tGAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAg    <  1:697947/140‑1 (MQ=255)
                                                                                                                             gAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGa   <  2:367167/140‑1 (MQ=255)
                                                                                                                             gAGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGa   <  1:260586/140‑1 (MQ=255)
                                                                                                                              aGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGAt  <  1:435973/140‑1 (MQ=255)
                                                                                                                              aGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGAt  <  1:661878/140‑1 (MQ=255)
                                                                                                                              aGTATATTACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGAt  <  1:360817/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
CTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGAT  >  NC_000913/1211946‑1212212

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: