Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 50,280 T→C F153S (TTT→TCT)  folA → type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927350,2800TC98.9% 343.2 / ‑4.2 90F153S (TTT→TCT) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (33/56);  total (33/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TTCCCGGATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTCCTA  >  NZ_CP009273/50195‑50366
                                                                                     |                                                                                       
ttCCCGGATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTgag                                                                                     >  1:858746/1‑90 (MQ=255)
     ggATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTgag                                                                                     <  2:967045/85‑1 (MQ=255)
     ggATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTgag                                                                                     >  1:967045/1‑85 (MQ=255)
     ggATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCt                                                                                >  1:864020/1‑90 (MQ=255)
     ggATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCt                                                                                >  1:1686552/1‑90 (MQ=255)
         tACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAg                                                                            <  1:2060487/90‑1 (MQ=255)
          atgAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgc                                                                           <  2:1301236/88‑1 (MQ=255)
           cGAGTCGGATGACTGGGAACCGGTATTCAGCGAATTCCACTATGCTGATGCGCAGAACTCTCACAGCTATTGCGCTGAGATTCTGGAgcg                                                                          <  2:937848/90‑1 (MQ=255)
           cGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgcg                                                                          <  2:3153668/90‑1 (MQ=255)
           cGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgcg                                                                          <  2:1559085/90‑1 (MQ=255)
           cGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgcg                                                                          >  2:1666586/1‑90 (MQ=255)
            gAGCCGGATGACTGGGAATCGGTATTCCGCGAATTCCACGATGATGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgcgg                                                                         >  1:2853930/1‑90 (MQ=255)
            gAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgcgg                                                                         >  1:1401474/1‑90 (MQ=255)
             aGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgcggc                                                                        <  1:1331124/90‑1 (MQ=255)
             aGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgcggc                                                                        <  1:1933019/90‑1 (MQ=255)
               ccGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGcggcgg                                                                      >  1:397649/1‑90 (MQ=255)
                cGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGt                                                                     <  2:49384/90‑1 (MQ=255)
                cGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGt                                                                     >  2:2511713/1‑90 (MQ=255)
                 ggATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTa                                                                    >  1:147593/1‑90 (MQ=255)
                     gACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt                                                                <  1:1816199/90‑1 (MQ=255)
                     gACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt                                                                <  2:1026802/90‑1 (MQ=255)
                     gACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt                                                                <  1:2990684/90‑1 (MQ=255)
                     gACTGGGAACCGGTATTCAGCGAATTCCACGGTGCTGATGCGCAGAACTCTCACAGCCATTGCTTTGAGAATCTGGAGCGGCGGTAAttt                                                                <  2:3065534/90‑1 (MQ=255)
                         gggAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGta                                                            <  2:1489776/90‑1 (MQ=255)
                         gggAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGGAAttttttt                                                            >  2:2064392/1‑87 (MQ=255)
                           gAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGtata                                                          <  2:1848988/90‑1 (MQ=255)
                           gAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGtata                                                          >  1:1726155/1‑90 (MQ=255)
                           gAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGtata                                                          <  1:2822701/90‑1 (MQ=255)
                           gAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGGAATTTTTtata                                                          >  2:447864/1‑90 (MQ=255)
                            aaTCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAg                                                         <  2:1530479/90‑1 (MQ=255)
                             aTCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGa                                                        <  1:3160170/90‑1 (MQ=255)
                                ggTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTTTTTAAAAtt                                                     >  2:2072240/1‑90 (MQ=255)
                                ggTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAAtt                                                     >  1:3063716/1‑90 (MQ=255)
                                   aTTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTAc                                                  <  2:3046236/90‑1 (MQ=255)
                                    ttCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACg                                                 >  2:2023116/1‑90 (MQ=255)
                                      cAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGc                                               <  1:2083330/90‑1 (MQ=255)
                                       aGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCt                                              <  1:1185013/90‑1 (MQ=255)
                                        gCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTTTATAGAATTTACGGCTa                                             >  1:839330/1‑90 (MQ=255)
                                        gCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTa                                             <  1:17522/90‑1 (MQ=255)
                                        gCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTa                                             >  1:2646442/1‑90 (MQ=255)
                                        gCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGGAATTTTTTTTTAAATTTTCGGCTa                                             >  2:2471335/1‑90 (MQ=255)
                                         cGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAg                                            <  2:1542583/90‑1 (MQ=255)
                                          gAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGGAATTTTGTTTAGAATTTTCGGCTAgg                                           >  1:3074261/1‑89 (MQ=255)
                                             ttCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTTTATAAAATTTACGGCTAGCGcc                                       >  1:1198844/1‑90 (MQ=255)
                                              tCCACGATGCTGATGCGCAGAACTCTCACAGCTAATGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCg                                      <  2:3243104/90‑1 (MQ=255)
                                                 aCGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGAt                                   <  1:2023116/90‑1 (MQ=255)
                                                     tGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAACGCCGGATGCGa                               >  1:1881439/1‑90 (MQ=255)
                                                         gATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAATTTGTATAAAATTTACGGCTAAGGCCCGAGGCGACGcc                           >  1:1628051/1‑90 (MQ=255)
                                                             cgcAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTGTAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTc                       <  2:2847490/90‑1 (MQ=255)
                                                             cgcAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTc                       <  2:2345050/90‑1 (MQ=255)
                                                             cgcAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTc                       <  1:891202/90‑1 (MQ=255)
                                                             cgcAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTc                       <  2:3063716/90‑1 (MQ=255)
                                                             cgcAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTc                       <  2:1899475/90‑1 (MQ=255)
                                                              gcAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTTTATAGAATTTACGGCTAACGCCGGATGCGACGCCGGGcg                      >  1:601993/1‑90 (MQ=255)
                                                              gcAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGGAATTTTTTTTAAAAATTACAGCTAACGCCCGGGTGGACCCCGGTcg                      >  2:2139268/1‑90 (MQ=255)
                                                                 gAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGt                   <  1:427318/90‑1 (MQ=255)
                                                                 gAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGGAATTTTTTATAGAATTTTCCGCTAGCCCCCGGTGCGGCGCCCGGcgcgg                   >  2:2641088/1‑89 (MQ=255)
                                                                 gAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGGAATTTTGTATAGAATTTTCGGCCAGCGCCGGAAGCGACACCGGTCGCGt                   >  1:2350767/1‑90 (MQ=255)
                                                                  aaCTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTc                  <  1:733922/90‑1 (MQ=255)
                                                                         acaGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCg           <  1:2797773/90‑1 (MQ=255)
                                                                           aGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTTTATAGAATTTACGGCTAGCGCCGGGTGCGACCCCCGGCGCGTCTTTTCCCGc         >  1:1268298/1‑90 (MQ=255)
                                                                           aGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGc         <  1:305226/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCTGGcc        <  2:86121/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  1:2484922/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  1:1088959/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  1:2713582/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  1:26919/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  2:3060799/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  2:864020/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  1:1596728/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  2:3336082/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  2:1054756/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  2:799640/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  2:591540/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  1:2236557/90‑1 (MQ=255)
                                                                            gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc        <  2:1881439/90‑1 (MQ=255)
                                                                             cTATTGCTCTGAGATTCTGGAGCGGCGGGAATTTTGTATTGAATTTACGGCTAACGGCGGGGGGGGGGCCGGTCGCGTCTTTTCCGGGCt       >  2:8983/1‑90 (MQ=255)
                                                                              tATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGt                                                             >  2:1577901/1‑36 (MQ=38)
                                                                              tATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGt                                                             <  1:1577901/36‑1 (MQ=38)
                                                                              tATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGGCGCGTCTTTTTCCGCCtt      >  1:1074407/1‑90 (MQ=255)
                                                                               aTTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAACGCCCGGTGGGGCGCCCGGCGCGTCTTTTCCGGCCTTc     >  1:1362369/1‑90 (MQ=255)
                                                                                ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTcc    <  2:642621/90‑1 (MQ=255)
                                                                                ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTcc    <  2:1229284/90‑1 (MQ=255)
                                                                                ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTcc    <  1:2048815/90‑1 (MQ=255)
                                                                                ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTAcc    <  2:2064376/90‑1 (MQ=255)
                                                                                ttGCTCTGAAATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTcc    <  1:3175833/90‑1 (MQ=255)
                                                                                  gCTCTGAGATTCTGGAGCGGCGGTAATTTTTTATAGAATTTACAGCTAGCGCCGGCTGCGTCGCCGGCCTCGGCTTATTCGGCCTTcccc  >  1:1570707/1‑88 (MQ=255)
                                                                                  gCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTCCta  <  2:2165751/90‑1 (MQ=255)
                                                                                  gCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTCCta  <  1:8983/90‑1 (MQ=255)
                                                                                  gCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTCCta  <  1:2595298/90‑1 (MQ=255)
                                                                                     |                                                                                       
TTCCCGGATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTCCTA  >  NZ_CP009273/50195‑50366

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