Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,827,095 T→C 54.3% V1365V (GTA→GTG pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,827,0950TC54.3% 15.4 / 70.9 59V1365V (GTA→GTGpvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (15/12);  new base C (17/15);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCAGCGGCCAGAAGCATTCCCACACCGAGACGTCGAAACTGAT  >  NC_002947/4826949‑4827241
                                                                                                                                                    |                                                                                                                                                  
cGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTa                                                                                                                                                   >  5:38988/1‑148 (MQ=255)
      tgCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTg                                                                                                                                            <  3:220731/149‑1 (MQ=255)
        cAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTg                                                                                                                                            >  6:104884/1‑147 (MQ=255)
         aGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAg                                                                                                                                         >  5:83685/1‑149 (MQ=255)
          gTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGc                                                                                                                                        >  5:105993/1‑149 (MQ=255)
           ttGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGAACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGc                                                                                                                                        <  1:53000/149‑1 (MQ=255)
              cACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTa                                                                                                                                     <  5:104884/148‑1 (MQ=255)
              cACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTa                                                                                                                                     <  7:133738/148‑1 (MQ=255)
               aCCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAg                                                                                                                                         <  4:122051/143‑1 (MQ=255)
               aCCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAg                                                                                                                                         >  3:122051/1‑143 (MQ=255)
               aCCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGt                                                                                                                                   >  2:108441/1‑149 (MQ=255)
                         gcaCCTGCAGCACGCGGTCACGCAGCGTGTGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCAt                                                                                                                           >  1:27598/1‑148 (MQ=255)
                           aCCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGc                                                                                                                       <  6:38988/149‑1 (MQ=255)
                                              gCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGACCACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATa                                                                                                    <  1:108441/149‑1 (MQ=255)
                                                aGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACg                                                                                                  >  8:258795/1‑149 (MQ=255)
                                                    tGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGc                                                                                              <  2:204350/149‑1 (MQ=255)
                                                       gCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCg                                                                                             <  2:34937/148‑1 (MQ=255)
                                                         ggACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAgg                                                                                                             <  6:32194/130‑1 (MQ=255)
                                                         ggACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAgg                                                                                                             >  5:32194/1‑130 (MQ=255)
                                                         ggACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGTATCGCTGccctg                                                                                          >  3:183624/1‑145 (MQ=255)
                                                           aCAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTGCCCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTcg                                                                                        >  4:349890/1‑149 (MQ=255)
                                                            cAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTcgc                                                                                       >  6:260375/1‑149 (MQ=255)
                                                               cGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGAt                                                                                    >  1:229482/1‑149 (MQ=255)
                                                               cGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGAt                                                                                    >  1:234653/1‑149 (MQ=255)
                                                                gCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATg                                                                                   >  6:91458/1‑149 (MQ=255)
                                                                         cccGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAgg                                                                          <  5:260375/149‑1 (MQ=255)
                                                                         cccGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATaa                                                                                                                         <  6:254648/102‑1 (MQ=255)
                                                                         cccGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATaa                                                                                                                         >  5:254648/1‑102 (MQ=255)
                                                                          ccGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAggg                                                                         <  6:105993/149‑1 (MQ=255)
                                                                          ccGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAggg                                                                         >  3:211136/1‑149 (MQ=255)
                                                                           cGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAGAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGc                                                                        >  3:305227/1‑149 (MQ=255)
                                                                               aaCACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGc                                                                    <  2:27598/149‑1 (MQ=255)
                                                                                  acCCAGCGCAAGCTACTGCATGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAg                                                                 <  6:340069/149‑1 (MQ=255)
                                                                                  acCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAg                                                                 >  2:32478/1‑149 (MQ=255)
                                                                                    ccAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCa                                                               >  1:71223/1‑149 (MQ=255)
                                                                                        cgcAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTc                                                                                <  4:261068/128‑1 (MQ=255)
                                                                                        cgcAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACATGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTc                                                                                >  3:261068/1‑128 (MQ=255)
                                                                                              cTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTg                                                       >  5:74154/1‑147 (MQ=255)
                                                                                                   gCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGAt                                                                                    <  6:53233/113‑1 (MQ=255)
                                                                                                   gCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGAt                                                                                    >  5:53233/1‑113 (MQ=255)
                                                                                                    cAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCg                                                 <  2:71223/147‑1 (MQ=255)
                                                                                                      ggCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGa                                             <  4:211136/149‑1 (MQ=255)
                                                                                                        cACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCAACAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGAc                                            >  6:349163/1‑148 (MQ=255)
                                                                                                             gCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCAGCg                                       <  7:258795/148‑1 (MQ=255)
                                                                                                                      tCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCAGCGGCCAGAAGCa                             >  4:68950/1‑149 (MQ=255)
                                                                                                                         tGGACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCAGCGGCCAGAAGCATTc                          >  1:90315/1‑149 (MQ=255)
                                                                                                                          ggACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCAGCGGCCAGaa                                <  7:196152/142‑1 (MQ=255)
                                                                                                                          ggACGAACACCTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCAGCGGCCAGaa                                >  8:196152/1‑142 (MQ=255)
                                                                                                                               gAACACCTGCAACAGCGGCTGCACAAAGTTCAGCGTAGTCACAACATAAGAATGCAACAGGGCAGCGATACGCTGCGGGTCGCGATGTTCTAAAGGGCAGGCGAGCACCAGCTTGCAGACGGTGACAAg                                         >  3:67228/1‑129 (MQ=255)
                                                                                                                               gAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCAGCGGCCAGAAGCATTCCcaca                     >  2:243181/1‑149 (MQ=255)
                                                                                                                               caaCACCTGCAACAGCGGAGGCACAAAGTGAAGCGTAGTCACCCCATAAGCCTTCACGAGTGCAGCTATACGCTGCGTGTCGCGATGTTCACCAGTGCCTGCTAGAACCATCTTTCAGCCTTTGACCAg                                         <  4:67228/128‑1 (MQ=255)
                                                                                                                                 acacCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCg                                                 >  4:248847/1‑119 (MQ=255)
                                                                                                                                 acacCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCg                                                 <  3:248847/119‑1 (MQ=255)
                                                                                                                                     cTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACTCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGa                                                                  <  4:293117/98‑1 (MQ=255)
                                                                                                                                     cTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAg                                                    <  4:162469/112‑1 (MQ=255)
                                                                                                                                     cTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAg                                                    >  3:162469/1‑112 (MQ=255)
                                                                                                                                     cTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCAGCGGCCAGAAGCATTCCCACACCgaga               >  7:89214/1‑149 (MQ=255)
                                                                                                                                     cTGCAACAGCGGCGGCACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGTCCAGCGGCCAGAAGCATTCCCACACCgaga               <  1:32478/149‑1 (MQ=255)
                                                                                                                                                  ggTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCAGCGGCCAGAAGCATTCCCACACCGAGACGTCGAAACTGAt  <  6:74154/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGGTTGTGCAGTTGCACCTGCGGCAGCACCTGCAGCACGCGGTCACGCAGCGTGGCGGACAACGCTTCGCCACCCGAGAACACCCAGCGCAAGCTACTGCAGGCACTTGCCTGCGGCTCCTGGACGAACACCTGCAACAGCGGCGGTACAAAGTGCAGCGTAGTCACCCCATAAGCCTGCACCAGGGCAGCGATACGCTGCGGGTCGCGATGTTCACCAGGGCCGGCGAGCACCAGCTTGCAGCCGGTGACCAGCGGCCAGAAGCATTCCCACACCGAGACGTCGAAACTGAT  >  NC_002947/4826949‑4827241

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