Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,511,541 (G)5→6 coding (874/957 nt) fecD ← ferric citrate ABC transporter membrane subunit FecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,511,5361.G94.2% 164.4 / 2.7 52coding (879/957 nt)fecDferric citrate ABC transporter membrane subunit FecD
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base G (20/29);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 8.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGA  >  NC_000913/4511400‑4511665
                                                                                                                                          |                                                                                                                                    
cgtTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTTAGCACGCCAAACGGGGGGTCCCGTggg                                                                                                                                 >  2:115563/1‑140 (MQ=255)
agacaggTACCTTGTCTTTCCCCTAACTGACTGTCCGATATTCAGTTCCTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTggg                                                                                                                                 <  2:230534/134‑1 (MQ=255)
   tAAGTACCTTGTCTGTCCCGTAACTTACCTTCAGATTGTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTgggggg                                                                                                                             <  1:300335/140‑2 (MQ=255)
   tAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTgggggg                                                                                                                             <  2:147129/140‑2 (MQ=255)
            tgtctgtcCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGaataat                                                                                                                     >  2:843729/1‑140 (MQ=255)
                 gtcCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGTCAACCGGGAGCGCCGGTGGGGGGATGaataactgtct                                                                                                                >  1:230534/1‑134 (MQ=255)
                 gtcCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCg                                                                                                                >  2:584961/1‑140 (MQ=255)
                    ccGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGGGGGGGGATGAAAAATTCTCGCca                                                                                                             >  1:856165/1‑140 (MQ=255)
                     cGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCCCCAACGCGGGGATCAATAATTCTCGCcag                                                                                                            >  1:419638/1‑140 (MQ=255)
                       tAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCCGGAGCTCCAGTGGGGGGATGAATAATTCTCGCcagca                                                                                                          >  1:371734/1‑140 (MQ=255)
                              ccGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTTCAGTGGGTGGAGGAATAATTCACGCCAGCAGATCGGc                                                                                                   >  1:265319/1‑140 (MQ=255)
                                gTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCaa                                                                                                 <  1:209613/140‑1 (MQ=255)
                                gTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCaa                                                                                                 <  2:203374/140‑1 (MQ=255)
                                gTCAGATTTTCAGTTCGGAAAGTCATTTATCCCATCCTCTCAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCaa                                                                                                 <  2:371734/140‑1 (MQ=255)
                                     aTTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAaccacc                                                                                            >  1:125431/1‑140 (MQ=255)
                                       tttCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCCCCAGTGGGGGGGTGAATAATTCTCGCCAGCAGATCGGCAACCACCaa                                                                                          >  1:523684/1‑140 (MQ=255)
                                               cGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGGCAGAAGCTCGGCAACCACCAACAGCAACg                                                                                  >  1:333674/1‑140 (MQ=255)
                                               cGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACg                                                                                  >  2:46637/1‑140 (MQ=255)
                                               cGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACg                                                                                  >  1:459792/1‑140 (MQ=255)
                                                 tAAAGTCATTTATCGTATTCTCACATGCAGCCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCa                                                                                <  1:371085/140‑1 (MQ=255)
                                                 tAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCa                                                                                >  1:720820/1‑140 (MQ=255)
                                                       cATTTATCGCATTCTCACAAGCTACCAGACATACCTCGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTc                                                                          <  2:333674/140‑1 (MQ=255)
                                                              cGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTTTCGCCAGCAGATTGGCTACCACCAAGAGCAACGCACCTGTCAGGGCCg                                                                   >  2:283118/1‑140 (MQ=255)
                                                                cATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCCCCAACAACAACGCCCCTGTCAGGGCCGaa                                                                 >  1:336813/1‑140 (MQ=255)
                                                                cATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACg                                                                                  >  2:677481/1‑123 (MQ=255)
                                                                cATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACg                                                                                  <  1:677481/123‑1 (MQ=255)
                                                                cATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGaa                                                                 >  1:38261/1‑140 (MQ=255)
                                                                       acaAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGgca                                                          <  1:46637/140‑1 (MQ=255)
                                                                        caAGCAACCAGACAAACGACGGGGCACCGATTATGGCGGGCAGCACGCTAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGgcag                                                         <  2:336813/140‑1 (MQ=255)
                                                                        caAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGgcag                                                         <  2:38261/140‑1 (MQ=255)
                                                                            cAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGt                                                     <  1:845884/140‑1 (MQ=255)
                                                                             aaCCAGACAAACCACGGTGCACCGATAATGGCGGTCAGCTCGCCAACCTGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTc                                                    <  2:265319/140‑1 (MQ=255)
                                                                               ccAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTg                                                  >  2:482058/1‑140 (MQ=255)
                                                                                 aGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCg                                                <  1:668756/140‑1 (MQ=255)
                                                                                     aaaCCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCGGGGACGGAACAGGCAGCAGTCTGCGGTGa                                            >  1:287833/1‑140 (MQ=255)
                                                                                       aCCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACg                                          <  1:482058/140‑1 (MQ=255)
                                                                                          aCGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTcc                                       <  2:459792/140‑1 (MQ=255)
                                                                                          aCGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTcc                                       <  2:125431/140‑1 (MQ=255)
                                                                                          aCGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTcc                                       <  2:856165/140‑1 (MQ=255)
                                                                                            ggCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCAc                                     <  1:584961/140‑1 (MQ=255)
                                                                                             gcgcACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTg                                             >  2:692490/1‑131 (MQ=255)
                                                                                             gcgcACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTg                                             <  1:692490/131‑1 (MQ=255)
                                                                                             gcgcACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCAcc                                    >  1:770514/1‑140 (MQ=255)
                                                                                               gcACCGATGATGGCGGCGAGCACGCCAACCAGGAGATCCGGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCgg                                  <  2:419638/140‑1 (MQ=255)
                                                                                               gcACCGATAATGGCGGTCAGCACGCCAACCGTGAGCTCCTGTGGGGGGGTGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCgg                                  <  2:310801/140‑1 (MQ=255)
                                                                                                 aCCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGACGAAACATGCATCAGTCTGCGATGACGTCCACCGGTg                                >  1:409022/1‑140 (MQ=255)
                                                                                                   cGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCCCCAGCAGGTCGGCAACCACCAACAGCAACGCACCTGTCAAGGCCGAAACAGGCAGCAGGCTGCGGGGACGCCCACCGGTgct                              >  2:757320/1‑138 (MQ=255)
                                                                                                     aTAATGGCGGTCAGCTCGCCAATCGTGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGc                            <  1:283118/140‑1 (MQ=255)
                                                                                                        aTGGCGGTTAGCACGCCAACCGTGAGCTCCCCTGGCGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTAc                         <  2:409022/140‑1 (MQ=255)
                                                                                                         tGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACg                        <  1:843729/140‑1 (MQ=255)
                                                                                                           gCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGGGACGTCCACCGGTGATGCTACGca                      >  1:840625/1‑140 (MQ=255)
                                                                                                             ggTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGt                                 <  1:597516/127‑1 (MQ=255)
                                                                                                             ggTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGt                                 >  2:597516/1‑127 (MQ=255)
                                                                                                             ggTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCt                           >  1:702018/1‑133 (MQ=255)
                                                                                                             ggTCAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCt                           <  2:702018/133‑1 (MQ=255)
                                                                                                              gTCAGCACGCCAACCGGGAGCTCCAGGGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGcatca                   >  2:97959/1‑140 (MQ=255)
                                                                                                                cAGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTAGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAGACAGGCAGCAGTGCGCGGTGACGTCGACCGGTGAGGCTACCCATCata                 >  2:319366/1‑140 (MQ=255)
                                                                                                                 aGCACGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCatat                <  2:770514/140‑1 (MQ=255)
                                                                                                                 aGCACGCCAACCGGGAGCTCCACTGGGGGGATGAAAAATTCTCGCCAGAAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCatat                <  2:287833/140‑1 (MQ=255)
                                                                                                                     cGCCAACCGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCgg            <  1:757320/140‑1 (MQ=255)
                                                                                                                           ccGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGTTGATGCTACGCATCATATGCGGcaccac      <  2:523684/140‑1 (MQ=255)
                                                                                                                           ccGGGAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGcaccac      <  1:115563/140‑1 (MQ=255)
                                                                                                                               gAGCTCCAGTGGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACgaga  <  2:840625/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                    
CGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGA  >  NC_000913/4511400‑4511665

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